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The Hubbard model extended by either nearest-neighbour Coulomb correlation and/or nearest 
neighbour Heisenberg exchange is solved analytically for a triangle and tetrahedron. All eigen- 
values and eigenvectors are given as functions of the model parameters in a closed form. The 
groundstate crossings and degeneracies are discussed both for the canonical and grand-canonical 
energy levels. The grand canonical potential T, h) and the electron occupation N (n, T, h) 
of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the 
three interaction constants. In the pure Hubbard model we found various steps in N(fj., T = 0, h) 
higher than one. It is shown that the various degeneracies of the grand-canonical energy levels 
are partially lifted by an antiferromagnetic exchange interaction, whereas a moderate ferromag- 
netic exchange modifies only slightly the results of the pure Hubbard model. A repulsive nn 
Coulomb correlation lifts these degeneracies completely. The relation of the cluster gas results 
to extended systems is discussed. 



1 Introduction 

Besides the fact that the analytical solution of non-trivial models of strong correlation is a task of 
high pedagogical value, the renewed interest in the analytical solution stems mainly from cluster 
methods, which were successfully developed within the context of strong electron correlation 
during the last decade [1, 2, 3] and applied to problems of high-T c superconductivity and from 
modelling of electron transport through quantum dots, see e.g. [4] and the references therein. 
Both topics have in common that a detailed knowledge of the cluster physics is a key ingredient. 
The by far most employed model of strong correlated electrons is the Hubbard model. Despite its 
simplicity, the only analytical solutions available up to now are the groundstate for n-site rings 
by help of the Bethe ansatz or transfer matrix method [5, 6] and the complete solution on small 
clusters, i.e. the triangle, the tetrahedron [7] and the square [8]. For models containing five and 
more sites only approximate or numerical solutions exist. A further topic of increasing interest 
is the interplay between frustration and strong electron correlation - a topic which acquired a lot 
of interest due to the discovery of a superconducting phase in the hydrate of Na x Co02- For a 
theoretical survey to that topic I refer to [9] and the literature therein. Furthermore, the latter 
paper demonstrates explicitly the use one can make from the exact knowledge of small clusters 
in understanding the results of more elaborate many body methods. Although I have up to now 
not dealt personally with attractive Hubbard models this point seems to be also of some interest 
[10, 11]. Thus researchers in this field may probably also benefit from our work, since the results are 
valid for arbitrary values of the model parameters. The present paper is very related to our former 
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analytical solution of the Hubbard model on a square and the related cluster gas [8] , thus we will 
extensively refer to it. One of the main points made in [8] was the discovery of the degeneration 
of cluster states with different particle numbers and spin, which has surely a great influence of the 
coupling of spin- and charge degrees of freedom, and indeed, the degeneracy of the groundstates 
for the 4-site cluster occupied either with two or with four electrons, is responsible for the charge 
separation and inhomogeneitics in the pseudogap phase of the cuprates, whereas their bosonic 
character accounts for the superconducting phase [12, 13]. Also the d-wave symmetry of the gap 
seems to be a direct consequence of the symmetry of the 4-site cluster groundstates [14]. Usually 
the addition of further terms to the pure Hubbard Hamiltonian destroys this degeneracy and makes 
the system unsolvable due to a reduced number of global symmetries, with the consequence that 
we have a reduced number of analytical solutions and exact numerical results [15]. In the present 
paper we study the influence of additional nearest-neighbour (nn) Coulomb interaction, exchange 
term in the case of both the triangle and the tetrahedron or next-nearest neighbour (nnn) hopping 
in the case of the square respectively, which are added to the standard Hubbard model. The latter 
is used in the notation of [8] . A common feature of the considered models is frustration, which is 
present due to geometry in the triangle and tetrahedron, even without the additional terms, and in 
the square due to the nnn-hopping, what becomes evident since for t 1 = t the model is equivalent 
to the pure Hubbard model on the tetrahedron. For the triangular and tetrahedral geometry the 
spatial and spin symmetries are enough to get a complete solution. For the square geometry any of 
the additional terms breaks the global pseudospin symmetry, what prevents closed form solutions 
for all eigenvalues and eigenstates. Nevertheless, there exist non-trivial parameter sets, where 
analytical solutions exist at least for the canonical partition sums. In the following we abstain 
to discuss the square geometry since the most important solvable case t = t' is contained in our 
results for the tetrahedron. In the next section we shall present the extensions of the models 
and its closed form solutions, which we shifted to the appendices due to its unusual large page 
volume. In the third section the low lying levels and its parameter dependencies are discussed, 
where emphasis is given to their dependence on the additional interaction (hopping) parameters. 
A lot of degeneracies and ground state level crossings are presented, which may give rise to a 
likely rich physics in the "cluster gas" or extended systems built from these clusters as it was 
observed in the relation of the 4-site Hubbard model on a square to the square lattice. In a 
concluding section I will discuss the results. At least I want to beg for the readers pardon for the 
unusual volume of this article, which is due to the closed form expressions of the eigenstates of the 
models. The reason is that in response to Ref. [8] a lively interest arose in the closed form of the 
eigenstates, which I have published on my website only, whereas the thermodynamical results of 
the model contained in the Ref. [8] obviously did not find as much interest. Unfortunately I can 
not guarantee for the future existence of this website, therefore I decided to include the eigenstates 
as appendix to provide a hard copy to the interested reader this time. Whereas preliminary results 
on thermodynamics and magnetics were published in Ref. [16], more detailed data on these topics 
will be published in a forthcoming paper. 

2 The Hamiltonian 

In the following we consider the model 



H = Hh + He + H / 



(1) 



with Hh being the pure Hubbard model. 




(2) 
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He is the nn-Coulomb repulsion 
U C = W n ^ ■ 

Furthermore a nn-Hciscnberg exchange term Hj 
Hj = J SiSj 



(3) 
(4) 



was added. Here cf a and Ci a are the creation and destruction operators of electrons at site i 
with spin a, nj CT = c+c i(J and = J2ia n i<r- &i indicates the local spin operator at site i. 
For a detailed physical reasoning of the additional terms in the Hamiltonian see e.g. [17]. Our 
Hamiltonian is an extension of the Hamiltonian used e.g. in [18] (Although up to now an erata 
was not published, we think that the a in the hopping term of cq. 1 of that paper is simply a 
misprint), since we treat the nearest-neighbour correlation W and the nearest-neighbour exchange 
constant J as independent parameters. This is done mainly to cover pure Heisenberg type models 
by setting W/t — 0, J/t ^ 0. The chemical potential /x and the magnetic field h in z-dircction are 
introduced to take into account the effects of doping and applying external magnetic fields. Please 
note, that the signs in front of the hopping and exchange term are positive, thus one has to be 
carefully while comparing with papers where other conventions are used. Regarding the method 
of calculating the eigenvalues and eigenvectors we refer to [7, 8]. The analytical solution, i.e. the 
closed form expressions for the eigenvalues and eigenvectors in dependence on the parameters U, 
J, and W are shifted to the appendices, where they are numbered. These numbers are occasionally 
used for referencing in the following sections. Since I did not found any analytical solution of that 
model, it was hard to prove the results. For W/t = J/t = we recover the results of [7], i.e. the 
characteristical polynomials and characteristics of the eigenstates, which themself are not given 
explicitely. For the triangle we recover, again setting W/t = J/t = 0, the analytical result for the 
groundstate energy and the wave function with 4 electrons given in [9]. Furthermore for small 
electron occupation or some limiting cases the results could be proved by a pencil calculation. 
Nevertheless, the main confidence stems from our experience with the Hubbard model on the 
square. A very difficult point was also the transcription of the analytical computer results into a 
human readable form. Since this had to be done automatically also, additional progamming was 
necessary, and finally I had to manhandle some of the expression to format it for publication, a 
further chance to introduce mistakes. 



3 The triangle and the related cluster gas 

3.1 The spectrum of the triangle 

In Fig. 1 we study the dependence of the canonical eigenvalues on the parameters V, J', and W' 
separately. The new parameters are defined according to 

u ' = uTW\ < B > 
r = TW\ (6) 

W = ^— (7) 

t+\w\ y) 

in order to map the parameter space from the interval (— oo,+oo) to (—1,1). The mapping 
of the on-site correlation was chosen in a different way, since we are mainly interested in the 
strong interacting case, and thus, roughly speaking, we achieved that the weak interaction region 
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is limited to the interval (—0.5, 0.5) and the strong correlated case outside. To avoid infinite 
energies we map the energy in the same way 

E 

E ' = — (8) 
t+\E\ y ' 

A multitude of level crossings in dependence of the different interaction parameters is found. 
Although level crossings of higher states may be important for spectroscopy, we focus to the 
level crossings of the canonical groundstates, due to its relation to quantum phase transitions in 
extended systems. First of all there is only one groundstate level crossing in dependence of U if 
J and W is set to zero. This occurs if the cluster is occupied with two electrons and U changes 
sign. In order to study the dependence of the groundstate on the other parameters we fix U/t to 
four. The dependence on W is also very simple, since there is one levelcrossing for two electrons 
and W = U only. Otherwise, the dependence on the exchange parameter is very interesting. 
In Fig. 2 we show all the groundstate levelcrossings. The insets show the J dependence of the 
complete canonical spectra. The levels are denoted in the legend by their quantum numbers, 
if the cluster is occupied with two electrons we find an U independent crossing from a twofold 
degenerate groundstate with spin one and spin projection one (no. 8, 20) to the non-degenerate 
eigenstate (no. 13) with spin one and spin projection zero. These states belong both to the 
one-dimensional irreducible representation r 2 . At the second crossing point, a crossover to two 
eigenstates of spin zero (no. 14, 17) occurs. The symmetry changes to T 3 resulting again in 
a twofold degenerate groundstate, due to the two-dimensionality of the irred. representation. 
Whereas the first transition happens always at J c /t = 0, the second one is dependent on the 
on-site correlation according to 

J - = — (9) 
t 2t + U • { ' 

For three electrons there is only one ground state level crossing, changing from fourfold degenerate 
level (no. 27, 30, 36, 39) with spin one half and irred. representation T3 to the two-fold degenerate 
states(no. 23, 46) with spin three half and irred. representation T 2 . The crossing point J c is 
implicitely given as function on U by 

= 11 J -4 £7 + 4^27*2 + J 2 + 2 JU + U 2 cos ( I arccos( ( J + U ) 5-)) (10) 

V 3 (27i 2 + J 2 + 2JU + U 2 y J 

Fig. 3 shows the solution for the primed values. A crossover from a fourfold degenerated ground- 
state (no. 44, 45, 56, 57) of symmetry T3 with spin one and spin projection ±1 to a singlet 
(no. 46) with Ti and spin zero takes place for negative values of J (remember: J < favourites 
parallel spins in our notation) if the electron occupation is four. The U dependence of the critical 
exchange parameter is given by 

J c = \(2t + 3U- ^112t 2 + 24Ut + 9U 2 ). (11) 



3.2 The cluster gas 

In the following the most interesting aspects of an ensemble of triangular clusters (clustergas) will 
be discussed. In the left panel of Fig. 4 we show the dependence of the electron occupation on the 
chemical potential for the pure model for different values of the on-site electron correlation. It is 
obvious, that the 16-fold degeneracy of states of one, two, three, and four electron states is lifted 
by the on-site electron correlation. But the more interesting feature is the jump of the occupation 
number from to 2. The reason is, that a N — 1 state is never the (grand-canonical) ground 
state. This feature is persistent even for arbitrary large values of U, since the N — 2 groundstate 
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avoids correlation. This situation is analog to the four-site model, where states with two, three, 
and four electrons are degenerated for certain values of the chemical potential. It was conjectured, 
that this may be one of the main issues for both the superconduction and the stripe structure 
if the high-Tp-cuprates are modeled by a square lattice Hubbard model. Since the triangular 
cluster gas is related in the same sense to the triangular lattice, it seems reasonable to conjecture 
that we have analogous peculiarities in systems, which are modeled by a Hubbard model on a 
triangular lattice. The recently reported superconductivity in Na x Co02 may be a hint in this 
direction. Although we are mainly interested in the repulsive model we show also the curves for 
the attractive case (U < 0) for seek of completeness. The right panel depicts the groundstate 
"phase diagram", i.e. the eigenvalues of the groundstates in dependence on the chemical potential 
and the electron correlation. In order to cover the whole parameter range we used again eq. (5) 
and introduced /x' according to 

p! = T^r- r (12) 
t+ \n\ 

The meaning of the colours is explained by help of the palette shown in Fig. 5, whereas the left 
picture shows the groundstate of the complete fi'-h' plane, with 

b! = — \- . (13) 
t+\h\ v ' 

and U/t fixed to four. This U value was chosen since it is the region, where in extendeded systems 
the cross-over from weak correlation to strong correlation happens, what means that perturbation 
methods can not be applied . The influence of an ferromagnetic (J < 0) and antifcrromagnetic 
exchange (J > 0) interaction is studied in Fig. 6. This interaction is not able to lift the mentioned 
degeneracy of the groundstates. And the same holds for an attractive nearest-neighbour Coulomb 
interaction, whereas the more realistic case of repulsive interaction does, as may be seen in Fig. 
7. Regarding the grand canonical spectrum of the triangle the most interesting physics results for 
(U = —1/2. In the left part of Fig. 8 we depicted the lowest part of the grand-canonical spectrum. 
At this point, an empty state is degenerated with the four single-electron states belonging to 
irreducible representation T3 and furthermore with a two-electron triplet of symmetry T2, if the 
nn-interactions are neglected. A small external magnetic field lifts this degeneracy partially. If 
the exchange interaction is taken into account this degeneracy is also lifted partially as may be 
seen from the right panel of Fig. 8, where we depicted the grand-canonical level system for both 
a ferromagnetic and an antifcrromagnetic exchange integral. 



4 Results for the Tetrahedron and the related cluster gas 

4.1 The spectrum of the tetrahedron 

In Fig. 9 we depict the dependence of the canonical eigenvalues on the parameter U' . The degener- 
acy can be seen from the legends, since the irreducible representations Ti, T 2 are one-dimensional, 
T3, is two dimensional, and T^, are three-dimensional. Furthermore these degeneracies have to 
be multiplied by the normal spin degeneration. For the pure Hubbard model there is a ground- 
state crossing at U = for an electron occupation of 2, 3 and 4 but more interesting is its absence 
for N = 5 where the groundstate T 4 , S = 1/2 is sixfold degenerated and for N — 6 where we 
find a nondegenerated Ti, S = groundstate for the whole parameter range. The groundstate in 
dependence on U and J is plotted for the cases, where level crossings in the ground state appear. 
Note that we used again primed parameter values to cover the whole plane. The meaning of 
the colours is explained in the palette shown at the down-right corner. The dependence on the 
exchange-parameter J' is shown in Fig. 10. Due to the exchange interaction the degeneration is 
lifted and the degeneration factor reduces to two, besides the r-multiplicity. Whereas the ground- 
state crossing for the N = 2 system was rather expected, since it is coupled to the sign change 
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of J, we find a more complex behaviour for three and four electrons on the cluster. As may be 
seen from Fig. 11 we find a change from a sixfold degenerated groundstate with Ts, 5 = 1/2 to 
a twofold degenerated groundstate with T 2 , 5 — 1/2 at a positive value of J for reasonabel U 
values. 

10 t + 3 U - VlOO t 2 + 48 t U + 9 U 2 
Jc ~ 3 

what amounts for U = At to J c /t = (22 — 2 \/109) /3 = 0.373129. Another crossover happens when 
J changes the sign, since for negative values a 5 = 3/2 state with S z — 3/2 is groundstate. Since 
in the nature the high degeneration of the T 5 groundstate will be lifted by a Jahn- Teller distortion, 
one has to expect the system to undergo a structural transition at positive values of J followed 
by a magnetic one when J changes sign. A similar situation happens, if the cluster is occupied 
with four electrons, see Fig. 12. Now, we have transitions for negative values of J only the first 
one, which is structural and magnetic, happens for J c ,i(U = At) = — 0.19106 1 and the second one, 
which is purely magnetic, for J C ,2(U = At) = —0.575885 t. The dependence of J c from the Hubbard 
parameter U one gets from equating the eigenvalues No. 115 with No. 102 and No. 102 with 
No. 94 respectively (for the eigenvalues see appendix B). These equations contain the parameter 
J in an essentially non-algebraic way and have to be solved numerically. The degeneracy of the 
groundstate steps from four to six in the first crossover and remains six in the second one. For 5 
electrons on the cluster we find only one crossing point J C (U = At) = —3. 02201 1 from equating 
the eigenvalues No. 164 with No 174. If the cluster is occupied with six electrons a crossover from 
a non-degenerated non-magnetic state (T\, 5 = 0) to a sixfold degenerated groundstate happens 
at 

At + 3U - V208t 2 + A8tU + 9U 2 
Jc ~ 3 

what yields -2.44127 1 if U = At is used. 
In Fig. 13 we plotted the canonical levels in dependence on the nearest neighbour Coulomb 
interaction W scaled to the primed value W for J = 0. Again, we have chosen U = At. Here we 
find one levelcrossing for W c = U/2 in each of the cases of 2, 3, and 4 electrons on the cluster and 
none for the other occupation numbers. At the first glance this seems to be not very interesting, 
but if one looks closer, one recognises a transition from a singlet state with Ti, 5 = to an 
eleven-fold (!) degenerated groundstate (nine with T 5 , 5=1 and two with r 3 , 5 — 0) for N = 2 
and W < U/2. For N = A the transition is to a doublet with T 3 and 5 = 0. For three electrons 
we have for W > U/2 a Ti 5=1/2 doublet as groundstate and below the transition a T 3 spin 
quartet. 

4.2 The tetrahedral cluster gas 

The most interesting part of this paper is surely the ensemble of tetrahedral clusters (clustergas) , 
which are coupled to a common electron reservoir. In Fig 14 we study the electron occupation in 
dependence of the the chemical potential for different values of the Hubbard correlation energy. 
In the trivial case of U = we have a two step function only, due to the huge degeneration of the 
grand-canonical groundstate. Of course this degeneration is partially lifted if the on-site electron 
correlation is switched on. Nevertheless there remains a huge degeneration of grand-canonical 
groundstates with 0, 1, 2, and 3 electrons at fj, = —t. Indeed the empty state is degenerated with 
the six N = 1 T 4 states, further with the two N = 2 5 = T 3 and nine N = 2 5 = 1 T 5 states 
and finally with the four N = 3 5 = 3/2 T 2 states. Of course this behaviour is due to the fact, 
that for such low numbers of electrons uncorrelated many-particle states exist without correlation. 
This situation is drawn in detail in Fig. 16. For other values of the correlation energy U this 
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degeneration happens at 

3U VmWTsWTTp Vl6i 2 + C/ 2 (l , 12 yftpu A ,„„, 

u = —t H - - 1 = cos - arccosl T ) (16) 

4 4 V3 \3 \wt* + urf J ) 

The situation is completely different for fi/t = 2.4056, where we find a degeneration of the two 
N = 4 5 = T 3 states with the single TV = 6 S = Ti state, as shown in Fig. 17. A N = 5 state is 
never a grand-canonical groundstate. The situation resembles closely the Hubbard square, where 
we also found a step from N — 2 to N = 4 and from N = 4 to N — 6 [8]. In Fig. 18 we depicted 
the mentioned grand-canonical groundstate energies over the chemical potential in very detail. 
Now, taking an exchange interaction into account the situation is changed according to Fig. 19. 
Whereas a ferromagnetic exchange does not change the qualitative picture, an antifcrromagnetic 
does. For that case the system avoids an odd number of electrons to stay with total spin zero. 
For an extended system this means that for an odd mean electron occupation only clusters with 
even electron occupation will exist in the groundstate. Of course \J/t\ = 1 is usually a to large 
value. We have chosen it to make the plots clearer. Nevertheless the picture is not qualitatively 
changed if smaller values are considered as one can see from the right panel in Fig. 19. As one 
can see from Fig. 20 a repulsive nearest-neighbour interaction lifts the degeneracy of levels with 
different electron occupation completely Otherwise, a large attractive interaction would result in 
a one-step function. 



5 Discussion 

Since my original interest for these clusters stems from extended systems I will focus the dis- 
cussion to the consequences of the cluster behaviour for crystals. Although the eigensystem and 
the shown figures cover the whole parameter space from — oo to +oo for all model parameters I 
will not discuss the cases with attractive Coulomb interaction, since I believe they are of minor 
interest with respect to real electron systems only. Nevertheless papers dealing with attractive 
Hubbard models are published from time to time, thus the presented results may be interesting 
to that community also. 

The two (or three if the square with diagonal hopping t' = t is taken as a separate model) clus- 
ters are considered usually as "frustrated" systems. The reason is mainly that for large on-site 
coulomb interaction the Hubbard model may be mapped onto an anti-ferromagnetic Hcisenberg 
model with the consequence that the localised magnetic moments can not be arranged in such a 
configuration that all localised spins are antiparallel. The reader may wonder that the conception 
of frustration was never touched in the paper. That is due to the fact, that thinking in localised 
spins means thinking in one-particle states, which do not play any role here, since the basis states, 
e.g. |2, u, 0, d), are many particle states from the beginning. Of course, if for instance the electron 
or spin density as explicite functions of the position are wanted, one has to specify the action of 
the creation operators, what is done usually by means of determinants or one-particle-orbitals in 
clusters or Wannier functions in crystals. Nevertheless, in my opinion the consequences of frustra- 
tion are reflected in some sense by the some times large degeneracy of the ground states and the 
degeneration of ground states with different occupation numbers. The "frustration" increases the 
groundstate energy for some occupation numbers in such a manner that adjacent "non- frustrated" 
ground states have lower grand canonical energies, i.e. exhibit higher (or equal) weights in the 
grand-canonical partition sum. This feature, first observed in the Hubbard model on the square, 
is much more pronounced in the Hubbard model on the triangle and tetrahedron. In a periodical 
arrangement of such clusters with a weak inter-cluster hopping, these degeneration will result 
in patterns built up from clusters with different electron occupation, alike the stripe structures 
observed in the copper-oxide planes. But in the case of the triangular or tetrahedral array these 
patterns will be much more richer, since we have for the latter e.g. steps from zero to three (for 
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the pure Hubbard model at T = 0) in the underdoped region or from four to six in the overdoped 
region. Indeed, an material which can be described by an array of tetrahedra with dominating 
on-site correlation will show very intersting features if it is doped to five electrons per tetrahedron. 
If there is a small intercluster hopping, the three groundstates (see Fig. 17) are broadened to small 
bands with the bandwidth of the order of the (small!) intercluster hopping. Since the N = 4 and 
N = 6 bands will be magnetically inactive for instance a small magnetic field will split the N = 5 
groundstatc shifting the lower one below the TV = 4 or N = 6 groundstates. An equal effect could 
be reached by a Jahn- Teller scenario, but since a distortion of the tetrahedral symmetry will also 
split the two fold degeneration of the N — 4 groundstate, whereas the non-degenerated ground 
state with six electrons will remain in position, it is not clear, which state will be the groundstate, 
but this could be decided easily by an experiment determining the ground-state spin density, or by 
including the elastic energy into the Hamiltonian and treating it perturbatively. Before one deals 
with such subtle features, we had to answer the question, to wether degree a nearest-neighbour 
exchange or coulomb correlation changes the spectra of the pure Hubbard model. A repulsive 
nearest neighbour coulomb interaction destroys the degeneration of groundstates with different 
electron occupation. But if it is smaller or in the order of the width of the bands formed from the 
cluster states, if an array is considered, the subtle interplay of groundstates with different occu- 
pation numbers survives probably. But very interesting is the fact that an exchange interaction 
will not necessarily lift these degeneracies. For example a ferromagnetic exchange does not change 
the step curve of the pure Hubbard model qualitatively, whereas a relative big antiferromagnetic 
exchange results in a step function avoiding odd number occupation. Furthermore, a detailed 
view of the spectra in dependence of the parameters offers several ground state level crossings in 
dependence on the parameters U, J and W. For extended systems these level crossings will result 
in "quantum phase transitions". If the chemical potential is far from the "degeneration points" 
a weak-coupled cluster array will contain clusters with a well fixed integer number of electrons. 
The "quantum phase diagrams" of such a system will be very alike to the related "groundstate 
diagram" of the cluster, calculated in the same way as it was done for Figs. 14, 19 and 20. 

To relate the results of the present paper to arrays of triangular or tetragonal clusters it may 
be helpful to recall what we learned from the analytical solution of the square-cluster gas ap- 
plied to the Hubbard model on a square lattice, since for that case numerical cluster methods 
provided additional insight [1, 2, 20, 21]. Due to the inter-cluster hopping the three degenerated 
groundstates with N=2,3,4 gain a bandwidth and will be probably shifted slightly in position, 
nevertheless an overlap in some region remains. If the chemical potential is within this region 
the consequence is that we have degenerated non-magnetic and magnetic band states stemming 
from the N=2,4 and the N=3 cluster states. Since I see no reason, why the effective mass of these 
states should be the same, it will be not astonishing that the two charged (and bosonic!) states 
(remember: for a dotation of N=3 to the cuprates neutrality is restored by the ions not considered 
in the square Hubbard model) and the uncharged magnetic (and fermionic!) state will propagate 
with different velocity (of course this is the spinon-holon picture extended to the cluster states) 
resulting in the so-called spin-charge separation [1]. This feature is to expect within triangular- 
or tetragonal cluster arrays also, since the main reason, i.e the degeneration of spinless charged 
and uncharged magnetic grand canonical groundstates, exists as well. 

Finally I want to speculate about the relation of my results to (high-T c -) superconductivity in 
arrays of triangular or tetragonal clusters. It was shown [12, 13, 14], that the clustergas already 
accounts for a lot of the features of the phase diagram of the high-T c cuprates. Both the position 
of the superconducting domes (for electron and hole dotation) for low temperature and the stripe 
phase at higher temperatures are roughly spoken in the vicinity of the degeneration points of 
the related cluster gas, i.e. for 3 or 5 electrons per 4-site cluster. The slight asymmetry in the 
phase diagram for hole and electron doping within the cuprates hints to a broken electron-hole 
symmetry, a feature which can not be modeled by help of the pure Hubbard model. Otherwise 
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the models considered here lack this symmetry at all. The most interesting physics for electron 
dotation will be e.g. in an array of tetrahedrons if five electrons per cluster are present, whereas 
for hole dotation it will be a much broader region between one and three electrons per cluster. 
Regarding superconductivity I believe that a good approximation for the superconducting ground- 
state is a covering of the two dimensional plane by a mixture of degenerated cluster states. For 
instance for an array of squares coupled to an electron bath a mean occupation of five electrons 
per cluster would result in a ground state containing x % N = 4, y % N = 6, and 1 — x — y 
% N=5 cluster states. Since these three cluster ground states are degenerated, interchanging of 
an N = 4 cluster and an adjacent N = 5 or TV = 6 cluster will need no additional energy. Of 
course this interchange is forbidden if the clusters are isolated. But introducing a small hopping, 
as it is done in the checkerboard model [14] will then, besides selecting one state out of the vast 
number of degenerated coverings, allow for that process, thereby propagating one electron (or two 
if an N=4/N=6 interchange is also allowed) from one cluster to another. In the square lattice 
such an interchange is only possible in (1,0) and (0,1) direction, but forbidden over the diagonal 
(as may be seen easily for the checkerboard covering). If one adopts this simple picture the d 
symmetry of the charge propagation turns out to be a simple geometrical effect. Also the fact 
that the results calculated within the numerical cluster methods by help of a 4-site cluster do 
not show any qualitative new features if the cluster size is increased hints to the central role of 
the four-site square cluster [1, 2]. If an array of triangles or tetrahedrons is considered the same 
simple mechanism will take place if the system is doped to a chemical potential in the vicinity of 
the relevant degeneration points, i.e. where N(fi) steps at least by two, of course different geo- 
metrical constraints will be effective. Apart from the introductorily mentioned but not discussed 
topics, where the presented cluster solution may be useful, there is, due the fast development of 
nanotechnology, also some chance that two-dimensional arrays of triangular or tetragonal clusters 
with various geometries will be available and utilised for experiments in the near future. 

Acknowledgements I am grateful to K. Kikoin and D. Efremov for discussion on Kondo tunneling 
through small clusters. 
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A Eigensystem of the triangle 



In the following tables the first column gives the number of the state for referencing in the text. The 
second column labels the eigenstate by the eigenvalues of the U-independcnt symmetry operators, 
i.e. the electron occupation number N e , the spin-projection S z in z-direction m s , the eigenvalues 
s(s + 1) of S 2 and the spatial symmetry is indicated by T^j, where the first index labels the 
irreducible representation of the tetrahedral group and the second numbers the partner. The 
notation is based on Rcf. [19]. The third column gives the enery eigenvalues in abbreviated 
form. The abbreviations are listed subsequently to the tables. In the last column a numerical 
value is given for example. For comparison wc have chosen the same parameters as in Ref. [8] 
corresponding to a pure Hubbard-model with U = 5t, W = J = Ot. If the grand-canonical energy 
levels in an applied magnetic field are needed one has to substract \i N e + hm s . 



A.l The eigenvalues of the triangle 



Eigenkets and eigenvalues for N e =0 and m s = 0. 



No 


Eigenstate 


Energy 


Example 


1 


10,0,0,11) 





0. 


Eigenkets and eigenvalues for N e =l and m s = 


No 


Eigenstate 


Energy 


Example 


2 
3 
4 


|lj |) r 3,l) 
|l,-2' 4' r 3,2) 


It 
-t 

-t 


2. 

-1. 
-1. 


Eigenkets and eigenvalues for N e =l and m s = \. 


No 


Eigenstate 


Energy 


Example 


5 
6 
7 


!> J> |' r 3,l) 
|1, 3) 1^3,2) 


It 
-t 

-t 


2. 

-1. 
-1. 


Eigenkets and eigenvalues for N e =2 and m s = —1. 


No 


Eigenstate 


Energy 


Example 


8 
9 

10 


|2,-i,2,r 2 ) 

|2,-i,2,r 3 ,i) 

|2,-i,2,r 3;2 ) 


^ -2t + W 

+ + t + w 
i + t + w 


-2. 
1. 
1. 
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Eigenkets and eigenvalues for N e =2 and m s = 0. 



No 


Eigenstate 


Energy 


Example 


11 


|2,o,o,ri) 


— y/A B — J+2 t+U+W 
2 


0. 


12 


|2,o,o,ri) 


^M-J+2t+U+W 


6. 


13 


|2,0,2,r 2 ) 


-2t + W 


-2. 


14 


|2,o,o,r 3 ,i) 


-\/M-J-t+U+W 
2 


-1.37228 


15 


|2,o,o,r 3 ,i) 


y/M-J-t+U+W 


4.37228 


16 


|2,0,2,r 3il ) 


t + W 


1. 


17 


|2,o,o,r 3 , 2 ) 


-s/~A~[-J-t+U+W 
2 


-1.37228 


18 


|2,0,0,I\ 2 ) 


</A~^-J-t+U+W 


4.37228 


19 


2,0,2,I\ 2 ) 


t + W 


1. 



Eigenkets and eigenvalues for N e =2 and m s = 1. 



No 


Eigenstate 


Energy 


Example 


20 
21 
22 


|2,i,2,r 2 ) 

|2,i,2,r 3 ,i) 

|2,i,2,r 3 , 2 ) 


±-2t + W 

+ + t + w 
i + t + w 


-2. 
1. 
1. 


Eigenkets and eigenvalues for N e =3 and m s = — §. 


No 


Eigenstate 


Energy 


Example 


23 


|3, -f , X' r 2) 


3 (J+2 W) 
2 


0. 


Eigenkets and eigenvalues for N e =3 and m s = — |. 


No 


Eigenstate 


Energy 


Example 


24 
25 
26 
27 
28 
29 
30 
31 
32 


|3, — 2 ' \ ' ^i) 
3, — o; 41 ^2} 
3> _ 2' T> 

3, — 4; r 3 ,i) 

3, 41 T3,l) 

3, -31 41 T3,i) 
3, —31 41 T3, 2 ) 

3, —91 Z' ^3,2) 
3, - 2 l 4' r 3,2) 


£7 + 214^ 
U + 2W 

J+6W 
2 

-J+2 C/+7 W-2 -JM cos(0i) 

I 

I 
3 

-J+2 U+7 W-2 cos(Oi) 

I 

I 
3 


4. 
4. 
0. 

-1.27492 
6.27492 
3. 

-1.27492 
6.27492 
3. 



14 



Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =3 and m s = ^ 


No 


Eigenstate 


Energy 


Example 


33 
34 
35 
36 
37 
38 
39 
40 
41 


3, 47 ^2) 

3, ^, -3-, r 2 ) 

3, \, 4,r 3j i) 

3, \ , 4^3,1} 

3, ^, f , r 3: i) 
3, \, 4,r 3)2 ) 
3, |, 4,r 3 , 2 ) 
|3, 5, f ,r 3j2 ) 


U + 2W 
U + 2W 

J+6W 
2 

-J+2 U+7 W-2 cos(Oi) 
I 
I 

3 

-J+2 U+7 W-2 cos(0i) 

I 

I 
3 


4. 
4. 
0. 

-1.27492 
6.27492 
3. 

-1.27492 
6.27492 
3. 


Eigenkets and eigenvalues for N e =3 and m s = |. 


No 


Eigenstate 


Energy 


Example 


42 


|3, §, x> r 2) 


3 (J+2 W) 
2 


0. 


Eigenkets and eigenvalues for N e =4 and m s = —1. 


No 


Eigenstate 


Energy 


Example 


43 
44 
45 


|4,-i,2,r 2 ) 

|4,-i,2,r 3 ,i) 

|4,-i,2,r 3 , 2 > 


£ + 2t + U + 5W 
+ -t + U + 5W 
i-t + U + 5W 


6. 
3. 
3. 


Eigenkets and eigenvalues for N e =4 and m s = 0. 


No 


Eigenstate 


Energy 


Example 


46 
47 
48 
49 
50 
51 
52 
53 
54 


|4,o,o,ri) 
|4,o,o,r!) 

|4,0,2,r 2 ) 

|4,o,o,r 3 ,i) 
|4,o,o,r 3 ,i) 

|4,0,2,r 3 ,i) 

|4,o,o,r 3 , 2 ) 
|4,o,o,r 3i2 ) 

|4,0,2,r 3 , 2 ) 


-VA3-J-2 t+3 U+9 W 
2 

s/7Q-J-2t+3U+9W 

2t + U + hW 

-^M- J+t+3 U+9 W 
2 

■/M-J+t+3U+9W 

-t + U + 5W 

-\fAl- J+t+3 U+9 W 
2 

sfAi- J+t+3 U+9 W 

-t + U + 5W 


0.876894 
9.12311 
6. 

4.43845 
8.56155 
3. 

4.43845 
8.56155 
3. 


Eigenkets and eigenvalues for N e =4 and m s = 1. 


No 


Eigenstate 


Energy 


Example 


55 
56 
57 


|4,i,2,r 2 ) 

|4,i,2,r 3 ,i) 

|4,i,2,r 3 , 2 > 


± + 2t + U + hW 
±-t + U + hW 
^-t + U + hW 


6. 
3. 
3. 
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Eigenkets and eigenvalues for N e =5 and m s = — ^ 

. * i 



No 


Eigenstate 


Energy 


Example 


58 
59 
60 


|5,-i|,ri) 
|5,-|,|,r 3 ,i) 


-2t + 2U + 8W 
t+2U+8W 
t+2U+8W 


6. 
9. 
9. 


Eigenkets and eigenvalues for N e =5 and m s = \. 


No 


Eigenstate 


Energy 


Example 


61 
62 
63 


5' Y |' r 3,i) 

|5, 5 , 3,r 3i2 ) 


-2t + 2U + 8W 
t+2U+8W 
t+2U+8W 


6. 
9. 
9. 


Eigenkets and eigenvalues for N e =6 and m s = 0. 


No 


Eigenstate 


Energy 


Example 


64 


|6,o,o,ri) 


3 {U + AW) 


12. 



A. 2 List of abbreviations 

List of abbreviations 



A 2 
A 3 
Ai 
A 5 
A e 
A 7 

A* = - 



4 (2t 2 + tU + U {J-W)) + (J + t-U-W) 2 
J 2 + 27 1 2 + 2 J (U - W) + (U - W f 
J 2 + 36 1 2 + At (U - W) + (U - W) 2 + 2 J (2t + U - W) 
J 2 + 9t 2 -2t (U-W) + (U-W) 2 -2J (t-U + W) 
4 (8t 2 ~2tU + U (J-W)) + (-J + 2t + U + W) 2 
J -2U -7W + 2 cos(0i) 

(cos(6i) - V3 sin(6i)) 

(cos(6i) + V3 sin(6i)) 



J + 2U + 7W+ V / A 2 
J + 2U + 7W+ y/M 



(17) 



0i = 



arccos (^±^)!) 

A 2 2 



A. 3 The eigenvectors 

In the following the unnormaliscd eigenvectors are given in abbreviated form, whereby the indices 
give the number of the state in accordance with the numbering of the energy levels given above. 
For getting the normalised eigenvectors the coefficients have to be divided by the normalisation 
constant Ni, which is given for the vectors containing more than one coefficient. 
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Unnormalized eigenvectors for N e = and m s = 0. 



= 10,0,0,^) 
= i (1000)) 

Unnormalized eigenvectors for N e = 1 and m s = — |. 

I*>2 = I 1 ' |' ri > 

= (|00d) + |0d0) + |d00)) 



l*>3 = |i,-^,r 3 ,i> 

= ±= (|00d) - |0d0» 



l*> 4 - |i,-^,|,r 3 ,2> 

= C 4 _i (|00d) + |0d0» 
+ C 4 _ 2 (|d00» 



C4-1 — — 



1 

71 



G 



4-2 — 



1-2 2 



iV 4 = ^/2C4-i 2 + C 4 

Unnormalized eigenvectors for N e = 1 and m s = 



l*>5 - 



(\00u) + \0uQ) + \u00)) 

V3 



We 



|!> 2' i' 1 " 3 ' 1 ^ 
1 



V2 



(|00u) - |0u0)) 



l*> 7 = |i,|,|,r 3 , 2 > 

= C 7 _i (I 00m) + |0u0» 
+ C 7 _ 2 (|m00)) 
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) 



N 7 = 



\J 2 C7-1 2 + Cj- 2 2 



Unnormalized eigenvectors for N e = 2 and m s = — 1. 



|*) 8 = |2,-i,2,r 2 ) 

= (\0dd) - \dOd) + \ddO}) 
V3 

|*) 9 = |2,-i,2,r 3 ,i) 

= C 9 _i(|0dd» 

+ C9-2 (\dOd) - \ddO)) 



Unnormalized eigenvectors for N e = 2 and m s = 0. 



|*)n = |2,o,o,ri) 

= Cn_i (|002) + |020) + |200)) 

+ Cn_2 (\0du) - \0ud) + \dOu) + \duO) - \uOd) - \udO}) 



|*)i 2 = |2,0,0,ri) 

= C12-1 (|002) + |020) + |200» 

+ C12-2 (|0du) - \0ud) + \dOu) + \duO) - \uOd) - \udO)) 




|*>io = 



|2,-i,2,r 3 , 2 > 

4= (\dOd) + \ddO)) 
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C ia -! = 2yj\t 

C 12 -2 = ( K ^/A-J + 2t-U + W^j 



N 12 = ^/3Ci2-i 2 + 6Ci2-2 2 



|*)i 3 = |2,0,2,r 2 ) 



(\0du) + \0ud) - \dOu) + \duO) - \uOd) + \udO)) 

V6 



|*)i4 = |2,o,o,r 3 ,i) 

= Cu-i (|002) - |020}) 

+ C14-2 (\dOu) - \duO) - \uOd) + \udO}) 

C14-1 = t 

C14-2 = \ (\fM + J + t + u-w) 

N u = v /2Ci4-i 2 +4C 14 - 2 2 

|*) 18 = |2,o,o,r 3l i) 

= Cib-i (|002) - |020» 

+ Ci 5 _2 (|dOu) - |duO) - |uOd) + |udO» 

Cl5_l = t 

Cl 5 _2 = J (-V^+J + t + tf- 



= V /2Ci 5 _i 2 +4C: 



15-2 2 



|*)i 6 - |2,0,2,r 3 ,i) 

= C 16 _! (|0d«) + |0ud)) 

+ C*i 6 _2 (|dOu) - |duO) + |uOd) - |udO)) 



C16-1 — — 
-1 



1 

V3 



Ci6 



-2 — 



2V3 

AT 16 = ^2C 16 _ 1 2 +4C 16 



2 



|*)i 7 - |2,o,o,r 3 , 2 ) 

= C 17 _! (|002) + |020>) 

+ Ci 7 _2 (|0du) - |0ud)) 

+ Ci 7 _ 3 (|200» 

+ C17-4 (\dOu) + \duO) - \uOd) - \udO)) 
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C17-1 = 


1 / 
1 ( 


C17-2 = 




C17-3 = 




Ci7_4 = 





-1) (VAT-J-t-^ + w) 



7Vi 7 = yj 2 Ci 7 _! 2 + 2 Ci7_ 2 2 + C 17 _ 3 2 + 4 C17- 

|*)i 8 = |2,o,o,r 3 , 2 ) 

= Cis-i (|002) + |020)) 

+ C18-2 (|0d«) - |0ud)) 

+ Ci 8 _ 3 (|200)) 

+ Ci8_ 4 (|d0u) + \du0) - \u0d) - \ud0}) 

Cl8_2 = 

1 

r 



2 

4 



Ci 8 _3 = -p V Ai + J + * + £7 - W 



C 



18-4 



f 



Nis = \J2 Cis-i 2 + 2 C18-2 2 + Cis-3 2 + 4Ci8-4 2 

|*) 19 - |2,0,2,r 3 , 2 ) 

= i (|d0u) + |du0) + \u0d) + \ud0)) 

Unnormalized eigenvectors for N e = 2 and m s = 1. 

|v&) 20 = |2,l,2,r 2 ) 

= (\0uu) - |u0u) + \uu0)) 

v3 

|*)ai = |2,l,2,r 3 ,i) 
= C21-1 (|0«u)) 
+ C21-2 (|«0u) - |««0)) 
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|*> 22 = |2,i,2,r 3 , 2 ) 

= -J= (\u0u) + \uu0)) 



Unnormalized eigenvectors for N e = 3 and m s = — |. 



I*> 23 = |3,-|^,r 2 ) 

= 1 (\ddd)) 



Unnormalized eigenvectors for N e = 3 and m s = — \. 



I*>24 = l 3 '-|.|» r i> 

= -j^ (|02d) + |0d2) + |20d) + |2d0) + |d02) + |d20» 
V6 



l*>26 = \^-\\^) 

= —= (|02d) - |0d2) - |20d) + |2d0) + |d02) - |d20» 
V6 



l*> 26 - |3,-i^,r 2 ) 

= — = (\ddu) + | dud) + |udd)) 
v3 



l*> 27 = |3,-^,r 3il ) 

= C 27 _i (|02d) - |0d2)) 

+ C 27 _ 2 (|20d) - |2d0» 

+ C 27 _ 3 (|d02) - |d20» 

+ C 27 _4 (|ddu) + | dud)) 

+ C 27 _ 5 (\udd)) 
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C27-1 
C27-2 
C27-3 
C27-4 



c, 



N 27 

|*>28 



a 



28-5 



27-5 — 



^J-3t + U-W + 2^y%cos(Q 1 )^ 
(j + 3t + U-W + 2y/A 2 ~cos(G 1 )^ 

(V6i 2 ) 

(A 6 2 - 27 t 2 + 6 J £7 - 3 U 2 + 12 J W - 30 U W - 48 W 2 
+ 12 U cos(6i) + 24 V^T^ cos(9i)) 

( A 6 2 - 27 i 2 + 6 J [/ - 3 U 2 + 12 J W - 30 f7 - 48 W 2 
+ 12 yfA~ 2 U 008(9!) + 24 cos(9i)) 



9V6 



— \j 2 C27-1 2 + 2 C27-2 2 + 2 C27-3 2 + 2 C27-4 2 + C; 



27-5 



+ 
+ 

C*28-l = 
C28-2 = 
C28-3 = 
C28-4 = 



' 3 ' _ 2' I' 1 " 3 ' 1 ^ 
C 28 -i (|02d) - |0d2» 
C 28 -2 (|20d) - |2d0» 
C 28 -3 (|d02) - |d20)) 
C 28 -4 (|ddu) + |dud)) 
C 28 _ 5 (|udd)) 

-t 



_ (^\/3J-3\/^f + V3C/-\/3VK-\/3v / ^cos(9i)-3v / ^2"siri(9i)) 
^= (\/3 J + 3V3t + V3U - V3 W - V3 \[M cos(9i) - 3 sin(9i)) 

- (Vet 2 ) 

(A 8 2 - 27 1 2 -6 AsU 



9V6 



iV 28 = 

|*>29 = 



+ 
+ 
+ 



+ 9 U 2 - 12 A 8 W + 36 U W + 36 W 2 ) 



(V-27t 2 -6A 8 C7 
+ 9 C/ 2 - 12 A 8 W + 36 U W + 36 W 2 ) 

\j 2 C28-1 2 + 2 C28-2 2 + 2 C28-3 2 + 2 C28-4 2 + C28- 



I 3 , 4' r 3>i) 
C29-1 (|02d) - |0d2)) 
C29-2 (|20d) - |2d0)) 
C29-3 (|d02) - |d20)) 
C29-4 (|ddw) + I dud)) 
C29-5 (|wdd)) 
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C29-I 
C29-2 
C29-3 
C29-4 



-t 



a 



29-5 



- (^\/3J-3\/3t + \/3^-\/3VK-\/3v / ^cos(ei) + 3v / ^'sin(ei)) 
j-^= (VSJ + 3V3t + V3U -V3W - V3 y/A~ 2 cos(6i) + 3 V^b sin(6i)) 
- (Vei 2 ) 

(A 7 2 - 27t 2 -6A 7 [7 



9^6 



7V 29 = 



+ 9 U 2 - 12 A 7 W + 36 C/ + 36 W 2 ) 



(A 7 2 - 27t 2 -6A 7 U 
+ 9 U 2 - 12 A 7 W + 36 17 W + 36 W 2 ) 
^/ 2 C29-1 2 + 2 C29-2 2 + 2 C29-3 2 + 2 C29-4 2 + C*29-5^ 



I*) 



30 





|3,- 


1 3 r ) 






C30- 


-1 (|02d) + 


|0rf2}) 


+ 


C30- 


-2 (|20d) + 


|2d0)) 


+ 


C30- 


-3 (|d02) + 


|d20)) 


+ 


C30- 


-4 (|c£du) — 


\dud)) 



C30-1 
C30-2 



30-3 — 



C30-4 
iV 30 



6 
1 

18 



(-2J + 9t-2J7 + 2VK + 2 i/M cos(6i) 
(-A 6 2 + 3 J 2 - 12 Jt + 27t 2 -l2tU + 6U 2 -18 JW + 12tW + 24UW 



18 



+ 51 VK 2 + 6 V^b J cos(6i) + 12 yJMt cos(6i) 
-6 V^b t/ cos(0i) - 30 yiaV cos(6i)) 

(A 2 -3 J 2 + 6 Jf + 6t?7-6[/ 2 + 18 JW -6tW -24UW 



-51 W 2 - 6 y/%J cos(6i) - 6 cos(6i) 
+ 6 v^C/ cos(6i) + 30 y/T 2 W cos(0i)) 

Y (-J + 3* — C7 + W -2 v/As" cos(Gi)) 

+ C30-2 + C30-3 + C30-4 



I*) 



31 



I 3 ' 7' r 3,2/ 



1 3 

2'4' 

= C31-1 (|02d) + |0d2» 
+ C 3 i- 2 (|20d) + |2d0» 
+ C31-3 (|d02) + |d20» 
+ C 31 -4(\ddu)-\dud)) 
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C31-1 
Cm_2 



O. 



31-3 



C31-4 
N 31 



t 

6 
1 

36 



(2J-9t + 2[/-2VK + cos(6i) + V3 V^b sin(0i)) 

-2A S 2 + 6J 2 -9A s t-33Jt + 54:t 2 + 6A S U + 6JU -6tU 
+ 30 A 8 W - 6 JW + 87 1 W - 54UW 
-108 W 2 -6 y/A^J cos(0i) - 3 cos(0i) 
-6\/3 -/^J sin(0i) - 3V3 yf~A 2 t sin(0i)) 

= — (2A 8 2 -6J 2 + 9A 8 t + 21Jt-6A 8 U -6JU -6tU 
36 

-30 A s W + 6 J W - 75 1 W + 54 J7 W + 108 VK 2 + 6 V^b J cos(©i) 
-3 \f~A 2 t cos(0i) + 6V3 sin(0i) - 3\/3V^M sin(©i)' 

= ^ (-J + 3t-£/ + W + v /l^cos(ei) + V3 v /I^sin(ei)) 

^ C31-1 2 + C*3i_2 2 + C3I-3 2 + C*3i_4 2 



I*) 



32 



— I 3 i _ 2' i'^ 3,2 ^ 

= C32-1 (|02d) + |0d2)) 

+ C32-2 (|20d) + |2d0» 

+ C32-3 (|d02) + |d20» 

+ C32-4 (|ddu) - I dud)) 



C32-1 
C32-2 



32-3 



C32-4 
^32 



^ (2J-9t + 2C/-2VK + V^b cos(©i) - V3 sin(0i)) 

— (— 2 A 7 2 + 6 J 2 — 9 A 7 1 — 33 J t + 54t 2 + 6 A 7 U + 6 J U 
36 

-6 i U + 30 A 7 TV - 6 J TV + 87 1 W - 54 U W - 108 TV 2 - 6 yf~A 2 J cos(0i) 
—3 \fA~2~t cos (0i) + 6 a/^ J sin(0i) + 3 V3 ^A~ 2 t sin(©i) 

(2 A 7 2 -6 J 2 + 9 A 7 1 + 21Jt-6 A 7 U — 6JU — 6tU — 30 A 7 W 



36 



+6 J TV - 75 1 TV + 54 U W108 W 2 + 6 ^f~A 2 J cos(©i) 
-3-/<M cos(0i) - 6>/3\/^J sin(9i) +3V3 sin(0i)) 

Y (-J + 3t-C/ + TV + v / ^ 2 "cos(ei)-\/3v / A2"sin(ei)) 



\Fl \ C32-1 2 + C32-2 2 + C32-3 2 + C32- 



Unnormalized eigenvectors for N e = 3 and m s = ^. 



33 



(|02u) + |0u2) + |20u) + |2uO) + |u02) + |u20)) 

v6 



24 



Schumann: Extended Hubbard model on small clusters 



l*>34 - |3,|,|,r 2 > 



2' 4' 

(|02u) - |0u2) - |20u) + |2u0) + |u02) - |u20)) 

V6 



l*> 35 - |3,^,r 2 ) 



—= (\duu) + \udu) + \uud}) 
v3 



I*) 



30 











3 ' 2' 4' r 3' 1 ' > 






C 36 ^ (|02u) - 


|0u2» 


+ 


C 3 6-2(|20 M )- 


\2u0)) 


+ 


C 36 -3 




+ 


C 36 -4(|w02)- 


\u20)) 


+ 


C 36 _ 5 (|udu) + 


\uud)) 



C 36 -i = — (-2 J + 9t-2J7 + 2VK + 2VXcos(ei) 

C 36 - 2 = 4 (-2A 6 2 + 6A J-9A 6 t + 33 Jt 
36 v 

+ Mt 2 -6A 6 U +18 JU + 6tU -30A 6 W + 36 JW -87tW -54UW 
-108W 2 -6^/%t cos(0i)) 

C 36 _ 3 = ^ (j + 9t + U-W + 2^008(0^ 

C 36 _4 = ^- (-2A 6 2 + 6A 6 J -9A 6 t + 21Jt 
36 

+ 108t 2 -6 A 6 U +18 JU -6tU -30 A 6 W + 36 JW -75tW -54UW 
-IQ8W 2 + 6^fA^t cos(Qi) 



C 3 6-5 - ^ (j + 9i + C7- 1^ + 2^003(61) 
^36 = \/ 2 C36-1 2 + 2 C 3 6-2 2 + C 3 6- 3 2 + 2 C36-4 2 + 2 C36-5 2 



l*>37 = |3,|,|,r 3 ,i> 

= C37-1 (|02u) - |0u2)) 

+ C 37 _ 2 (|20w)- |2«0)) 

+ C 3 7-3(|d««» 

+ C 37 ^(\u02) - \u20)) 

+ C 37 _ 5 (\udu) + \uud}) 
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C37-I - 

C37-2 = 

C37-3 = 
C37-4 = 

C37-5 = 
N37 = 

|*>38 = 

+ 
+ 
+ 
+ 

C38-I = 

C38-2 = 

C*38-3 = 
C38-4 = 

C38-5 = 
^38 = 



- (2.J -9t + 2U -2W + yfM cos(6i) + V3 ^JT 2 sin(6i)) 

— (-2 A s 2 -6 A 8 J + 9 A 8 t + 33 Jt 
36 

+ 54 1 2 + 6^ 8 U+ 18 JU + 6tU + 30 A 8 W + 36 J W - %7tW - 54 UW 
-108 W 2 + 3 ^/%t cos(6i) + 3 V3y/A^t sin(6i)) 

f - (- J - 9* - U + W + \/M cos(6i) + V3 yj% sin(6i)) 

4 (-2 As 2 -6A s J + 9A s t + 2lJt 
36 v 

+ 108 t 2 + 6 A S U + 18 J U — 6tU + 30 A$W + 36 J W — 75tW — 54UW 
-108 VF 2 - 3 V^M cos(6i) - 3 V3^/%t sin(9i)) 

Y ( y -J-9t-U + W + y / A^cos(Q 1 ) + V3y / %sm(e 1 )^ 

Y 2 C37_i 2 + 2 C37_2 2 + C < 37_3 2 + 2 C37_4 2 + 2 C37_5 2 



I 3 ; i' 1 " 3 ' 1 ^ 
Css-i (|02u) - |0u2» 
C 38 - 2 (|20 U )- |2«0)) 
C 38 - 3 (\duu)) 
C 38 -4(|w02) - |«20)) 
C 38 -5 (\udu) + \uud)) 



— (2J-9t + 2U-2W+ y/A~ 2 cos(6i) - V3 ^/M sin(6i)) 

— (-2A 7 2 -6A 7 J + 9A 7 t + 33Jt 
36 v 

+ 54 1 2 + 6 A 7 U + 18 J U + 6 i U + 30 A 7 W + 36 J - 87 1 W - 54 U W 
-108 VK 2 + 3 V^M cos(6i) - 3 V3^/A^t sin(6i)) 

| (- J - 9 1 - J7 + W + ^/M cos(6i) - V3 V^b sin(6i)) 

— (-2A 7 2 -6^ 7 ^ + 9A 7 t + 21Jt 
36 v 

+ 108 1 2 + 6 A 7 U + 18 J U - 6 i [/ + 30 A 7 W + 36 J W - 75 1 W - 54 f7 W 
-108 W 2 - 3 y/A~ 2 t cos(0i) + 3 \/3y/A^t sin(9i)) 

Y (-J-9t-C/ + T4^ + v / ^cos(ei)-\/3v / ^sin(ei)) 
2 C38-1 2 + 2 C38-2 2 + C38-3 2 + 2 C38-4 2 + 2 C38-5 2 



26 



Schumann: Extended Hubbard model on small clusters 



|*>39 = 



|3, -, ^,r 3 ,2) 

C39-1 (|02u) - 
C 39 - 2 (|20 U )- 
C 39 -3(|w02)- 
C39-4 (\udu) - 



|0u2)) 
|2u0» 
|u20)) 
\uud)) 



C39-1 

C3Q-2 



Q 



39-3 



G39-4 
^39 



= -y (-2 J + 9t-2J7 + 2^ + 2v/l2"cos(ei)) 



6 

4 (-A 6 2 + 3 J 2 - 12 Jt + 27t 2 -l2tU + 6U 2 -18 JW + 12tW + 24UW 
18 

+ 51 VK 2 + 6 V^b J cos(0i) + 12 y]~K~ 2 t cos(6i) 
-6 V^b t/ cos(6i) - 30 V^b^ cos(9i)) 

— U 6 2 -3 J 2 + 6 Jt + 6t?7-6t/ 2 + 18 JV^-6tVF-24[/iy 
18 

-51 - 6 i/I^ J cos(6i) - 6 y/A~ 2 t cos(6i) 
+ 6 V^bt/ cos(0i) + 30 j~A 2 W cos(0i)) 

Y (-J + 3t-C/ + T^-2 v /A 2 "cos(ei)) 



= V2 

\/ C39-1 2 + C*39_2 2 + G39-3 2 + C39_ 



l*>4o = |3,|,|,r 3)2 > 

= C40-1 (|02u) + |0u2>) 

+ C 4 o- 2 (|20u) + |2 U 0)) 

+ C i0 -3(\u02) + \u20)) 

+ C40-4 (|udu) — \uud)) 



C40-1 
Go -2 



G 



40-3 



G40-4 
^40 



^ (2J-9i + 2[/-2W + V^b cos(6i) + V3 V^b sin(6i)) 

— (-2A S 2 + 6J 2 -9A 8 t-33Jt + 54:t 2 + 6A 8 U + 6JU 
36 

-6 1 U + 30 A 8 W - 6 J W + 87 1 W - 54 U W - 108 W 2 - 6 J cos(Gi) 
-3 \/~A 2 t cos(6i) - 6V3 V^b >/ sin(0i) - 3\/3 \[A~ 2 t sin(6i)) 

(2A 8 2 -6 J 2 + 9 A 8 t + 21Jt- 6 A 8 U -6 JU -6tU -30 A 8 W 



36 



+6 JW - 75tW + 54UW + 108 W 2 + 6 ^/T 2 J cos(6i) 

—3 \fA~2~t cos (6i) + 6-^/3-/^ J sin(6i) - 3\/3 sin(6i)) 

Y (-J + 3t-C/ + T^ + v / Xcos(ei) + \/3v / A2"sin(ei)) 
V2 \l G40-1 2 + C40-2 2 + Go-3 2 + C40-4 2 
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|*>41 = 



|3, -, ^,r 3 ,2) 



= C41-1 (|02u) + |0u2» 

+ C 4 i-2(|20u) + |2u0)) 

+ G 4 i_ 3 (|u02) + |u20)) 

+ 641-4 (\udu) — \uud)) 



C41-1 
C41-2 



f 

6 
1 

36 



G 



41-3 — 



36 



C41-4 
N 41 



(2J-9t + 2U-2W+ -/A* cos(6i) - V3 ^~T 2 sin(0i)) 

(-2 A? 2 + 6 J 2 - 9 A 7 1 - 33 J i + 54 1 2 + 6 A 7 U + 6 J U - 6 t [7 + 30 A 7 W 
-6 JW + 87 1 W - 54UW - 108 TU 2 - 6 yJ~A~ 2 J cos(0i) 
-3 cos(Gi) + 6 V3 J sin(0i) + 3 V3 V^b * sin(6i; 

(2A 7 2 -6J 2 + 9^ 7 i + 21Jt-6^ 7 ?7-6J[/-6tC/-30A 7 T^ + 6JV^ 
-75 i W + 54 £7 TU + 108 W 2 + 6 y]~A 2 J cos(6i) 
-3 cos(Gi) - 6\/3 -/^ J sin(9i) +3\/3 sin(©i)) 

(- J + 3t-U + W + ^/A~ 2 cos(6i) - V3 sin(6i)) 



2 ' 



42 



I*) 



43 



I*) 



44 



C44-I 
C44-2 

N44 

|*>48 



— V2 y G41-1 2 + C41-2 2 + C41-3 2 + G41-4 2 

Unnormalized eigenvectors for N e = 3 and m s 

is 3 15 r > 

= 1 

Unnormalized eigenvectors for N e = 4 and m s = — 1. 

= |4,-i,2,r 2 ) 

= {\2dd) - \d2d) + \dd2)) 

v3 

= |4,-l,2,r 3 ,i) 
= C44-1 (|2dd)) 
+ C44-2 (|d2d) - |dd2)) 

^2 

3 
1 

V6 

= \J C44-1 2 + 2 C44-2 2 

= |4,-i,2,r 3 ,2> 

= (|d2d) + |dd2» 
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Schumann: Extended Hubbard model on small clusters 



Unnormalized eigenvectors for N e = 4 and m s = 0. 

|*>46 = 14,0,0,^) 

= C46-1 (| 022) + 1 202) + 1 220)) 

+ C 46 - 2 (|2du) - \2ud) + \d2u) + \du2) - \u2d) - \ud2)) 



C46-1 — 2 y - t 

C i6 - 2 = ^= + J + 2t + U- W^j 



N i6 = ^3C 4 6-l+6C 4 6-2' 

|*>4 7 = |4,o,o,ri) 

= C47-1 (|022) + |202) + |220>) 

+ C47-2 (\2du) - \2ud) + \d2u) + \du2) - \u2d) - \ud2}) 

C47-1 = 2j|i 



a 



47-2 — 



2V6 



)= (VM- J -2t-U + W^j 



N i7 = ^dr-! 2 + 6C 47 . 



2 

2 



|*>48 = |4,0,2,r 2 ) 

= —= (\2du) + \2ud) - \d2u) + \du2) - \u2d) + \ud2)) 
V6 

|*) 49 = |4,o,o,r 3 ,i> 

= C49-1 (| 202) - |220» 

+ C49-2 (\d2u) - \du2) - \u2d) + \ud2}) 



C49-1 = 7^= (yA 4 - J + t-u + w) 
C49-2 = — ' — 



\V2 



N49 = \/2C49_i 2 +4C 49 _2 2 



|*>5o = |4,0,0,r 3 ,i) 

= C50-1 (|202) - |220» 

+ C50-2 (\d2u) - \du2) - \u2d) + \ud2)) 

Cbo-i = (y/Al + J-t + U- W) 



t 
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I*>B1 = 
+ 

C51-1 = 
C5I-2 = 

N 51 = 

I*) 52 = 

+ 
+ 
+ 

C*52^1 = 
C52-2 = 
C52-3 = 
C*52-4 = 

N 52 = 

|*>53 = 

+ 
+ 
+ 

C53-I = 
C53-2 = 
C53-3 = 
C53-4 = 

N 53 = 

|*>54 = 



14,0,2,1V) 

C51-1 (|2du) + |2ud» 

C51-2 (|d2u) - |du2) + \u2d) 

1 

1 



id2)) 



2V3 

y / 2C 5 i-i 2 +4C 5 i-2 2 

|4,o,o,r 3 , 2 ) 

C52-1 (| 022)) 
C52-2 (| 202) + |220» 
C52-3 (|2d«) - |2ud)) 
C52-4 (|d2u) + |du2) - |u2d) 



id2)) 



2t 

71 



f 

7s 



A 4 + J -t + U 



- ^(^I + j-t + cz-tt-) 

^/ C52-1 2 + 2 (C52-2 2 + C52-3 2 + 2 C52-4 2 ) 

|4,o,o,r 3 ,2> 

C53-1 (1 022)) 

C53-2 (1 202) + |220» 

C 53 - 3 (|2d«) - \2ud)) 

C53-4 (|d2u) + |du2) - |u2d) - |ud2)) 



2t_ 

71 



71 



^ (yfc-J + t-U + W) 



2V3 
-1 



4 V3(^- J + t - t/ + ^) 
^/ C53-1 2 + 2 (C 53 _ 2 2 + C53-3 2 + 2 C 53 _4 2 ) 



|4,0,2,r 3 , 2 ) 

i (|d2u) + |du2) + |u2d) + \ud2)) 
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Schumann: Extended Hubbard model on small clusters 



Unnormalized eigenvectors for N e = 4 and m s = 1. 



I*) 



55 



= |4,i,2,r 2 
i 



V3 



(\2uu) - \u2u) + |uu2» 



|*) B6 = |4,l,2,r 3 ,i) 
= C 56 _i(|2uw)) 
+ C 56 _ 2 (\u2u) - \uu2)) 



C56-1 
C56-2 

^56 



'2 
3 

1 

7i 



C. 



56-1 



2 c. 



56-2 



|*) 57 = |4,i,2,r 3 , 2 > 



(|u2u> + |im2» 



Unnormalized eigenvectors for N e = 5 and m s = — |. 



-t= (|22d) + |2d2) + \d22)) 
v3 



l*>59 - |5,-^,r 3 ,i) 

= ^ (|22d) - |2d2» 

l*>eo - |5,-^,r 3 , 2 > 

= Ceo-i (|22d) + |2d2» 
+ Ceo-2 (|d22)) 

1 

7! _ 

\J 2 Cqq-1 2 + C$0-2 2 



C'eo-i = 

C*60-2 = 

N 60 = 
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Unnormalized eigenvectors for N e = 5 and m s = -| 



2 ' 



l*>ei = |5,W' r i) 



2'4' 

(|22u) + \2u2) + \u22}) 

v3 



l*> 62 = I5,|,|,r 3 ,i> 

= (|22«) - |2u2» 



l*>63 = 15,-^^3,2) 

= Cea-i (|22u) + |2«2)) 
+ C 6 3- 2 (|w22)) 



63-1 



1 

^2 



63-2 — 



^63 — \/2 C 6 3-l 2 + C63- 2 Z 



Unnormalized eigenvectors for N e = 6 and m s = 0. 



|*) 64 - |6,o,o,ri) 

= 1(|222» 



B Eigensystem of the tetrahedron 

In the following tables the first column gives the number of the state for referencing in the text. The 
second column labels the eigenstate by the eigenvalues of the U-independent symmetry operators, 
i.e. the electron occupation number N e , the spin-projection S z in z-direction m s , the eigenvalues 
s(s + 1) of S 2 and the spatial symmetry is indicated by T^j, where the first index labels the 
irreducible representation of the tetrahedral group and the second numbers the partner. The 
notation is based on Rcf. [19]. The third column gives the enery eigenvalues in abbreviated 
form. The abbreviations are listed subsequent to the tables. In the last column a numerical value 
is given for example. For comparison with we have chosen the same parameters as in Ref. [8] 
corresponding to a pure Hubbard-model with U = 5t, W = J = Ot. If the grand-canonical energy 
levels in an applied magnetic field are needed one has to substract fj,N e + h*m s . 



B.l The eigenvalues 

Eigenkets and eigenvalues for N e =0 and m s = 0. 



No 


Eigenstate 


Energy 


Example 


1 


10,0,0,^) 





0. 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =l and m s = — |. 



No 


Eigenstate 


Energy 


Example 


2 
3 
4 
5 


1,-2' |, r i) 
|1,-|, |, r 4 : l> 

1, -5, 4, ^4,2) 
I 1 ' _ |) |, r 4,3) 


3t 
-i 


3. 
-1. 
-1. 
-1. 


Eigenkets and eigenvalues for N e =l and m s = |. 


No 


Eigenstate 


Energy 


Example 


6 
7 

8 
9 


1, 2' |, r i) 
1, \i 4^4,1} 
|1, 5, |, r 4,2) 
|1, 5> |, r 4,3) 


3t 
-t 
-t 
-i 


3. 
-1. 
-1. 
-1. 



Eigenkets and eigenvalues for N e =2 and m s = — 1. 



No 


Eigenstate 


Energy 


Example 


10 


|2,-i,2,r 4 ,i) 


± + 2t + 2W 


2. 


11 


|2,-i,2,r 4>2 ) 


+ + 2t + 2W 


2. 


12 


|2,-i,2,r 4; 3> 


j + 2t + 2W 


2. 


13 


|2,-i,2,r 5 ,i) 


j -2t + 2W 


-2. 


14 


|2,-i,2,r 5;2 ) 


j -2t + 2W 


-2. 


15 


|2,-i,2,r 5;3 ) 


j -2t + 2W 


-2. 
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Eigenkets and eigenvalues for N e =2 and m s = 0. 



No 


Eigenstate 


Energy 


Example 




16 




2, 0,0, 


^a-|+2t + f + W 






1. 




17 




2,0,0,1^) 


j + 2t+f + W 






8. 




18 




2,o,o,r 3 ,i) 


-J-2t + 2W 






-2. 




19 




2,o,o,r 3 , 2 ) 


-J-2t + 2W 






-2. 




20 




2,o,o,r 4 ,i) 


^i-l + IL + W 






-0.701562 




21 




2,o,o,r 4 ,i) 


^~j + ^ + W 






5.70156 




22 




2, 0, 2, r 4i i) 


2t + 2W 






2. 








2,o,o,r 4 , 2 ) 








-0 701 562 




24 




2,o,o,r 4 , 2 ) 


^L-l + ll + W 






5.70156 




25 




2,0,2,r 4 , 2 ) 


2t + 2W 






2. 




26 




2,o,o,r 4i3 ) 


+ + W 






-0.701562 




27 




2,o,o,r 4 , 3 ) 








5.70156 




28 




2,0,2,r 4 , 3 ) 


2t + 2W 






2. 




29 




2,0,2,r 5 ,i) 


-2t + 2W 






-2. 




30 




2,o,2,r 5 , 2 ) 


-2t + 2W 






-2. 
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2,o,2,r 5 , 3 ) 


-2t + 2W 






-2. 






Eigenkets and eigenvalues for N e =2 and m s = 


1. 




No 


Eigenstate 


Energy 


Example 




32 


|2,i,2,r 4 ,i) 


± + 2t + 2W 






2. 




33 


|2,i,2,r 4 , 2 ) 


+ + 2t + 2W 






2. 




34 


|2,i,2,r 4 , 3 > 


± + 2t + 2W 






2. 




35 


|2,i,2,r B ,i) 


4 - 2t + 2W 






-2. 




36 


|2,i,2,r 5 , 2 ) 


4 - 2t + 2W 






-2. 




37 


|2,i,2,r B , 3 > 


4 - 2t + 2W 






-2. 






Eigenkets and eigenvalues for N e =3 and m s = 


3 
2" 




No 


Eigenstate 


Energy 


Example 




38 


|3,-l,f,r 2 ) 


^ -3i + 6W 






-3. 




39 


3, -f , X' r 5,i) 


^ +t+6W 






1. 




40 


3, — | , 


kl+t + 6W 






1. 




41 


3,-§, X' r 5,3) 


^-+t + 6W 






1. 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =3 and m s = — |. 



No 




Eig 


enstate 


Energy 


Example 


42 


|3, 


i 

~ 2 


|.Ti> 




t + U + 4W 


6. 


43 


|3, 


~2 






^ — 61 + W 


-3. 


44 


|3, 


1 

— 2 


45 1 3,1/ 




_ J _ t + £ + 5iy 


-0.791288 


45 


|3, 


~ 2 


f,r 3 ,i> 




^-■f-i+f + 5V^ 


3.79129 


46 


|3, 


1 
2 


45 r 3,2) 






-0.791288 


47 


|3, 


1 
2 


3 r \ 
4: 1 3,2/ 




^ - 4 - * + T + 5 w 


3.79129 


48 


|3, 


1 
2 


3 F \ 
41 1 4,1/ 


-,/ 
3 


1 j 2 1/ | 14 2V^5Cos(e 4 ) 


-0.248211 


49 


|3, 


1 
2 


3 r \ 
4: 1 4,1/ 




3 


8.51895 


50 


|3, 


1 
2 


3 r \ 

4: 1 4,1/ 




A20 
3 


4.72926 


51 


|3, 


1 
2 


3 r \ 
41 1 4,2/ 


—,/ 
3 


| 2 u 14 Vv ^ V ^-5 COS(t^4j 


-0.248211 


52 


|3, 


1 
2 


41 F 4,2) 




^121 

3 


8.51895 


53 


|3, 


— 2 


3 T 1 \ 

4' 1 4,2/ 




3 


4.72926 


54 


|3, 


1 
2 


41 F 4,3) 


—J 
3 


| £ 2U 1 14 W ^ V cos^4j 


-0.248211 


55 


|3, 


1 
2 


3 r \ 
4: 1 4,3/ 




A21 
3 


8.51895 


56 


|3, 


_ 2 


3 T 1 \ 

4; 1 4,3/ 




A20 
3 


4.72926 


57 


|3, 


1 

2 


45 r 5,l) 






-2.40512 


58 


|3, 


1 

— 2 


§5 F 5,l) 




- 1 - 1 + T + 5 W 


5.40512 


59 


|3, 


1 
2 


f,r 5ll } 






1. 


60 


|3, 


1 
2 


3 r \ 
45 1 5,2/ 






-2.40512 


61 


|3, 


1 
2 


3 r \ 
45 1 5,2/ 






5.40512 


62 


|3, 


1 
2 


f,r 5 , 2 ) 




f +t + 6T^ 


1. 


63 


|3, 


1 
2 


3 r \ 
45 1 5,3/ 




+ 5V^ 


-2.40512 


64 


|3, 


1 
2 


3 r \ 
45 1 5,3/ 






5.40512 


65 


|3, 


1 
2 


15 r \ 

— 5I 5,3/ 






1. 



adp header will be provided by the publisher 



35 



Eigenkets and eigenvalues for N e =3 and m s = |. 



No 


Eigenstate 


Energy 


Example 


66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 


I 3 ) 5' l' r l) 

M,f,r 2 ) 
|3, 5, |,r 3) i) 
|3, |, f ,r 3) i) 
|3, 2, i,r 3 , 2 ) 
|3, 5 , f.r 3 ,2) 
|3, 25 1) r i 4 i i) 

|3, 5) i ) r.4,1) 

|3, 5, f ,r 4 ,i) 
|3, 5, |,r 4)2 ) 
|3, |, §,r 4)2 ) 

3, 2' l' r 4,2> 

3, 5, 1^4,3) 

3, 5. 4) r 4, 3 ) 
3, 2> f j ^4,3) 

3, 5, f. r 5,i) 

|3,5,f,r 5 ,i} 
|3,|,^,r 6 ,i> 

|3, 5, 4^5,2) 

3, j) 4) r 5i2 } 
|3,if,r 5i2 ) 

|3, 5, |,r 5)3 ) 

3, 5, I' r 5,3) 
I 3 , |,X' r 5.3> 


t + U + AW 
+ -3t + 6W 
I t + u 2 + 5 W 

^¥ - i - 1 + f + 5 w 

^- + -t+^ + 5W 

-,/ | . | 2U I 14 W 2V^Icos(e 4 ) 
3 ~r 1 ' 3 ' 3 3 

3 
A20 
3 

— J | . | 2U | 14 W 2 V / ^I cos(0 4 ) 
3 ' 1 ' 3 ~r 3 3 
A21 

3 
A20 
3 

-,/ | . | 2 [7 | 14 W 2V^Jcos(e 4 ) 
3 ~T L ' 3 ' 3 3 
A21 

3 
A20 
3 

2 - i - 1 + T + 5 ^ 
4 +t + 6W 

j+t+6W 

- 4 - 1 + T + 5 W 
j+t+6W 


6. 

-3. 

-0.791288 
3.79129 

-0.791288 
3.79129 

-0.248211 
8.51895 
4.72926 

-0.248211 
8.51895 
4.72926 

-0.248211 
8.51895 
4.72926 

-2.40512 
5.40512 
1. 

-2.40512 
5.40512 
1. 

-2.40512 
5.40512 
1. 


Eigenkets and eigenvalues for N e =3 and m s = |. 


No 


Eigenstate 


Energy 


Example 


90 
91 
92 
93 


|3, f, — ,T 2 ) 

3, 5> — ^5,1) 
3, j) X > r 5,2) 

I 3 , §,— ,r 5 ,3> 


^ -3i + 6T4^ 
^ +t + 6W 
^-+t + 6W 
^-+t + 6W 


-3. 
1. 
1. 
1. 


Eigenkets and eigenvalues for N e =4 and m s = —2. 


No 


Eigenstate 


Energy 


Example 


94 


|4,-2,6,r 2 ) 


3 J+ 12W 


0. 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =4 and m s = — 1. 



No 


Eigenstate 


Energy 


Example 


95 




4—19 Fo\ 

L t z j 1 2/ 


| + U + 10W 


5. 


96 




4—16 r-,) 


2 i 


0. 


97 




4-i 2 r^i 

^1 i l ^7 1 3,1/ 


| + C/ + low 


5. 


98 




4 -i 2 i\ 9) 


4 + [/ + low 


5. 


99 




4 -1 2 Tn) 

i ; J-j 1 4,1/ 


{ + U + 10W 


5. 


100 




4,-i,2,r 4 , 2 > 


i + u + iow 

2 


5. 


101 




4,-i,2,r 4>3 ) 


i + u + iow 

2 


5. 


102 




4,-i,2,r B ,i) 


An 
6 


-1.5136 


103 




4,-i,2,r B ,i) 


A ls 
3 


8.34789 


104 




4,-i,2,r B ,i) 


dm 

3 


3.16571 


105 




4,-i,2,r 5;2 ) 


An 
6 


-1.5136 


106 




4,-i,2,r B ,2> 


^18 

3 


8.34789 


107 




4,-i,2,r B , 2 > 


^16 

3 


3.16571 


108 




4,-i,2,r B , 3 > 


An 
6 


-1.5136 


109 




4,-i,2,r B;3 ) 


dia 

3 


8.34789 


110 




4,-i,2,r B;3 ) 


Al 6 

3 


3.16571 
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Eigenkets and eigenvalues for N e =4 and m s = 0. 



No 


Eigenstate 


Energy 


Example 


111 


1 A 

|4 







- 4 + +9W 


2.78301 


112 


|4 





0, Ti) 


2 - 1 + ¥■ + 9W 


12.217 


1 1 Q 


4 





2,r 2 ) 


U + W W 


r 
0. 


114 


1 A 

|4 





6, r 2 ) 


J + 12 W 


0. 


115 


|4 





o,r 3 ,i) 


7- . 7-7- inn? 2 a/A7 cos(Oi ) 

- J + (7 + 10 TV y — ^ ' 


-1.84429 


116 


|4 





o,r 3 ,i) 


J+U + 10W+ ^™ s ( e i) + VATsin(ei) 


10.8443 


117 


|4 





o,r 3 ,i> 


J+[/ + 10VK+ ^™ s(6l) V^sin(6i) 


6. 


118 


|4 





2,r 3 ,i) 


U + 10W 


5. 


119 


|4 





o,r 3 , 2 ) 


j + u + iow 2 v^r^(0i) 


-1.84429 


120 


|4 





o, r\ o) 

? o,z/ 


j + [/ + iovk+ v^Tcosoo + ^7 sin ( e ) 

V3 


10.8443 


121 


Izl 


n 
u 


U, 1 3,2/ 


J + [/ + 10VK+ y7rrc ° s(6l) A/ATsinfei) 


6. 


122 


IzL 


n 
u 


z i 1 3,2/ 


rj + 10 w 


5. 


123 


Izl 


n 
u 


u i 1 4,1/ 


-2 Am 
3 


1.65211 


124 


Izl 


n 
u 


n f , , \ 

U, 1 4,1/ 


A14 

3 


11.5136 


125 


1/1 


U 


U, 1 4,1/ 


A13 

3 


6.83429 


126 


|4 





2,r 4 ,i) 


rj + 10 w 


5. 


127 


|4 





o, r 4 . 2 ) 


-2A W 
3 


1.65211 


128 


4 


o 


o,r 4 ■>) 


3 


11.5136 


129 


4 


o 


o r 4 9) 


Ms 
3 


6.83429 


130 


|4 





2,r 4 , 2 ) 


c/ + 10 w 


5. 


131 


|4 





0, r 4 , 3 > 


-2Aio 
3 


1.65211 


132 


|4 





0, r 4 , 3 > 


A14 

3 


11.5136 


133 


|4 


o, 


0, r 4 , 3 ) 


A13 
3 


6.83429 


134 


|4 





2,r 4 , 3 > 


C/ + 10 w 


5. 


135 


|4 





o,r 5 ,0 


-J+ 17+ 10 


5. 


136 


|4 





2,r 5 ,i) 


— J | 2C/ | 32 W 2v / A^cos(0 3 ) 
3 ~T 3 ' 3 3 


-1.5136 


137 


|4 





2,r 5 ,i) 


A i7 
3 


8.34789 


138 


|4 





2,r B ,i) 


A15 
3 


3.16571 


139 


|4 





o,r 5 , 2 ) 


-J+ C/ + 10W 


5. 


140 


|4 





2,r 5 , 2 ) 


-J 1 2C/ 1 32 W 2^Mcos(e 3 ) 
3 + 3 + 3 3 


-1.5136 


141 


|4 





2,r 5 , 2 > 


An 
3 


8.34789 


142 


|4 





2,r 5 , 2 ) 


A15 
3 


3.16571 


143 


|4 





o,r 5 ,3> 


-j + [/ + ioiy 


5. 


144 


|4 





2,r B , 3 > 


— J , 2t/ I 32 W 2 v / 3Icos(0 3 ) 
3 " r 3 " r 3 3 


-1.5136 


145 


|4 





2,r B , 3 > 


A17 
3 


8.34789 


146 


|4 





2,r 5 , 3 ) 


AlB 

3 


3.16571 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =4 and m s = 1. 



No 


Eigenstate 


Energy 


Example 


147 




4 1 9 Fo\ 

-Lj z l 1 2/ 


^ 4- TJ 4- 10 VF 




5. 


148 




^1 X > U l 1 2/ 


^ 4- 12 W 

2 ^ 1Z " 




0. 


149 




4i 2 r 5 ii 

^! x j x 3,1/ 


- 4 f7 + 10 




5. 


150 




x ? ^) - 1 - 0,2/ 


- 4- 77 4- 10T4 7 




5. 


151 




4 1 2 L,) 

i ; J- 1 5 -i 4,1/ 


^ + fJ 4- 10 W 




5. 


152 




4, i,2,r 4 , 2 > 


i + u + iow 

z 




5. 


153 




4,i,2,r 4 , 3 > 


+ + U + 10W 




5. 


154 




4,i,2,r B ,i) 


An 
6 




-1.5136 


155 




4,i,2,r 5 ,i) 


Ais 
3 




8.34789 


100 




4,i,2,r 5 ,i) 


Al6 

3 




3.16571 


157 




4,i,2,r 5 , 2 ) 


An 
6 




-1.5136 


158 




4,i,2,r 5 , 2 ) 


^18 

3 




8.34789 


159 




4,i,2,r 5 , 2 ) 


^16 

3 




3.16571 


160 




4,i,2,r B , 3 > 


An 
6 




-1.5136 


161 




4,i,2,r 5l3 ) 


^18 

3 




8.34789 


162 




4,i,2,r B , 3 > 


A le 
3 




3.16571 




Eigenkets and eigenvalues for N e =4 and m s = 


2. 


No 


Eigenstate 


Energy 


Example 


163 


|4,2,6,r 2 ) 


3 J+ 12W 


0. 




Eigenkets and eigenvalues for N e =5 and m s = 


3 


No 


Eigenstate 


Energy 


Example 


164 


|5,-l,f,r 2 ) 


^ + 3t + [/+ 18T4 7 




8. 


165 


I 5 ' ~f ' T' r M) 


Y -t + tf + is w 




4. 


166 


5 ; — f i X' 1 " 5 ^) 


-t + U + 18W 




4. 


167 


|5,-|, X' r 5,3> 


-t + U + 18W 




4. 
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Eigenkets and eigenvalues for N e =5 and m s = — |. 



No 



Eigenstate 



Energy 



-t + 2U + 16W 
4 + 3t + U+ 18W 

^-{+t+^- + 17W 
=^-j+t+^ + 17W 
17W 



Example 



4 
3 
4 
3 
4 
3 
4 
3 
4 

3 
4 

3 
4 
3 
4 

3 
4 

3 
4 

3 
4 

3 
4 

3 
4 

3 
4 
3 
4 
3 
4 
15 

4 
3 
4 
3 
4 
15 

4 
3 
4 
3 
4 
15 



168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 



IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
IV 
15,- 



,r 2 ) 
r 3 ,i 
r 3 ,i 
r3,2 
r3,2 
r 4 ,i 
r 4,i 
r 4 ,i 

^4,2 
^4,2 
^4,2 

r 4 , 3 
r 4 , 3 
r 4 , 3 
r 5 ,i 
r 5 ,i 
,r 5 , 
r 5 , 2 
r 5 , 2 
.r 6 , : 



r 5 , 3 ; 
r 5 , 3 / 

X' r 5,3) 



3 U 



^2^ 2 t 2 
—L _ i _|_ 5£ i 50 W _ 2 cos(e 5 ) 
3 1 " r 3 " r 3 3 
-,/ _ j. I 5U_ , 50 W , V^e cos(6 5 ) , sin(6 5 ) 
3 L ^ 3 ^ 3 ^ 3 ~T v/3 

— J _ j. I 5J7 , 50 TV , cos(e 5 ) _ sin(65) 

3 ' 3 3 Q ^ 



3 1 ~r 3 



3 V3 
50 W _ 2 yC57 cos(B 5 ) 



3 " 1 3 1 3 3 

— J + i 5U i 50 TV , V'As" cos(e 5 ) , sin(6 5 ) 

3 1 " r 3 " r 3 " r 3 " r 

-.7 _ j. I 5U_ , 50 TV , v^6 cos(e 5 ) _ sin(6 5 ) 
3 1 + 3 ^ 3 ^ 3 ^3 

— J j. , 5U , 50 _ 2 yG47 cos(e 5 ) 

3 " r 3 " r 3 3 

=Z — 7 4- £H 4- 50 TV , y^e cos(e 5 ) , y^e sin(6 5 ) 

3 t+ 3^3 + 3 ~T ^3 



50 W i \Me cos(6 5 ) _ V^6 sin(e 5 ) 
3 + 3 y/3 



*4 



7 



-t+^f + 17W 



2 2 

±+t+^ + 17W 

± - t + U +18W 

=^ - j +t+2£ + 17W 

j+t+^ + 17W 

{ -t + U+ 18 W 

=^±-{+t+^ + 17W 

^ - ±+t+^ + 17W 
j 

2 



2 2 

±-t + U + 18W 



9. 
8. 

4.59488 
12.4051 
4.59488 
12.4051 
1.54341 
12.2755 
8.18113 
1.54341 
12.2755 
8.18113 
1.54341 
12.2755 
8.18113 
6.20871 
10.7913 
4. 

6.20871 
10.7913 
4. 

6.20871 
10.7913 
4. 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =5 and m s = |. 



No 




Ei 


genstate 


Energy 


Example 


192 


15 


i 

2> 


4J 1 1/ 






_i + 2C/+ 16 


9 


193 


|5, 


1 

2 ' 


¥.r 2 ) 






4 + 3t + U+ 18W 


8. 


194 


5, 


1 

2 ' 


3 r \ 

4 5 o,l/ 






J 2 +t+ 3 " +17W 


4.59488 


195 


|5, 


_L 

2 ' 


3 p \ 
4 ' 0,1/ 






j+t+^ + 17W 


12.4051 


196 


|5, 


_L 

2 ' 


3 p \ 
4) 1 3,2/ 






=^-- {+t+^ + 17W 


4.59488 


197 


|5, 


_L 

2 ' 


3 p \ 

4 > L 3,2/ 






^ - i + t + ^f + 


12.4051 


198 


|5, 


_L 

2 ' 


3 p \ 
Ji -L 4 1/ 




— J 
3 


j. 1 5(7 , 50 W 2V^ecos(e 5 ) 
' 3 ' 3 3 


1.54341 


199 


|5, 


1 

2 ' 


3 p \ 
4' 1 4,1/ 


-,/ 
3 


-t + 


5(7 1 50 1 V^e cos(6 5 ) , v^i" sin(6 5 ) 

3 3 3 ■ /q 


12.2755 


200 


|5, 


1 

2 ' 


3 p \ 
7,14 1 / 

4' 


—,/ 

3 


-t + 


5(7 | 501V | V^"cos(0 5 ) sin(e 5 ) 
3 3 1 3 ^3 


8.18113 


201 


5 ; 


1 

2' 


3 p \ 

4 > 4,2/ 




-,/ 
3 


j. 1 5(7 | 50 W 2V^cos(e 5 ) 
' 3 ' 3 3 


1.54341 


202 


5 5 


1 

2 ' 


3 p \ 

4 ' 4,2/ 


-.7 
3 


- t + 


5U 1 501V 1 VAe cos(6 5 ) , sin(6 5 ) 
3 1 3 3 


1 2 2755 


zUo 


1 K 
\0, 


1 

2 ' 


3 p \ 
4! 1 4,2/ 


— J 
3 




5(7 I 50 W I cos(6 5 ) sin(0 5 ) 
3^3 1 3 ^3 


o.lollo 




1 ^ 


1 

2 • 


41 1 4,3/ 




— J 
3 


4. 1 5(7 1 50 W 2VA^cos(e 5 ) 
1 " r 3 " r 3 3 


1 .o L ±o L ±i. 


205 


|5, 


\ 

2 • 


3 \ 

4! 1 4,3/ 


— J 
3 


- t + 


5(7 1 50 W 1 V^6"cos(e 5 ) | \/^6 sin(6 5 ) 
3 3 3 _^ 3 


12.2755 


206 


|5, 


! 

2 ' 


3 T 1 \ 

4' 1 4,3/ 


— J 
3 


- t + 


5(7 1 50 W I v^e'cos(05) v^6 sin(65) 
3 3 1 3 ^ 


8.18113 


207 


|5, 


1 


3^5,1) 






=^-i + t + ^ + nw 


6.20871 


2 ' 








208 


|5, 


2 ' 


41^5,1) 






+ +t+^ + 17W 


10.7913 


209 


|5, 


1 

2 ' 


f,r 5il ) 






± - t + U + 18W 


4. 


210 


|5, 


1 

2 ' 


3 p \ 

4! 1 5,2/ 






=^ - j +t+^f + 17 W 


6.20871 


211 


|5, 


1 

2 ' 


3 p \ 
4 : 1 5,2/ 






j+t+^ + 17W 


10.7913 


212 


|5, 


1 

2' 


f,r 5 , 2 > 






{ -t + U+ 18W 


4. 


213 


|5, 


1 

2 ' 


3 p \ 

4' 1 5,3/ 






j+t+^ + 17W 


6.20871 


214 


|5, 


1 

2 ' 


3 p \ 






^-±+t+^f + 17W 


10.7913 


215 


|5, 


1 

2 ' 


— , F 5,3) 






% -t + U+ 18W 


4. 



Eigenkets and eigenvalues for N e =5 and m s = |. 



No 


Eigenstate 


Energy 


Example 


216 


|5,|,f,r 2 > 


^ + 3t + U+ 18W 


8. 


217 


5) 2> T'^m) 


^- -t + U + 18W 


4. 


218 


5) 2' 


^--t + u + mw 


4. 


219 


5, 2' T'^5,3) 


If -t + U + 18 W 


4. 



adp header will be provided by the publisher 



41 



Eigenkets and eigenvalues for N e =6 and m s = — 1. 



No 


Eigenstate 


Energy 


Example 


220 


|6,-l,2,r 4 ,i) 


± - 2t + 2U + 26W 




8. 


991 


|6,-l,2,r 4>2 ) 


4 -2t + 2U + 26W 




8. 


222 


|6,-l,2,r 4>3 ) 


4 -2t + 2U + 26W 




8. 


223 


|6,-i,2,r 5 ,i) 


j + 2t + 2U + 26W 




12. 


224 


|6,-i,2,r 5;2 ) 


j + 2t + 2U + 26W 




12. 


225 


|6,-i,2,r 5 , 3 > 


4 + 2t + 2U + 26W 




12. 




Eigenkets and eigenvalues for N e =6 and m s = 


0. 


No 


Eigenstate 


Energy 


Example 


226 




6,o,o,ri) 


- v 2 As J 2 2t+ 5 !/ +25W 




4.82109 


227 




6,0,0,^) 


^J*- j-2t+^f + 25W 




16.1789 


228 




6,o,o,r 3 ,i) 


-J + 2t + 2U + 26W 




12. 


229 




6,o,o,r 3 , 2 ) 


-J + 2t + 2U + 26W 




12. 


230 




6,o,o,r 4 ,i) 






9.29844 


231 




6,o,o,r 4 ,i) 


^-j + ^-+25W 




15.7016 


232 




6,0,2,r 4 ,i) 


-2t + 2U + 26W 




8. 


233 




6,o,o,r 4 , 2 > 


=^-{ + ^ + 25W 




9.29844 


234 




6,o,o,r 4 , 2 ) 


±+ r ^ + 25W 




15.7016 


235 




6,0,2,r 4 , 2 ) 


-2t + 2U + 26W 




8. 


236 




6,o,o,r 4 , 3 > 


+ + ^ + 25W 




9.29844 


16 ( 




6,o,o,r 4 , 3 > 


^¥ - i+ r± T + 25^ 




15.7016 






6,0,2,r 4 , 3 > 


-2t + 2U + 26W 




8. 






6,0,2,r 5 ,i) 


2t + 2U + 26W 




12. 






6,o,2,r 5 , 2 > 


2t + 2U + 26W 




12. 


OA 1 




6,0,2,r B , 3 > 


2t + 2U + 26W 




12. 




Eigenkets and eigenvalues for N e =6 and m s = 


1. 


No 


Eigenstate 


Energy 


Example 


242 


|6,i,2,r 4 ,i) 


£ -2t + 2U + 26W 




8. 


243 


|6,i,2,r 4 , 2 ) 


4 - 2t + 2U + 26W 




8. 


244 


|6,i,2,r 4 , 3 > 


j - 2t + 2U + 26W 




8. 


245 


|6,i,2,r 5 ,i) 


j + 2t + 2U + 26W 




12. 


246 


|6,i,2,r 5 , 2 ) 


j + 2t + 2U + 26W 




12. 


247 


|6,i,2,r B , 3 > 


j + 2t + 2U + 26W 




12. 
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Schumann: Extended Hubbard model on small clusters 



Eigenkets and eigenvalues for N e =7 and m s = — |. 



No 


Eigenstate 


Energy 


Example 




7,-5, 


-3t + 3U + 36W 




1 9 




I 7 , f,r 4i i) 


t + 3U + 36W 




10. 


9^0 


7, -5, 4, r 4i2 ) 


t + 3U + 36W 




1 fi 

±U. 


251 


I 7 ' - §! |> r 4,3) 


t + 3U + 36W 




16. 




Eigenkets and eigenvalues for N e =7 and m s = 


1 

2 




No 


Eigenstate 


Energy 


Example 


252 


7, 5, 


-3t + 3U + 36W 




12. 


253 


7, 5, 4^4,1} 


t + 3[/ + 36T^ 




16. 


254 


|7, 5, f ,r 4i2 ) 


i + 3 £7 + 361^ 




16. 


255 


|7, ^, |,r 4j3 ) 


t + 3 J7 + 36 W 




16. 




Eigenkets and eigenvalues for N e =8 and m s = 


0. 


No 


Eigenstate 


Energy 


Example 


256 


18,0,0,^) 


4 [/ + 48 W 


20. 
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B.2 List of abbreviations 



A x = J 2 + Wt 2 + 2J (U-2W) + (U-2W) 2 

A 2 = J 2 +A8t 2 + 2J (U-2W) + (U-2W) 2 

A 3 = J 2 +6At 2 + 2J (U -2W) + (U -2W) 2 

A 4 = J 2 + 16t 2 + J (-At + 2U — A W) — At (U — 2 W) + (U — 2 W) 2 

A 5 = J 2 + A8t 2 + J (-3t + 2U — A W) — 'it (U — 2 W) + (U — 2 W) 2 

A 6 = J 2 + A8t 2 + J (3t + 2U - A W) + 'it (U - 2 W) + (U — 2 W) 2 

A 7 = J 2 + 16t 2 + At (U -2W) + (U -2W) 2 + 2J (2t + U -2W) 

A s = J 2 + 6At 2 + 8t (U -2W) + (U -2W) 2 + 2J {At + U -2W) 

A 9 = A (I2t 2 -AtU + U (J-2W)) + {-J + At + U + 2W) 2 

A w = J -2U -UW + ^/A 2 ~ cos(0 2 ) 

An = J + 4J7 + 64T^-4v/]4^cos(e3) 

A 12 = J -2U -32W + 2^/M cos(9 3 ) 

A 13 = -2J + AU + 28W+ yfM (cos(9 2 ) - V3 sin(6 2 )) 



A 20 = -J + 3t + 2U +1AW + \[M (cos(e 4 ) - V3 sin(9 4 ) 
A 21 = -J + 3t + 2U +1AW + a/A^ (cos(9 4 ) + V3 sin(6 4 ) 
A 22 = J + 3t-5U -50W + 2 yf~A % cos(9 5 ) 



A 23 = J + 3t-5U -50W + ( - cos(6 5 ) + \/3 sin(6 5 ) 



A 2i = J + 3t-5U-50W-^/% (cos{e 5 ) + V3sin(Q 5 ) 





(18) 
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Schumann: Extended Hubbard model on small clusters 



6! 



e 2 = 



/ 12 v / 3t 2 (J+U-2W)- 

arccos — — g '- 



, ( (j 2 -36t 2 +2.J (U-2W) + (U-2W) 2 ) (J+t/-2W)V 

arccos(- ^ 3 '- , 



, (j 2 -36t 2 +2J (U-2 W) + (U-2 Wf) (J+U-2 W) . 

arccos^ ^ j 

e 3 = ^ 



, (2.7 2 +36 t 2 + J (-9t+4U-8W)-9t(U~2W)+2 (U-2W) 2 ) (J+U-2 W) . 

arccos(- 3 ) 



e 4 2A 55 



arccos( 



3 

(2 ,7 2 +36 t 2 + J (9 t+4 (7-8 W)+9 t ((7-2 VK)+2 (U-2 W) 2 ) (J+U-2 W) , 



° 5 - 3 

B.3 The eigenvectors of the tetrahedron 

In the following the unnormaliscd eigenvectors are given in abbreviated form, whereby the indices 
give the number of the state in accordance with the numbering of the energy levels given above. 
For getting the normalised eigenvectors the coefficients have to be divided by the normalisation 
constant Ni, which is given for the vectors containing more than one coefficient. 

Unnormalized eigenvectors for N e = and m s = 0. 

- 10,0,0,^) 
= 1 (|0000)) 



Unnormalized eigenvectors for N e = 1 and m s = — | 



2 ' 



l*> 2 = Ii.-|.|.ri> 



i (|000d) + |00d0) + |0d00) + |d000)) 



l*> 3 = |i,-^,r 4 ,i> 



i (|000d) + |00d0) - |0d00) - |d000» 



l*> 4 = |i,-^,r 4;2 ) 



1 (|000d) - |00d0) - |0d00) + |d000)) 



2 



l*> 5 = |i,-^,r 4i3 ) 



i (|000d) - |00d0) + |0d00) - |d000)) 
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Unnormalized eigenvectors for N e = 1 and m s = ^. 

l*> 6 = IL^.ri) 

= - (|000u) + |00w0) + |0u00> + |uOOO)) 

= ^(\000u) + \00u0)-\0u00)-\u000)) 

l*> 8 - |i,^,r 4 , 2 > 

= i (|000u) - |00m0) - |0u00) + |u000)) 

l*> 9 = |i,|,|,r 4 , 3 > 

= i (|000u) - |00m0) + |0u00) - |uOOO}) 

Unnormalized eigenvectors for N e = 2 and m s = — 1. 

|*)io = |2,-l,2,r 4 ,i) 

= i (|0d0d) + \0dd0) + \d00d) + |dOdO}) 

|*)n - |2,-i,2,r 4 ,2) 

= i (\00dd) + \0d0d) - \dOdO) - \dd00)) 

= |2,-i,2,r 4 ,3) 

= i (\00dd) - \0dd0) + \d00d) + |dd00)) 

|*)i 3 = |2,-i,2,r 6 ,i) 

= i (|00dd) + |0dd0) - |d00d) + |dd00)) 

|*)i 4 = |2,-i,2,r 6 ,2> 

= i (|OdOd) - |0dd0) - |d00d) + |dOdO>) 

|*)i 8 = |2,-i,2,r 6 ,3> 

= i (|00dd) - |0d0d) + |dOdO) - |dd00)) 
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Schumann: Extended Hubbard model on small clusters 



Unnormalized eigenvectors for N e = 2 and m s = 



0. 



|*)i 6 = |2,o,o,ri) 

= Ci 6 _i (|0002) + |0020) + |0200) + |2000» 

+ C*i6_ 2 (|00du) - \00ud) + \0d0u) + \0du0) - \0u0d) - \0ud0) 

+ \d00u) + \d0u0) + \du00) - \u00d) - \u0d0) - \ud00)) 

C w -i = Vst 

C 16 -2 - jj= + J-4:t + U-2W~) 

N\e — 2 \J G\s-i 2 + 3 C16-2 2 

|*)i 7 = |2,o,o,ro 

= Cvt-1 (|0002) + |0020) + |0200) + |2000» 

+ Cvt-2 (|00du) - \00ud) + \0d0u) + \0du0) - \0u0d) - \0ud0) 

+ \d00u) + \d0u0) + \du00) - \u00d) - \u0d0) - \ud00)) 

C17-1 = Vst 

Ci7-2 = (y/M- J + 4,t-U + 2W^j 

N 17 = 2 V / C 17 _ 1 2 + 3C 17 _ 2 2 

|*)i 8 - |2,o,o,r 3 ,i) 

= C 1S -i (\00du) - \00ud) + \0du0) - \0ud0) + \d00u) + \du00) - \u00d) - \ud00}) 

+ C 1S -2 (\0d0u) - \0u0d) + \d0u0) - \u0d0)) 



2V6 



1 



N 1S = 



V6 

2 \J 2 Cis-i 2 + Cis-2 2 



|*>19 = 



2,o,o,r 3 , 2 ) 

(|00du) - \00ud) - \0du0) + \0ud0) - \d00u) + \du00) + \u00d) - \ud00)) 



|*>20 



|2,o,o,r 4 ,i) 

C20-1 (|0002) + |0020) - |0200) - |2000)) 
C20-2 (|00du> - |00ud) - \du00) + \ud00)) 



+ 




N 2 o = 



2 \J C20-1 2 + C20-2 2 
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|*) 2 i = |2,o,o,r 4 ,i) 

= C21-1 (|0002) + |0020) - |0200) - |2000)) 
+ C21-2 (|00du) - \00ud) - \du00) + \ud00)) 



C21-1 
C21-2 

N21 



= -t 



2 \/ C*2i— i 2 + C21-2 2 



|*)22 = |2,0,2,r 4 ,i) 

= — ^ (|0d0u) + |0dw0) + |0u0d) + |0ud0) + \d00u) + |dOuO) + \u00d) + \uOdO}) 
2 v 2 

|*)23 = |2,o,o,r 4 , 2 ) 

= C23-1 (|0002) - |0020) - |0200) + |2000)) 
+ C23-2 (|0du0) - \0ud0) - \dOOu) + \u00d)) 



C23-1 
C23-2 

^23 



t 
1 

4 



r M-J-U + 2W^j 



2 \/ C*23-i 2 + C23-2 2 



|*)24 - |2,o,o,r 4 , 2 ) 

= C* 24 -i (|0002) - |0020) - |0200) + |2000)) 
+ C* 24 - 2 (|OduO> - |OudO) - |d00u) + |u00d» 



C*2 4 -l — 

C 2 4-2 = 

iV 24 = 

I*) 25 = 



1 (v 7 ^- J-U + 2W^j 

2 \J C24-1 2 + C 24 -2 2 



|2,0,2,r 4 ,2) 

-^-j= (|OOdu) + |OOud) + |0d0u) + |OuOd) - |dOuO) - |dwOO) - |uOdO) - \udOO)) 



|*)26 = |2,o,o,r 4 , 3 > 

= C 26 -i (|0002) - |0020) + |0200) - |2000» 
+ C 26 -2 (|OdOu) - \OuOd) - \dOuO) + \uOdO)) 



C*26-l — 
C 2 6-2 = 

^26 = 



l - (-^/M- J-U + 2W) 
2 \J C 2 6-i 2 + C 2 6-2 2 
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Schumann: Extended Hubbard model on small clusters 



|*) 27 = |2,o,o,r 4 , 3 > 

= C 27 -i (|0002> - |0020) + |0200) - |2000)) 
+ C27-2 (|0rf0u) - \0u0d) - \d0u0) + \u0d0}) 

C27-1 = —t 

C 27 -2 = ^(-/Al- J -U + 2W^ 
N27 = 2 \J C27-1 2 + C27-2 2 

|*)28 - |2,o,2,r 4 ,3> 

= — !-= (|00du) + |00ud) - \0du0) - \0ud0) + \d00u) + \du00) + \u00d) + \ud00)) 
2 V2 

|*)29 = |2,0,2,r 6l i) 

= — !-= (|00du) + \00ud) + \0du0) + \0ud0) - \d00u) + \du00) - \u00d) + \ud00}) 
2 v 2 

|*)so - |2,0,2,r 6 , 2 > 

= — ^ (|0d0u) - |0dw0) + |0u0d) - |0ud0) - |d00u) + |dOuO) - |u00d) + \uOdO)) 
2 v 2 

(|00du) + |00wd) - |0rf0u) - \0u0d) + \dOuO) - \du00) + \uOdO) - \ud00)) 



|*>3i - |2,0,2,r 6l 3 

1 



2^ 

Unnormalized eigenvectors for N e = 2 and m s = 1. 



|*)32 = |2,i,2,r 4 ,i) 



i (\0u0u) + |0uw0) + \u00u) + \u0u0)) 



|*)33 - |2,l,2,r 4 ,2) 



= i (|00mi) + |0u0u) - |m0m0) - |mi00)) 

|*>34 = |2,l,2,r 4 ,3> 

= i (\00uu) - |0uw0) + |u00u) + |tiu00» 

|*)38 = |2,i,2,r 5jl > 

= i (|00mi) + |0uw0) - |u00u) + |mi00)) 

|*)36 = |2,i,2,r 5 , 2 > 

= i (\0u0u) - \0uu0) - \u00u) + \uOuO)) 

|*) 37 = |2,i,2,r 6 ,3> 

= i (\00uu) - \0u0u) + \uOuO) - \uu00)) 
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Unnormalized eigenvectors for N e = 3 and m s = — |. 

I*>38 = |3,-|^,r 2 ) 

= i (\0ddd) - \dOdd) + \ddOd) - \dddO)) 

l*>39 = |3,-|^,r 5il ) 

= 1(| Oddd) + | rfOdd) + | ddOd) + \ dddO) ) 

l*>40 = |3,-|^,r 5 , 2 ) 

= 1 (|0d<ta) - |dOdd) - \ddOd) + \dddO)) 

l*>4i = |3,-|,^,r 5;3 ) 

= 1 (|0d<ta) + |dOdd) - \ddOd) - \dddO)) 
Unnormalized eigenvectors for N e = 3 and m s = — |. 

I*>42 = |3,-|,|,ri> 

= — !-= (|002d) + |00rf2> + |020rf> + |02d0) + |0d02) + |0rf20) 
2 V3 

+ |200rf) + |20d0) + |2d00) + |rf002) + |d020) + |rf200>) 



l*)43 = |3,-l^,r 2 ) 

= — !-= (\0ddu) + \0dud) + \0udd) - \dOdu) - \dOud) + \ddOu) 
2 v3 

-\dduO) + \duOd) - \dudO) - \uOdd) + \udOd) - \uddO)) 

l*>44 - |3,-^,r 3 ,i) 

= C44-1 (|002d) + |00d2) + |02d0) + |0d20) + |200rf) + |2d00) + |d002) + |d200>) 

+ C44-2 (|020d) + |0rf02) + |20d0) + |d020» 

+ C44-3 (\0ddu) - \0udd) - \dOud) - \ddOu) + \dduO) + \duOd) + \uOdd) - \uddO}) 

6*44-1 



2V2 

C 44 _ 2 = ± 

C44-3 = -^-j=(^/M + J-2t + U-2W^j 

N44 = 2 \J 2 C44_i 2 + C44_ 2 2 + 2 C44_3 2 
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Schumann: Extended Hubbard model on small clusters 



l*>« - |3,-^,r 3 ,i) 

= C45-1 (|002d) + |00d2) + |02d0) + |0d20) + |200rf) + |2d00) + |d002) + |d200>) 

+ C45-2 (|020d) + |0rf02) + |20d0) + |d020)) 

+ ^45-3 (\0ddu) - \0udd) - \d0ud) - \dd0u) + \ddu0) + \du0d) + \u0dd) - \udd0}) 



C45-1 
C45-2 
C45-3 
^45 



-t 
t 

-1 
4V2 



[y/M- J + 2t- u + 2vy) 

2 \j 2 C45-1 2 + C 4 5_2 2 + 2 C 45 _ 3 2 



l*>« = |3,-^,r 3 ,2> 

= C 46 -i (|002d) + |00rf2> - |02d0) - |0d20) - |200rf) + |2d00> - |d002> + |rf200>) 

+ C 46 ^2 (\0ddu) + \0udd) - \d0ud) + \dd0u) - \ddu0) + \du0d) - \u0dd) - \udd0}) 

+ C i6 - 3 (\0dud) - \d0du) - \dud0) + \ud0d}) 



C46-1 
C46-2 

C46-3 

^46 



4\/fi V 



'A 4 + J-2t + U-2W 



) 



4V6 

— \= (JaI+J -2t + U -2W) 
2V6 V / 

2 \J 2 C46-1 2 + 2 C46-2 2 + C46-3 2 



I*) 



47 



1 3 

2'4' 



|3,--,-,r 3 , 2 ) 



= C47-1 (|002d) + |00d2) - |02d0) - |0d20) - |200d) + |2d00) - |d002) + |d200)) 
+ C47-2 (\0ddu) + \0udd) - \d0ud) + \dd0u) - \ddu0) + \du0d) - \u0dd) - \udd0}) 
+ C47-3 (\ odud ) - \ dGdu ) - \ dud0 ) + \ud0d)) 



C47-1 
C47-2 

C47-3 

^47 



-(#) 

2 

A 



M 4 - J + 2t-£/ + 2VK 



) 



4V6 

— (jAl - J + 2t- U + 2W) 

2 \J 2 C47_i 2 + 2 C*47_2 2 + C< 



^47-3 
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|*>48 = 

+ 
+ 

C48-I = 
C48-2 = 
C*48-3 = 



I*) 



' 3 ' _ 2' I' 1 " 4 ' 1 ^ 

C 48 -i (|002d) + |00d2) - |2d00) - |rf200}) 

ds-2 (|020rf) + |02d0) - |0d02) - |0d20) + |200rf) + |20d0) - |d002) - |d020)) 
C 48 -3 (\0ddu) - \0dud) + \dOdu) - \dOud) - \duOd) - \dudO) + \udOd) + \uddO)) 



-t 



3V2 
-t 

2\/2 
1 

18 V2 



J +12t + U -2W + 2 VA5 cos 



(©4)) 

(j - At + U -2W + 2 cos(9 4 )) 
(A 19 2 + 3 Jt - 45 t 2 + 6 J U- 15 1 U - 3 U 2 + 24 J W - 78 1 W - 60 U W 

-192 W 2 + 6 v^i cos(6 4 ) 
+ 12 y/A^U cos(6 4 ) + 48 y/%W cos(9 4 ) 



N 48 = 

|*>49 = 

+ 
+ 

6*49-1 = 
C49-2 = 
C49-3 = 

^49 = 



2 ^/ C48-1 2 + 2 (C48-2 2 + C48-3 2 ) 
' 3 ' ~2~' 4' r4 '^ 

C49-1 (|002d) + |00rf2> - |2d00) - |d200» 

C49-2 (|020d) + |02rf0) - |0d02) - |0d20) + |200rf) + |20d0> - |d002) - |d020>) 
C49-3 (\0ddu) - \0dud) + \d0du) - \d0ud) - \du0d) - \dud0) + \ud0d) + \udd0)) 



-t 

3V2 

-t 

2V2 
1 

18 V2 



(j +12t + U -2W - ^/A~ 5 cos(e 4 ) - VS \[A~ h sin(9 4 )) 
(j -At + U -2W - v / I^cos(e 4 ) - \/3~V^5 sin(6 4 )) 

(A 2 i 2 - z A 2 it - mt 2 - & A21U 

+ 9tU + 9U 2 -2AA 21 W + 36tW + 72UW + 1AAW 2 ) 
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+ 

C50-I = 
C5O-2 = 
C5O-3 = 

^50 = 



2 \J C'49-i 2 + 2 (C49_ 2 2 + C49_ 3 2 ) 

C50-1 (|002d) + |00rf2) - |2d00) - |d200>) 

C50-2 (|020d) + |02d0) - |0d02) - |0d20) + |200rf) + |20d0) - |d002) - |d020» 
C50-3 (\0ddu) - \0dud) + \d0du) - \dOud) - \du0d) - \dud0) + \udOd) + \udd0)) 



-t 

3V2 
-t 

2y/2 ' 
1 

18 V2 



(j +12t+U -2W - y/M cos(e 4 ) + V3 y/M sin(6 4 )) 

J -At + U -2W - yfM cos(9 4 ) + V3 V / Asin(e 4 )) 

(A 20 2 - 3 A20 i - 36 1 2 - 6 A 20 C/ 
+ 9tC/ + 9 U 2 -2AA 20 T^ + 36t W + 72 UW + 144 W 2 ) 



2 C50-1 2 + 2 (C50-2 2 + C50-3 2 ) 
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Schumann: Extended Hubbard model on small clusters 



I*>bi = |3,-^,r 4 , 2 ) 

= C51-1 (|002d) - |00d2) + |020d) - |0d02) - |20rf0) - |2d00) + |d020) + |d200» 

+ C51-2 (|02d0) + |0rf20) - |200d) - |d002)) 

+ C51-3 (\0dud) - \0udd) + \dOdu) + \ddOu) + \dduO) - \dudO) - \uOdd) - |udOd» 
C51-1 = [j - 4t + U -2W + 2 \f~A 5 cos(e 4 )) 



2V2 

C51-2 - (.J +12t + U -2W + 2^fM cos(e 4 )) 



3\/2 

C51-3 = — K= (^i9 2 + 3 Jt-45t 2 + 6 JC/- 15t ?7 - 3 C/ 2 + 24 JW - 78tW - 60UW 
18 v 2 



-192 VF 2 + 6 v^i cos(6 4 ) 

+ 12 U cos(9 4 ) + 48 W cos(6 4 ) 

JV51 - 2 v /2C 5 i-i 2 + C5i-2 2 + 2C 5 i 



2 

-3 



l*>B2 = |3,-^,r 4i2 ) 

= C52-1 (|002d) - |00rf2> + |020d) - |0d02) - |20rf0) - |2d00) + |d020) + |d200)) 

+ C52-2 (|02d0) + |0d20) - |200d) - |d002)) 

+ C52-3 (\0dud) - \0udd) + \dOdu) + \ddOu) + \dduO) - \dudO) - \uOdd) - \udOd)) 

C52-1 = ^= (j-4t + U -2W - \/A~ 5 cos(9 4 ) - V3 \/T 5 sin(9 4 )) 

C52-2 = -A= (j +12t+U -2W - y/M cos(9 4 ) - V3 sin(6 4 )) 



3V2 

C 52 - 3 - (A 21 2 -SA 21 t-36t 2 -6A 21 U 

+ 9tU + 9U 2 -24A 21 W + 36tW + 72UW + U4W 2 ) 



N 52 = 2^2C52-l 2 +C*52-2 2 + 2C 52 - 



2 

3 



l*>53 = |3,-^,r 4 , 2 ) 

= C53-1 (|002d) - |00rf2> + |020d) - |0d02) - |20rf0) - |2d00> + |d020> + |d200)) 

+ C53-2 (|02d0) + |0rf20> - |200d) - |d002)) 

+ C53-3 (\0dud) - \0udd) + \dOdu) + \ddOu) + \dduO) - \dudO) - \uOdd) - \udOd)) 

C53-1 = ^= {j -At + U -2W - x/T 5 cos(6 4 ) + V3 x/% sin(9 4 )) 

C53-2 = ^= ( J + 12 1 + U - 2 W - jA~ h cos(e 4 ) + V3 y/M sin(6 4 )) 

C53-3 = — V ( A 2o 2 -3A 20 t-36t 2 -6A 2a U 
18 V2 

+ 9tU + 9U 2 -24A 20 W + 36tW + 72UW + U4W 2 ) 



N53 — 2 \J 2 C53_i 2 + 653-2 2 + 2 C53_3^ 
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l*>54 = 13,-^,^4,3) 

= C54-1 (|002rf) - |00d2) - |02d0) + |0d20) + |200rf) + |2d00) - |d002) - |d200>) 

+ C54-2 (|020rf) + |0rf02) - |20d0) - |rf020>) 

+ ^54-3 (\0ddu) - \0udd) + \dOud) - \ddOu) - \dduO) + \duOd) - \uOdd) + \uddO}) 

C54-1 = -^j= [j-4t + U-2W + 2 V^cos(6 4 )) 

C54-2 = -^= (j+12t + U -2W + 2^ cos(e 4 )) 



3\/2 

C54-3 = {A w 2 + ZJt-Abt 2 + &JU -IbtU -ZU 2 + 2AJW -7%tW -&QUW 

18 v 2 



-192 W 2 + 6 v^i cos(6 4 ) 

+ 12 U cos(9 4 ) + 48 W cos(6 4 ) 

iV 54 = 2 \J 2 C54-1 2 + C54-2 2 + 2 C54 



2 

-3 



l*>55 = |3,-^,r 4l3 ) 

= C55-1 (|002d) - |00rf2> - |02d0) + |0d20) + |200rf) + |2d00) - |d002> - |rf200>) 

+ C55-2 (|020d) + |0rf02> - |20d0) - |rf020» 

+ C55-3 (\0ddu) - \0udd) + \dOud) - \ddOu) - \dduO) + \duOd) - \uOdd) + \uddO}) 

C55-1 = ^= (j-4t + t/-2W - cos(9 4 ) - V3 sin(9 4 )) 

C55-2 = ^= (j+12i+{/-2W- \[M cos(9 4 ) sin(6 4 )) 

C55-3 = (A 21 2 -SA 21 t-36t 2 -6A 21 U 

+ 9tU + 9U 2 -24A 21 W + 36tW + 72UW + U4W 2 ) 



-1 + C55-2 +2C55- 



2 

3 



Was = |3,-^,r 4 , 3 ) 

= C 56 -i (|002d) - |00rf2> - |02d0) + |0d20) + |200rf) + |2d00> - |d002> - |d200>) 

+ C56-2 (|020d) + |0rf02> - |20d0) - |rf020» 

+ C 56 -3 (\0ddu) - \0udd) + \dOud) - \ddOu) - \dduO) + \duOd) - \uOdd) + \uddO)) 

C 56 -i = (j -At + U -2W - cos(6 4 ) + V3 y>% sin(6 4 )) 

C 56 -2 = ^= (j +12t+U -2W - ^TAr, cos(e 4 ) + V3 sin(6 4 )) 

C 56 -3 = —^(A 2a 2 -3A 20 t-36t 2 -6A 2a U 
18 V2 

+ 9tU + 9U 2 -24A 20 W + 36tW + 72UW + U4W 2 ) 



^56 — 2 ^/ 2 C56-1 2 + C*56-2 2 + 2 C56-3^ 
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Schumann: Extended Hubbard model on small clusters 



l*>57 = |3,-|,|,r 5;1 > 

= C57-1 (|002d) - |00d2) + |02d0) - |0d20) - |200d) + |2d00) + |d002) - |rf200>) 

+ C57-2 (|0ddu) + \0udd) + \dOud) + \ddOu) + \dduO) + \duOd) + \uOdd) + \uddO)) 

+ C57-3 (\0dud) + \dOdu) + \dudO) + \udOd)) 



C: 



57-1 — 



2 



C57-2 = -^-j=(^/T 7 + J + 2t + U -2W^j 
C57-3 = — i [y/T 7 + J + 2t + U -2W^j 



N§ 7 — 2 \J 2 C57-1 2 + 2 C*57_2 2 + C57- 



2 

3 



l*>58 - |3,-^,r B ,i) 

= C 58 -i (|002d) - |00d2) + |02d0) - |0d20) - |200d) + |2d00) + |d002) - |d200>) 

+ C58-2 (|0ddu) + |0udd) + |dOud) + |ddOu) + |ddwO) + \duOd) + |uOdd) + |uddO)) 

+ ^58-3 (|0dud) + |dOdu) + \dudO) + \udOd)) 



C: 



58-1 



m 

2 

-1 



C 58 _ 2 = j-j= (y/Ai- J-2t-U + 2w) 



C: 



58-3 — 



l -= [y/T 7 - J - 2t - U + 2W^ 



2V6 



^58 — 2a/2C58-1 +2C 58 _2 + C 58 _ 3 z 



1 1 ^ 

l*>69 = 13,-2^,^,0 



— ^= (\0ddu) + \0dud) + \0udd) + \d0du) + \d0ud) + \dd0u) 
2 y 3 

+ \dduQ) + \du0d) + \dud0) + \u0dd) + \ud0d) + \udd0)) 



l*) 60 = |3,-^,r 5 , 2 ) 

= C 6 o-i (|020d) - |02d0) - |0d02) + |0d20) - |200d) + |20d0) + |d002) - |d020» 

+ C 60 -2 (|0ddu) + |0dud) - |d0du) - |d0ud) - |du0d) + |dud0) - \ud0d) + \udd0}) 

+ C 60 _3 (|0udd) - |dd0u) + |ddu0) - \u0dd}) 
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£-60-1 
£-60-2 

C60-3 



-(#) 

2 

(V^r + J + 2t + U -2W) 
[^/Tt + J + 2t + U -2W^j 



+ 2 6*60-2 + ^60- 



2 

3 



l*>ei - |3,-^,r 5 , 2 ) 

= C 6 i_i (|020d) - |02d0) - |0d02) + |0d20) - |200rf) + |20d0) + |d002) - |d020)) 

+ C 6 i_2 (\0ddu) + \0dud) - \d0du) - \dOud) - \duOd) + \dud0) - \udOd) + \udd0)) 

+ C 6 i-3 (\0udd) - \dd0u) + \dduO) - \u0dd)) 



Cqi-1 — 
£-61-2 = 

a 



-(#) 



2 

-±=(VA 7 -J-2t-U + 2w) 
61-3 = (VT 7 -J-2t-U + 2\¥) 



N 61 = 2V2C 61 _i" + 2C 6 i_ 2 ^ + C 6 i_ 



l*> 6 2 - |3,-l^,r 5>2 ) 



— !-= (|0dc?M) + |0du<i) + \0udd) - \dOdu) - \dOud) - \ddOu) 
2 V3 

+ \dduO) - \du(H) + \dudO) - \uOdd) - \udOd) + \uddO)) 



l*>63 = |3,-|,|,r 5 , 3 ) 

= C 63 _i (|002rf> - |00rf2> - |020d) + |0d02) + |20rf0) - |2d00> - |d020> + |rf200>) 

+ C 63 - 2 (\0ddu) + \dOud) - \duOd) - \uddO)) 

+ C" 63 _ 3 (\0dud) + \0udd) + \dOdu) - \ddOu) - \dduO) - \dudO) + \uOdd) - \udOd)) 



C63-1 
C63-2 

C63-3 



— \= ( + J + 2t + U -2W) 
2^6 V v J 

(V^t + J + 2t + U -2W^j 

^63 = 2 \J 2 C63-1 2 + C63-2 2 + 2 C63-3 2 
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Schumann: Extended Hubbard model on small clusters 



<>4 



C&4-1 

C64-3 
^64 



65 



= I 3 ; ~2> 4' r 5,3) 

= C 64 _i (|002d) - |00d2) - |020d) + |0d02) + |20rf0) - |2d00) - |d020) + |d200» 

+ C 64 -2 (\0ddu) + \dOud) - \duOd) - \uddO}) 

+ C 64 -3 (\0dud) + \0udd) + \dOdu) - \ddOu) - \dduO) - \dudO) + \uOdd) - \udOd)) 

_ -_m 

2 

= — !-= (^/T 7 - J -2t-U + 2w\ 
2V6 V / 

= — \= (-JAt- J~2t-U + 2W) 
4V6 V v J 

= 2 \J 2 C64-1 2 + C64-2 2 + 2 Cqh-Z 2 

1 15 

1 



2 V3 



(\0ddu) + \0dud) + \0udd) + \dOdu) + \dOud) - \ddOu) 



-\dduO) - \duOd) - \dudO) + \uOdd) - \udOd) - \uddO)) 



l*>e 



I*) 



<>7 



I*) 



(>8 



Unnormalized eigenvectors for N e = 3 and m s 



|3,i,2,r.) 



2 ' 



1 



2 V3 



(|002u) + |00u2) + |020u) + |02u0) + |0u02) + |0u20) 



+ |200u> + |20u0) + |2w00) + |u002) + |u020> + |u200)) 
(\0duu) + \0udu) + \0uud) - \dOuu) + \duOu) - \duuO) 



^68-1 
C*68-2 
C*68-3 
^68 



2 V3 

— \uOdu) — \uOud) + \udOu) — \uduO) + \uuOd) — \uudO)) 
13 

= C 68 ^i (|002u) + |00u2) + |02u0) + |0u20) + |200u) + |2u00) + |u002) + |u200)) 
+ C 68 -2 (|020u> + |0u02) + |20u0) + |u020)) 

+ Ces-3 (\0duu) - \0uud) - \dOuu) + \duuO) + \uOdu) - \udOu) + \uuOd) - \uudO)) 

-t 



2V2 
t 

71 



l -j= [\fM + J- 2t + U - 2W^j 



4^2 



2 \J 2 C68-1 2 + C*68-2 2 + 2 C68-3 2 
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l*>69 = |3,|,|,r 3 ,i> 

= C 69 _i (|002u) + |00u2) + |02u0) + |0w20) + |200u) + |2u00) + |u002) + |u200)) 

+ Ceg-2 (|020u) + |0u02) + |20u0) + |w020)) 

+ C 69 _ 3 (\0duu) - \0uud) - \dOuu) + \duuO) + \uOdu) - \udOu) + \uuOd) - \uudQ)) 



-t 



2V2 
t 

C'69-3 = — y (VM- J + 2t-U + 2W^j 



Neg — 2 yj 2 Ceg-i 2 + Cqv-2 2 + 2 Ceg-3^ 



I*>to = |3,|,|,r 3 , 2 > 

= C70-1 (|002u) + |00u2) - |02u0) - |0w20) - |200u> + |2w00) - |u002) + |u200)) 

+ C70-2 (\0duu) + \0uud) - \dOuu) - \duuO) - \uOdu) + \udOu) + \uuOd) - \uudO}) 

+ C70-3 (\0udu) + \duOu) - \uOud) - \uduO)) 



C70-1 
C70-2 



I* 



2 

(V^4 + J -2t + U -2W) 
C70-3 = — i (y/M + J-2t + U-2W) 



N70 — 2 \J 2 C70-1 2 + 2 C70-2 2 + C70- 



2 

3 



l*>n = |3,|,|,r 3>2 > 

= C 7 i_i (|002u) + |00u2) - |02u0) - |0u20) - |200u) + |2u00) - |u002) + |u200)) 

+ C71-2 (\0duu) + \0uud) - \dOuu) - \duuO) - \uOdu) + \udOu) + \uuOd) - \uudO)) 

+ C71-3 (\0udu) + \duOu) - \uOud) - \uduO)) 



C71-1 — 
C71-2 = 



1* 



jj= (^/M- J + 2t-U + 2W^j 
C 71 -3 = -^-j=[^/M- J + 2t-U + 2W^j 
N-ji = 2 \J 2 C71-1 2 + 2 C7i_2 2 + CVi_3 2 
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Schumann: Extended Hubbard model on small clusters 



l*>72 = |3,|,|,r 4>1 > 

= C* 72 -i (|002u) + |00u2) - |2?i00) - |u200» 

+ C72-2 (|020u) + |02u0) - |0u02) - |0w20) + |200u) + |20u0) - |u002) - |u020)) 

+ C72-3 (\0udu) - \0uud) - \duOu) - \duuO) + \uOdu) - \uOud) + \udOu) + \uduO)) 

C72-1 = — ^= (j+12t+U -2W + 2y/A^ cos(e 4 )) 



3V2 

CV2-2 = ^= (j-4t+ U-2 ^ + 2^^008(64)) 

C72-3 = — ^= (A 19 2 + 3 Ji-45t 2 + 6 JU- 15t ?7 - 3 C/ 2 + 24 JW - 78tW - 60UW 
18 v 2 

-192 W 2 + 6 v^i cos(6 4 ) 
+ 12 V^t/ cos(6 4 ) + 48 y/%W cos(9 4 ) 

N72 — 2 \J C72-1 2 + 2 (C72-2 2 + C*72-3 2 ) 

l*>73 = I3,|,|,r 4>1 > 

= C73-1 (|002u) + |00u2) - |2u00) - |u200)) 

+ C73-2 (|020u) + |02u0) - |0u02) - |0m20) + |200u) + |20m0) - |u002) - |u020)) 
+ C 7 3_ 3 (\0udu) - \0uud) - \duOu) - \duuO) + \uOdu) - \uOud) + \udOu) + \uduO)) 

C73-1 = — J= (j + 12t + U -2W - y% cos(e 4 ) - V3 y% sin(6 4 )) 



3V2 

C73-2 = — ^= ( J '- 4t + U -2W - \/M cos(9 4 ) - V3 \/M sin(6 4 )) 



2V2 

9tU + 9U 2 -24A 21 W + 36tW + 72 UW + 144 W 2 ) 



C73-3 = 7= U21 2 ~3A 2 it-36t 2 -6A21U 

18 V2 



\J C73-1 2 + 2 (C73_ 2 2 + C73_ 3 2 ) 



N 73 = 2jC> 

l*>74 = |3,|,|,r 4il > 

= C74-1 (|002u) + |00u2) - |2u00) - |u200)) 

+ C 74 - 2 (|020it) + |02u0) - |0u02) - |0u20) + |200u) + |20u0) - |u002) - |u020)) 

+ C 74 _ 3 (\0udu) - \0uud) - \duOu) - \duuO) + \uOdu) - \uOud) + \udOu) + \uduO)) 

C74-1 = — J= (j +12t+U -2W - y/M cos(e 4 ) + V3 y/M sin(6 4 )) 



3V2 

C74-2 = —4= (j -At + U -2W - y>% cos(9 4 ) + V3 sfM sin(6 4 )) 



2^2 

C74-3 = ^= (A 20 2 - 3^20^-36 1 2 - 6 A 20 C/ 

+ 9t[/ + 9C/ 2 -24A 20 W + 36iW + 72f7W + 144 W 2 ) 

^74 = 2 ^/ C74_i 2 + 2 (C7 4 _2 2 + C7 4 _3 2 ) 
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l*>75 = |3,|,|,r 4)2 ) 

= C 75 _i (|002u) - |00u2) + |020?i) - |0u02) - |20u0) - |2u00) + |u020) + |u200)) 

+ C75-2 (|02u0) + |0u20) - |200«) - |w002)) 

+ C 75 _ 3 (\0duu) - \0udu) + \dOuu) + \duOu) - \uOud) + \uduO) - \uvttd) - \uudO)) 

C75-1 = — ^= (j-4i + U-2 1^ + 2^^003(64)) 



2^2 

C75-2 - — ^= (j+12t + U -2W + 2^ cos(e 4 )) 



3\/2 

C75-3 = (A 19 2 + 3Jt-45t 2 + 6JU-15tU-3U 2 + 24JW-78tW-60UW 

18 v 2 



-192 + 6 v^i cos(6 4 ) 

+ 12 C/ cos(9 4 ) + 48 y/A^W cos(6 4 ) 

iV 75 = 2 v /2C75-i 2 + C75-2 2 + 2C 75 



2 

-3 



I*>r6 = |3,|,|,r 4>2 ) 

= C 76 -i (|002u) - |00u2) + |020u) - |0m02> - |20u0) - |2u00) + |u020) + |u200)) 

+ C 76 -2 (|02u0) + |0u20) - |200u) - |u002)) 

+ C 76 -3 (|0dtm) - \0udu) + \dOuu) + \duOu) - \uOud) + \uduO) - \uuOd) - \uudO)) 

C 76 -i = ^= (j-4t + t/-2W - cos(9 4 ) - V3 sin(9 4 )) 

C 76 -2 = — ^= (j+12i+{/-2W- cos(9 4 ) sin(0 4 )) 



3V2 

C 76 -3 - (A 21 2 -3^ 21 i-36t 2 -6A 21 C/ 

+ 9tU + 9U 2 -24A 21 W + 36tW + 72UW + 144 W 2 ) 



-1 + C76-2 +2C76- 



2 

3 



W77 = |3,|,|,r 4>2 > 

= C77-1 (|002u) - |00u2) + |020u) - |0u02> - |20u0) - |2u00) + |u020) + |u200)) 

+ CV 7 _2 (|02u0) + |0u20) - |200u) - |u002)) 

+ C77-3 (\0duu) - \0udu) + \dOuu) + \duOu) - \uOud) + \uduO) - \uuOd) - \uudO)) 

C77-1 = ^= (j -At + U -2W - cos(6 4 ) + V3 x/% sin(6 4 )) 

C77-2 = -^j= (j+12t+U-2W- v /^cos(e 4 ) + \/3v / Asin(e 4 )) 

C77-3 = — V ( A 2o 2 -3A 20 t-36t 2 -6A 2a U 
18 V2 

+ 9tU + 9U 2 -24A 20 W + 36tW + 72UW + U4W 2 ) 



Ntj ~ 2 \J 2 C 77 _i 2 + Ctj- 2 2 + 2 C77- 



2 

3 
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Schumann: Extended Hubbard model on small clusters 



l*>78 - |3,|,|,r 4 , 3 ) 

= C 78 _i (|002u) - |00u2) - |02u0) + |0w20) + |200u) + |2u00) - |u002) - |u200)) 

+ C78-2 (|020u) + |0u02) - |20u0) - |w020)) 

+ C 78 -3 (\0duu) - \0uud) + \d0uu) - \duu0) - \u0du) - \ud0u) + \uu0d) + \uud0)) 
C 78 -i = (j-4t + U-2 1^ + 2^003(64)) 

C 78 -2 = — ^= (j+12t + U -2W + 2x/A^ cos(e 4 )) 



3\/2 

C 78 - 3 = — ^= (Ai 9 2 + 3 Jt-45t 2 + 6 JU- lbtU-ZU 2 + 24 JW - 78tW - 60UW 
18 v 2 



-192 W 2 + 6 v^i cos(6 4 ) 

+ 12 U cos(9 4 ) + 48 y/A^ W cos(6 4 ) 

-/V78 = 2 ^/ 2 C78-1 2 + C 7 s-2 2 + 2 C 7 s 



2 

-3 



I*)™ = |3,|,|,r 4 , 3 > 

= C79-1 (|002u) - |00u2) - |02u0) + |0m20> + |200u> + |2u00) - |u002) - |u200)) 

+ CV9-2 (|020u) + |0u02) - |20u0) - |u020)) 

+ C79-3 (\0duu) - \0uud) + \d0uu) - |duu0) - \u0du) - \ud0u) + \uu0d) + \uud0)) 

C 79 _i = -^-j= (j -At + U -2W - cos(9 4 ) - V3 sin(9 4 )) 

C79-2 = ^= (j+12t+U-2W- \[M cos(9 4 ) -V3VA sin(6 4 )) 

C79-3 - (A 21 2 -3A 21 t-36t 2 -6A 21 U 

+ 9tU + 9U 2 -24A 21 W + 36tW + 72UW + U4W 2 ) 



N 79 = 2^2Ct~ 9 -1 + C 7 9_2 +2C79- 



2 

3 



l*>8o = I3,|,|,r 4 , 3 > 

= C 80 _i (|002u) - |00u2) - |02u0) + |0m20> + |200u) + |2u00) - |u002) - |u200)) 

+ 0*0-2 (|020it) + |0u02) - |20u0) - |w020)) 

+ C 8 o-3 (\0duu) — \0uud) + \dOuu) — \duu0) — \uOdu) — \udOu) + \uu0d) + \uud0)) 

C 80 _i = (j -At + U -2W - cos(6 4 ) + V3 sin(6 4 )) 

C 80 -2 = ^= ( J + 12 1 + U - 2 W - jA~ h cos(e 4 ) + V3 y/M sin(6 4 )) 

C 80 -3 = —^(A 20 2 -3A 20 t-36t 2 -6A 20 U 
18 V2 v 

+ 9if7 + 9 U 2 -24A 20 T^ + 36t W + 72 UW + 144 W 2 ) 



^80 — 2 \J 2 C$o-i 2 + C$o- 2 2 + 2 C 8 o-3^ 
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l*>si = I3,|,|,r 5 ,i> 

= C 8 i_i (|002u) - |00u2) + |02u0) - |0w20) - |200u) + |2u00) + |u002) - |u200» 

+ Csi-2 (\0duu) + \0uud) + \dOuu) + \duuO) + \uOdu) + \udOu) + \uuOd) + \uudO}) 

+ C 8 i-3 (\0udu) + \duOu) + \uOud) + \uduO)) 



Csi-l 
C81-2 

C81-3 



— U ( \/%+J + 2t + U-2w) 

— \= ( \[A~i + J + 2t + U -2W 
2V6 V 

2 yj 2 Csi-i 2 + 2 Csi-2 2 + C) 



81-3 



l*>82 = |3,|,|,r 5il > 

= C* 82 -i (|002u) - |00u2) + |02u0) - |0w20) - |200u) + |2u00) + |u002) - |u200)) 

+ C 8 2- 2 (\0duu) + \0uud) + \dOuu) + \duuO) + \uOdu) + \udOu) + \uuOd) + \uudO}) 

+ C 82 -3 (\0udu) + \duOu) + \uOud) + \uduO)) 



C-S2-1 
C%2-2 

N 82 



2 

-1 



_ (JI^ - J -2t-U + 2W) 

— ^= (^/T 7 - J -2t-U + 2W\ 
2V6 V / 

2 \ 2 C$2-l 2 + 2 C's2-2 2 + C 8 2-3 2 



|*>83 = 



1 15 

I 3 ' 2' T' r5,1/ 



1 



2 V3 



(\0duu) + \0udu) + \0uud) + \dOuu) + \duOu) + \duuO) 



+ \uOdu) + \uOud) + \udOu) + \uduQ) + \uuOd) + \uudO)) 



l*>84 = |3,^,|,r 5 , 2 ) 

= C*84-i (|020u) - |02u0) - |0u02) + |0w20) - |200u) + |20u0) + |u002) - |u020)) 

+ C 8 4-2 (\0duu) - \dOuu) - \uuOd) + \uudO)) 

+ C$ 4 -3 (\0udu) + \0uud) - \duOu) + \duuO) - \uOdu) - \uOud) - \udOu) + \uduO)) 
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Schumann: Extended Hubbard model on small clusters 



C84-2 
^84 



— U ( + J + 2t + U -2w) 
— 4= f \/%+J + 2t + U-2w) 

2 \J 2 C84-1 2 + C*84-2 2 + 2 C84-3 2 



l*>«5 = |3, |,|,r 5)2 ) 

= C 85 _i (|020u) - |02?i0) - |0u02> + |0w20) - |200u) + |20m0) + |u002) - |u020)) 

+ C85-2 (\0duu) - \dOuu) - \uuOd) + \uudO)) 

+ C 8 5-3 (\0udu) + \0uud) - \duOu) + \duuO) - \uOdu) - \uOud) - \udOu) + \uduO}) 



C85-1 
C85-2 

<?85-3 

^85 



J -2t-U + 2W 



) 



2V6 

— ^= [yfM- J-2t-U + 2w) 
4\/6 V / 

2 \ 2 Cs5-1 2 + Cs5-2 2 + 2 Cs5-3 2 



|*>86 = 



I*) 



87 



I 3 ' 2> ~' F 5,2) 
1 



2 V3 



flOdmt) + \0udu) + \0uud) - \dOuu) - \duOu) + \duuO) 



-\uOdu) - \uQud) - \udOu) + \uduO) - \uuOd) + \uudO)) 



1 3 

' 3 ' 2' 4' ' 



= C 87 _i (|002u> - |00u2) - |020u) + |0m02> + |20u0) - |2u00) - |u020) + |u200)) 
+ C , 87 _2 (\0duu) + \0udu) + \dOuu) - \duOu) + \uOud) - \uduO) - \uuOd) - \uudO}) 
+ C 8 7-3 (\0uud) - \duuO) + \uOdu) - \udOu}) 



<?87-3 

N 87 



-JM 

2 

— ^ f + J + 2t + U -2W 
4\/6 V 

— U ( + J + 2t + U -2W) 
2V6 V / 

2 \J 2 Cs7-i 2 + 2 C$t-2 2 + C\ 



) 



^87-3 
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l*>*» - |3, |,|,r 5>3 > 

= C 88 _i (|002u) - |00u2) - |020?i) + |0u02) + |20u0) - |2u00) - |u020) + |u200» 

+ Cs8-2 (\0duu) + \0udu) + \dOuu) — \duOu) + \uOud) — \uduO) — \uuOd) — \uudO)) 

+ C 8 8-3 (\0uud) - \duuO) + \uOdu) - \udOu)) 

Css-i - g 

Css-2 = -A= (\/A^-J-2t-U + 2W 



4V6 

C 88 -3 = ^(/M-J-2t-U + 2W) 



N$g — 2 \J 2 C 8 8-i 2 + 2 C 8 8-2 2 + C; 



88-3 



l*>89 = |3,^,r 5 , 3 ) 



— !-= (\0duu) + \0udu) + \0uud) + \dOuu) - \duOu) - \duuO) 
2 V3 

+ \uOdu) + \uOud) — \udOu) — \uduO) — \uuOd) — \uudO}) 



Unnormalized eigenvectors for N e = 3 and m s = |. 



!*)«, = |3,|^,r 2 ) 



— (\0uuu) — \uOuu) + \uuOu) — \uuuO)) 



Q 1 K 

|*) 91 = |3,-, T ,r 5jl ) 

= i (|0miw) + \uOuu) + \uuOu) + \uuuO)) 
Q 1 e; 

|*)92 = |3,-, T ,r 5 , 2 > 

= — (|0uuii) — \uOuu) — \uuOu) + \uuuO)) 

l*>93 = |3,|f,r 5 , 3 ) 

= i (|0uuu) + |uOuu) — \uuOu) — \uuuO)) 
Unnormalized eigenvectors for N e = 4 and m s = —2. 

|*} 94 - |4,-2,6,r 2 ) 
= \{\dddd)) 
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Schumann: Extended Hubbard model on small clusters 



Unnormalizcd eigenvectors for N e = 4 and m s = — 1. 



|*} 95 - |4,-i,2,r 2 ) 

= — !-= (\02dd) - \Qd2d) + \Qdd2) - \20dd) + \2d0d) - \2ddQ) 
2 v3 

+ \d02d) - \d0d2) - \d20d) + \d2d0) + \dd02) - \dd20)) 



|*} 96 = |4,-i,6,r 2 ) 

= - (\dddu) + \ddud) + \dudd) + \uddd}) 



|*} 97 = |4,-i,2,r 3 ,i) 

= -^-j= {\Q2dd) - \Qdd2) - \20dd) + \2dd0) - \d02d) + \d20d) + \dd02) - \dd2Q)) 

|*} 98 = |4,-i,2,r 3 ,2) 

= C* 98 _i (|02dd) + \Qdd2) - \20dd) - \2dd0) + \d02d) - \d20d) + \dd02) - \dd20}) 

+ C 98 -2 (|0d2d) - \2d0d) + \d0d2) - \d2d0)) 

C98-1 = 



a 



2V6 

1 

N g8 = 2 V /2C 98 _ 1 2 + C< 



98-2 



|*)99 = |4,-i,2,r 4 ,i) 

= -^-j= {\M2d) + \0dd2) + \2d0d) + \2dd0) + \d02d) + \d0d2) + \d20d) + \d2d0)) 

|*)ioo - |4,-l,2,r 4> 2> 

= — !-= {\Q2dd) + \0d2d) + \20dd) + \2d0d) - \d0d2) - \d2d0) - \dd02) - \dd20}) 
2 v 2 

|*)ioi = |4,-l,2,r 4;3 ) 

= (\02dd) - \0dd2) + \20dd) - \2dd0) + \d02d) + \d20d) + \dd02) + \dd20)) 

2 v 2 

|*} 10 2 - |4,-i,2,r B ,i) 

= C102-1 (\02dd) + \0dd2) + \20dd) + \2dd0) - \d02d) - \d20d) + \dd02) + \dd20)) 
+ C102-2 (\0d2d) - \2d0d) - \d0d2) + \d2d0)) 
+ C102-3 (\dddu) — \ddud) + \dudd) — \uddd)) 
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C102-1 = 
C102-2 = 
C102-3 = 

N W2 = 

I*) 103 = 

+ 
+ 

Cl03-1 = 
C103-2 = 
Cl03-3 = 

^103 = 

I*) 104 = 

+ 
+ 

Cl04-1 = 
C104-2 = 
Cl04-3 = 

N W4 = 

I*) 105 = 

+ 
+ 



-9 J 2 + 288 t 2 + 12 An U - 36 J U 

-36 U 2 + 120 A n W - 360 J W - 720 U W - 3600 W 2 ) 

2 \J 2 C102-1 2 + C102-2 2 + C102-3 2 

|4,-i,2,r B ,i) 

C103-1 (|02dd) + |0rfd2) + |20dd) + |2dd0) - \d02d) - \d20d) + \dd02) + \dd20}) 
C103-2 (\0d2d) - \2d0d) - \d0d2) + |d2d0}) 
C103-3 (\dddu) — \ddud) + \dudd) — \uddd)) 

t -{2A w -'iJ-QU-QQW) 



-9 J 2 + 288 1 2 + 24 A\8 U — 36JU 

-36U 2 + 240 A 1S W - 360 J W - 720 U W - 3600 W 2 ) 

2 \J 2 C103-1 2 + C103-2 2 + C103-3 2 

|4,-i,2,r B ,i) 

C104-1 (|02dd) + |0<icZ2> + |20rfd) + |2drf0) - |d02d) - \d2Qd) + \dd02) + \dd20)) 
C104-2 (|0d2d) - |2rf0d) - \d0d2) + |d2d0» 
C104-3 (\dddu) — \ddud) + \dudd) — \uddd)) 

I (2 Aw - 3 J - 6 U - 60 W) 
-4t 2 

(-4 A i6 2 + 12 A w J 

—9 J 2 + 288 i 2 + 24 A w U — 36 JU 

-36 U 2 + 240 Ai 6 T4 7 — 360 JW — 720 UW - 3600 V^ 2 ) 



|4,-i,2,r B , 2 > 

C105-1 (|02dd) - \20dd) - \dd02) + |dd20)) 

C105-2 (|0d2d) - |0dd2) + |2d0d) - |2dd0) - \d02d) + \d0d2) - \d20d) + \d2d0}) 
C105-3 (\dddu) + \ddud) — \dudd) — \uddd)) 
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Schumann: Extended Hubbard model on small clusters 



C105-1 
Cl05-2 

C105-3 
^105 



4t 2 
t 



(j + U -2W + 2 ^/M cos(9 3 )) 



72 



{-Mi 2 + 6 An J 

—9 J 2 + 288 i 2 + 12 An U — 36JU 

-36 C/ 2 + 120 An - 360 J W - 720 U W - 3600 W 2 ) 

\J C105-1 2 + 2 C105-2 2 + C105-3 2 



|*)i06 = |4,-l,2,r 5 , 2 > 

= Cioe-1 (|02dd) - \20dd) - \dd02) + \dd20}) 

+ Cioe-2 (\0d2d) - \0dd2) + \2d0d) - \2dd0) - \d02d) + \d0d2) - \d20d) + \d2d0}) 

+ C106-3 (\dddu) + \ddud) — \dudd) — \uddd)) 



C106-1 
C106-2 



C 



106-3 — 



^106 



4r 
-t 

1 

72 



(2A 18 -3J-6f7-60W0 

(-4 A i8 2 + 12 Ai 8 J 

-9 J 2 + 288 t 2 + 24 Ai 8 U — 36JU 

-36 [7 2 + 240 Ai 8 T4 7 — 360 JW — 720 UW — 3600 V^ 2 ) 



2 C106-1 2 + 2 C106-2 2 + £106-3" 



|*)l07 = |4,-l,2,r B ,2) 

= C107-1 (|02dd) - |20<ta) - |rfd02> + |dd20» 

+ C107-2 (\0d2d) - \0dd2) + \2d0d) - \2dd0) - \d02d) + \d0d2) - \d20d) + \d2d0}) 

+ C107-3 (\dddu) + \ddud) — \dudd) — \uddd)) 



C107-1 
C107-2 

C107-3 



N 



107 



= At 1 

-t 



— (2Ai 6 -3 J-6Z7-60W) 



-9 J 2 



12 Ai6 J 

288t 2 + 24Ai 6 ?7-36 J{7 



36 /7 2 + 240Ai 6 VF-360 JW- 720J7V^-3600 VK 2 



= 2 



\/ C107-1 2 + 2 C107-2 2 + C107-3 2 



|*)ios = |4,-l,2,r 5 , 3 > 

= Cios-1 (\02dd) - \0d2d) + \20dd) - \2d0d) + \d0d2) + \d2d0) - \dd02) - \dd20)) 

+ C108-2 (|0dd2) - \2dd0) - \d02d) + \d20d)) 

+ C108-3 (\dddu) — \ddud) — \dudd) + \uddd)) 



adp header will be provided by the publisher 



67 



C108-2 = 4i 2 

Cios-3 = (An 2 — 6 An J 

+ 9 J 2 - 288 1 2 - 12 A U U + 36 JU 

+ 36U 2 - 120 An W + 360 JW + 720UW + 3600 VP 2 ) 



^Vio8 — 2 2 Cios-i 2 + C'108-2 2 + C; 



108-3 

|*) 109 = |4,-i,2,r 5 , 3 > 

= C109-1 (\02dd) - \0d2d) + \20dd) - \2d0d) + \d0d2) + \d2d0) - \dd02) - \dd20}) 
+ C109-2 (|0dd2) - \2dd0) - \d02d) + \d20d)) 
+ C109-3 (\dddu) — \ddud) — \dudd) + \uddd)) 

C109-1 - I (2Ai 8 -3 J-6U-60W) 


2 



C109-2 — 4t 
C109-3 = ^ (4A 18 2 -12A 18 J 



9 J 2 — 288 1 2 — 24 Ai 8 U + 36JU 

36 t/ 2 - 240 Ai S W + 360 JW + 720 UW + 3600 TP 2 ) 



-Wl09 — 2y2Ci09-l + Cio 9 -2 + Cio 9 -3 

|*)no = |4,-l,2,r B , 3 > 

= C110-1 (\02dd) - \0d2d) + \20dd) - \2d0d) + \d0d2) + \d2d0) - \dd02) - \dd20)) 

+ Cno-2 (|0dd2) - \2dd0) - \d02d) + \d20d}) 

+ C110-3 (\dddu) — \ddud) — \dudd) + \uddd)) 

C110-1 = I {2A m -3J-6U-60W) 


Cno-2 = 4i 2 

C110-3 = ^ (4A 16 2 -12A 16 J 

+ 9 J 2 - 288 1 2 - 24 A 16 U + 36 J U 

+ 36U 2 - 240 Ai 6 TP + 360 JW + 720 UW + 3600 TP 2 ) 



^110 — 2 \J 2 Cno-i 2 + Cno-2 2 + C110-3 2 



Unnormalized eigenvectors for N e = 4 and m s = 0. 

|*)m = |4,o,o,ri) 

= C111-1 (|0022) + |0202) + |0220) + |2002) + |2020) + |2200» 
+ C m _2 (|02du) - |02ud> + |0d2u) + \0du2) - \Qu2d) - \Qud2) 
+ |20dw) - \2Qud) + \2dQu) + \2du0) - \2u0d) - \2ud0) 
+ |rf02u) + |rf0u2) + |d20u) + \d2uQ) + \duQ2) + \du2Q) 
-\vJ02d) - \u0d2) - \u20d) - \u2d0) - \ud02) - \ud20}) 
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Schumann: Extended Hubbard model on small clusters 



Clll-1 
Cm-2 
Will 



2 V^ 



T*^ +J+U ~ 2W ) 



V /6C 111 _ 1 2 + 24C 111 _2 2 



|*)n2 - |4,o,o,ri) 

= C112-1 (|0022) + |0202) + |0220) + |2002) + |2020) + |2200>) 
+ C112-2 (|02du) - |02ud> + |0d2u) + \0du2) - \0u2d) - \0ud2) 
+ \20du) - \20ud) + \2d0u) + \2du0) - \2u0d) - \2ud0) 
+ |rf02u) + |rf0u2) + |rf20u) + \d2u0) + \du02) + \du20) 
-\u02d) - \u0d2) - \u20d) - \u2d0) - \ud02) - \ud20}) 



C112-1 
C112-2 
N 112 



4V6 



M 3 - J-U + 2W 



6C112-1 +24Ci 12-2 



I*) 



113 



|4,0,2,r 2 ) 



2 V6 



(\02du) + \Q2ud) - \0d2u) + \0du2) - \0u2d) + \0ud2) 

-\20du) - \2Qud) + \2d0u) - \2du0) + \2u0d) - \2ud0) 
+ |rf02u) - |rf0u2) - \d20u) + \d2u0) + \du02) - \du20) 
+ \u02d) - \u0d2) - \u20d) + \u2d0) + \ud02) - \ud20)) 



|*>ii4 = |4,0,6,r 2 ) 

= —= (\dduu) + \dudu) + \duud) + \uddu) + \udud) + \uudd)) 
V6 



|*)hb - |4,o,o,r 3 ,i) 

= C115-1 (|0022) + |0220) + |2002) + |2200» 

+ C115-2 (|0202) + 12020)) 

+ C115-3 (|02du) - \Q2ud) + \0du2) - \0ud2) + \20du) - \20ud) + \2du0) - \2ud0) 

+ \d02u) + \d20u) + \duQ2) + \du20) - \u02d) - \u20d) - \ud02) - \ud20)) 

+ C115-4 (\Qd2u) - \0u2d) + \2d0u) - \2u0d) + \d0u2) + \d2u0) - \u0d2) - \u2d0)) 

+ Cn5_5 (\dduu) — \duud) — \uddu) + \uudd)) 
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Cll5-1 — 

C115-2 = ^ (V3J + V3C/-2V3T^-2v^cos(ei)) 

C115-3 = — V (3 J 2 + 6 JU + 3U 2 -12 JW -12UW 
12 v3 

+ 12W 2 -4A lC os(ei) 2 ) 

C115-4 = -A= (3 J 2 + 6 JU + 3U 2 -12 JW -12UW 
6 V3 



+ 12 VP* -4Aicos(9ir 
C115-5 = * (V3 J + V3 [/ - 2 \/3 T4^ + 2 V^T cos(9i )) 

+ 2 Cn5_2 2 + 16 Cn5_3 2 + 8 Cn5_4 2 + 4 Cn5_5 2 

|*)na - |4,0,0,r 3 ,i) 

= Cii 6 _i (|0022) + |0220) + |2002) + |2200>) 
+ Cii6_ 2 (|0202) + 12020)) 

+ Ciie-3 (|02du) - |02ud> + |0du2) - \0ud2) + \20du) - \20ud) + \2du0) - \2ud0) 

+ \d02u) + \d20u) + \du02) + \du20) - \uQ2d) - \u20d) - \udQ2) - \ud20)) 
+ Ciie-4 (|0d2u) - |0w2d) + |2d0u) - |2u0d) + |rf0u2) + |d2u0) - |u0d2) - |u2d0» 
+ C116-5 (\dduu) — \duud) — \uddu) + \uudd)) 

C116-1 = ^ (y3J + V3U-2V3W+ V^7cos(6i) + \/3v / ^7siri(ei)) 

C116-2 = ^ J + V3U- 2 V3 1^ + ^003(61) + V3 V^isin(0i)) 

Cii 6 _3 = — f2V3A!-3\/3J 2 -6V3JC/-3\/3[/ 2 + 12\/3JT^ 
36 V 

+ 12 V3UW- 12V3VK 2 - V3Ai cos(2 6i) + 3 A 1 sin(2 6i)) 

Cii 6 _4 = (-2V3Ai +3V3 J 2 + 6V3 JU + 3V3U 2 - 12 V3 JW - 12V3UW 
18 V 

+ 12V3VK 2 + V3Ai cos(29i)-3Ai sin(26i)) 
C116-5 = t (V3 J + V3U- 2V3W - i/M cos(6i) - V3 V^i sin(6i)) 



iViie — y 4Cn6-i +2C116-2 +I6C116-3 +8C116-4 +4Cn6-5 

|*) 117 - |4,o,o,r 3 ,i) 

= Cii 7 _i (|0022) + |0220) + |2002) + |2200» 

+ C117-2 (10202) + 12020)) 

+ C117-3 (|02du) - |02wd) + |0du2) - |0ud2) + |20dw) - |20?id) + |2du0) - |2ud0) 

+ |d02u) + |rf20u) + |du02) + |du20) - |u02rf) - |u20rf) - |ud02) - |ud20)) 

+ C117-4 (|0d2u) - |0w2d) + |2d0u) - |2u0rf) + |rf0u2) + |d2u0) - |u0d2) - |u2dO)) 

+ Cn7_5 (\dduu) — \duud) — \uddu) + \uudd)) 
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Schumann: Extended Hubbard model on small clusters 



C117-1 = ^ (y3J + V3U-2V3W+y r A^cos{e 1 )-V3y / %sm(e 1 )^ 

C117-2 = ^ (V3J + \/3U-2\/3W + y/A^cos(e 1 )-\/3y/Msm{e 1 )j 

C117-3 = — (2V3A 1 -3V3J 2 -6V3JU -3V3U 2 + 12V3JW 
36 V 

+ 12V3UW - 12 V3W 2 - V3^i cos(2 6i) -3 Ax sin(2 6i)) 

C117-4 = ^ (-2V3Ai +3V3 J 2 + 6V3 JU + 3V3U 2 - 12V3 JW- 12V3UW 
18 V 

+ 12V3W 2 + V3Ai cos(29i) + 3Ai sin(26i)) 
C117-5 = i (V3 J + V3U- 2V3W - V^i cos(6i) + V3 V^i sin(6i)) 
Nin = \J 4 Cn7_i 2 + 2 Cn7_ 2 2 + 16 Cn7_3 2 + 8 Cn7_4 2 + 4 Cn7_ 5 2 

|*)ii 8 = |4,0,2,r 3 ,i> 

= i (|02du) + \02ud) - \0du2) - \0ud2) - \20du) - \20ud) + \2du0) + \2ud0) 

-\d02u) + \d20u) + \du02) - \du20) - \uQ2d) + \u20d) + \ud02) - \ud20)) 

|*)n9 = |4,o,o,r 3 , 2 ) 

= C119-1 (|0022) - |0220> - |2002) + |2200>) 

+ C lm - 2 (\02du) - \02ud) - \0du2) + \0ud2) + |20dw) - |20ud) - |2du0) + |2ud0) 

-\d02u) - |rf20u) + \du02) + |rfu20) + |u02rf) + \u20d) - |ud02) - \ud2Q)) 
+ Cn9_3 (\dduu) + \duud) + \uddu) + \uudd)) 
+ Cii9_4 (\dudu) + \udud)) 

C119-1 = 

— (-3J-3C/ + 6VK + 2V3 yfM cos(0i)) 
Cii9_ 2 = ^ (-3J 2 -6 JU-3U 2 + 12 JW + 12UW 

-12W 2 + AA 1 cos(9i) 2 ) 
Ciig-3 = ^ (-3 J-3U + 6W-2^/3^/A^cos(@ 1 j^ 
di9_4 = —7^- (-3 J - 3U + 6W — 2 V3 \[~A\ cos(9i)^ 

iVng = \j 4 Cn9_i 2 + 16 Ciig_2 2 + 4 Cng_ 3 2 + 2 Cng_4 2 

|*>i2o = |4,o,o,r 3 , 2 ) 

= C120-1 (|0022) - |0220) - |2002) + |2200>) 

+ C120-2 (\02du) - \02ud) - \0du2) + \0ud2) + \20du) - \20ud) - \2du0) + \2ud0) 

-\d02u) - \d20u) + \du02) + \du20) + \u02d) + \u20d) - \ud02) - \ud20}) 
+ C120-3 (\dduu) + \duud) + \uddu) + \uudd)) 
+ C120-4 (\dudu) + \udud)) 
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C120-1 — 

C120-2 = 

C120-3 = 

C120-4 = 



-t 

T 
1 

12 



(-3J-3C/ + 6VK-V3 cos(6i) - 3 sin(6i)) 
(2 Ai - 3 J 2 -6JU — 3 U 2 + 12 JW + 12 UW 
-12 W 2 - At cos(2 6i) + V3Ai sin(2 6i)) 
I (-3 7- 3^ + 61^ + ^3 cos(6i) + 3 sin(6i)) 
^ (-3J-3?7 + 6VK + \/3 V^" cos(0i) + 3 ^f~K[ sin(©i)) 



^120 = 


^/ 4 C120-1 2 + 16 C120-2 2 + 4 C120-3 2 + 2 C120-4 2 








|*>121 = 


|4,o,o,r 3 , 2 > 










C121-1 (|0022) - |0220) - |2002) + |2200>) 








+ 


C121-2 (\02du) - \02ud) - \0du2) + \0ud2) + \20du) - 


\2Qud) - 


|2du0) ^ 


- |2ud0) 




-\d02u) - \d20u) + \du02) + \du20) + \u02d) - 


h \u20d) - 


- |ud02) 


- |itd20)) 


+ 


C121-3 (\dduu) + \duud) + \uddu) + \uudd)) 








+ 


Ci2i_4 (\dudu) + \udud)) 








C12I-I = 


Y (-3J-3U + QW-V3 \fA~x cos(6i) + 3 


sin(0i)) 







C121-2 — 

C121-3 = 
C121-4 = 

#121 
|*>122 



12 



(2Ai - 3 J 2 - 6 J U — 3U 2 + 12 JW + 12UW 
-12 W 2 - A x cos(2 6i) - V3 Ai sin(2 6i)) 
(^-3J-3?7 + 6T^ + \/3 cos(0i) - 3 v^i sin(0i 
(-3J-3J7 + 6VK + V3 v/IT cos(0i) - 3 yfM sin(9i)) 



-2i 



= \/ 4 C'lai-i 2 + 16C 121 _ 2 2 +4C 12 i- 3 2 + 2 C_ 



+ 

C122-1 = 
C122-2 = 

#122 = 
I*) 123 = 

+ 
+ 



121-4 

|4,o,2,r 3 , 2 > 

C122-1 (\02du) + \02ud) + \0du2) + \0ud2) - \20du) - \20ud) - \2du0) - \2ud0) 

+ \d02u) - \d2Qu) + \duQ2) - \du20) + |u02rf) - \u20d) + |ud02) - |ud20)) 
C122-2 (|0d2u) + |0u2d) - |2d0u) - \2u0d) + |d0u2) - |d2u0) + |u0d2) - |u2d0)) 

-1 



Ay/3 

-1 

2 V3 



16 Ci22-1 2 + 8 Ci22-2 2 



|4,o,o,r 4 ,i) 

C123-1 (|0022) - |2200)) 

C123-2 (|02du) - |02wd) + |20du) - |20ud) - |(M)2) - |du20) + |ud02) + |ud20)) 
C123-3 (|0d2u) + |0du2) - |0u2d) - |0ud2) - |2rf0u) - |2du0) + |2u0d) + |2ud0) 

+ \d02u) + |rf0u2) - |d20u) - |d2u0) - |u02rf) - |u0rf2) + |u20d) + |u2d0» 
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Schumann: Extended Hubbard model on small clusters 



Cl23-1 
C123-2 

C123-3 
^123 



= -4\/2t 2 



9V2 



(2 A10 2 - 3 A w J - 36 t 2 + 9 A w U - 9 J U 



+ 9 1/ 2 + 54 A10 - 36 J W + 126 U W + 360 W 2 ) 
(V2t) 



( 



J + C/-2VF 



'A 2 cos 



<e 2 )) 



= ^/ 2 C123-1 2 + 8 C123-2 2 + 16 C123-3 2 



I*)i24 - |4,o,o,r 4 ,i) 

= C124-1 (|0022) - |2200)) 

+ C124-2 (|02du) - \02ud) + \20du) - \20ud) - \duQ2) - \du20) + \ud02) + \ud20)) 
+ C124-3 (\Qd2u) + \0du2) - \0u2d) - \0ud2) - \2d0u) - \2duQ) + \2u0d) + \2ud0) 

+ \d02u) + \d0u2) - \d20u) - \d2uQ) - \uQ2d) - \u0d2) + \u20d) + \u2d0)) 



C124-1 
C124-2 

C124-3 



-4V2t 2 



18 V2 



(A 14 2 + 3A 14 J- 72 1 2 -9A 14 U- 18 JU 



+ 18U 2 -54 A 14 W -72 JW + 252 UW + 720 W 2 ) 
^ (2J+2U-4W- cos(9 2 ) - V3 sin(6 2 ) 



3 V2 

2 C124-1 2 + 8 C124-2 2 + 16 C 



124-3 



|*)i2 B = |4,o,o,r 4 ,i) 

= C125-1 (10022) - |2200)) 

+ C125-2 (|02du) - |02wd) + |20du) - |20ud) - \du02) - \du20) + \ud02) + |ud20}) 
+ C125-3 (\0d2u) + \0du2) - \0u2d) - \0ud2) - \2d0u) - \2du0) + \2u0d) + \2ud0) 

+ \d02u) + \d0u2) - \d20u) - \d2uQ) - \u02d) - \u0d2) + \u20d) + \u2d0)) 



C125-1 
C125-2 

C125-3 

iVi25 



= -4\/2t 2 



18 V2 



-t 
3 V2 



(A 13 2 + 3A 13 J- 72 1 2 -9A 13 U- 18 JU 
+ 18U 2 -54 AwW -72 JW + 252UW + 720W 2 ) 
(2J+2U -AW - i/M cos(9 2 ) + V3 sin(6 2 )) 



\J 2 C125-1 2 + 8 C125-2 2 + 16 C125-3 2 



|*)i26 = |4,0,2,r 4 ,i) 

= i (\Qd2u) + \0du2) + \Qu2d) + \Qud2) + \2d0u) + \2duQ) + \2u0d) + \2ud0) 

+ \d02u) + \d0u2) + \d20u) + \d2u0) + \uQ2d) + \uM2) + \u20d) + \u2d0)) 



adp header will be provided by the publisher 



73 



I*)i27 = |4,o,o,r 4 , 2 ) 

= C127-1 (|0220) - |2002)) 

+ C127-2 (|02du) - \02ud) + \0d2u) - \0u2d) - \20du) + \20ud) - \2d0u) + \2u0d) 

-\d0u2) + \d2uQ) - \duQ2) + \du20) + \u0d2) - \u2d0) + \ud02) - \ud20)) 
+ C127-3 (|0du2) - |0wd2) + |2du0) - \2ud0) - \d02u) - \d2Qu) + \u02d) + \u20d)) 

C127-1 - 0= (4A W 2 -12A W J 

+ 9 J 2 - 72 t 2 + 12 A w U - 18 J U 

+ 9 U 2 + 120 A w W - 180 J W + 180 U W + 900 W 2 ) 



C127-2 = -j= (j+[/-2VK-2v^cos(e 2 )) 



3 V2 

C127-3 = — 2\/2i 2 

7V127 = \J 2 C127-1 2 + 16 C127-2 2 + 8 C127-3 2 

|*)i28 = |4,o,o,r 4 , 2 ) 

= C128-1 (|0220) - |2002)) 

+ C128-2 (|02du) - \Q2ud) + \0d2u) - \0u2d) - \20du) + \20ud) - \2d0u) + |2u0d) 

-|d0u2> + \d2u0) - |du02) + |du20) + |u0d2) - |u2d0) + |ud02) - |ud20)) 
+ C128-3 (|0du2) - |0wd2) + |2d?i0) - |2wd0) - |d02w) - |d20u) + |u02d) + |u20d» 

C128-1 = (A14 2 + 6 Au J 

+ 9 J 2 - 72t 2 -6AuU- 18 JU 

+ 9 U 2 - 60 A 14 W - 180 J W + 180 U W + 900 W 2 ) 

C128-2 = -A= (J+U-2W+ v / ^ 2 "cos(e 2 ) + V3 V^2"sin(e 2 )) 



3\/2 

C128-3 = —2\[2t 



2 



-^128 — Y 2 Ci 2 8-l 2 + 16 Ci28-2 2 + 8 Ci28-3 2 

|*)i29 = |4,o,o,r 4 , 2 ) 

= Ci 29 _i (|0220) - |2002)) 

+ C129-2 (|02du) - \Q2ud) + \0d2u) - \0u2d) - \20du) + \20ud) - \2d0u) + |2u0d) 

-\dOu2) + \d2uQ) - \duQ2) + |du20) + \u0d2) - |u2d0) + |ud02) - |ud20)) 
+ C129-3 (|0du2) - |0wd2) + |2du0) - |2wd0) - |d02w) - |d20u) + |u02d) + \u20d)) 

C129-1 = (A 13 2 + 6 A 13 J 

7 2 70 +2 



9 J z - 72 r -6A 13 U- 18 JU 
lis 



9 U 2 - 60 A 13 W - 180 J W + 180 U W + 900 W 2 ) 



C129-2 = jj= (j+t/-2VK+ cos(6 2 ) - V3 V^2"sin(e 2 )) 

-2 



3\/2 

C129-3 = — 2 a/2 t 2 
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Schumann: Extended Hubbard model on small clusters 



|*) 130 = |4,0,2,r 4 , 2 ) 

= j (|02du) + \02ud) + \Qd2u) + \0u2d) + \20du) + \2Qud) + \2d0u) + \2u0d) 

-|d0u2) - \d2u0) - \duQ2) - \du2Q) - \u0d2) - \u2d0) - \ud02) - \ud20}) 



|*)i3i = |4,o,o,r 4 , 3 > 

= C131-1 (|0202) - |2020)) 

+ C131-2 (\02du) - \Q2ud) + \0du2) - \0ud2) - \20du) + \20ud) - \2du0) + \2ud0) 

-\dQ2u) + \d2Qu) + \du02) - \du20) + \uQ2d) - \u20d) - \ud02) + |ud20)) 
+ C131-3 (|0d2u) - |0w2d) + |2d0u) - |2u0d) - \d0u2) - \d2uQ) + \u0d2) + \u2d0)) 



C 



131-1 — 



Ci3i_2 
Cl31-3 



-1 



9 V2 



(4A 10 2 



12 A w J 



+ 9 J 2 - 72 t 2 + 12 A w U - 18 J U 

+ 9 U 2 + 120 A w W - 180 J W + 180 U W + 900 W 2 ) 



^— (j+U -2W-2 ^A~ 2 cos(e 2 )) 



-2V2t 2 



\J 2 Ci3i_i 2 + 16 Ci3i_2 2 + 8 Ci3i_3 2 



|*)i32 = |4,o,o,r 4 , 3 > 

= C132-1 (10202) - |2020)) 

+ C132-2 (\02du) - \02ud) + \0du2) - \0ud2) - \20du) + \20ud) - \2du0) + \2ud0) 

-\d02u) + \d20u) + \du02) - \du20) + \u02d) - \u20d) - \ud02) + \ud2Q)) 
+ Ci32-3 (\0d2u) - \0u2d) + \2d0u) - \2u0d) - \d0u2) - \d2u0) + \u0d2) + \u2d0)) 



Cl32-1 — 

C*132-2 = 

Cl32-3 = 

^132 = 



-1 

9 V2 



(A 14 2 + 6A U J 



9 J 2 - 72 1 2 - 6 A 14 , U - 18 JU 



+ 9 U 2 - 60 A u W - 180 J W + 180 U W + 900 W 2 ) 
jj= ( J + U - 2 W + ^T 2 cos(6 2 ) + V3 sin(9 2 )) 
-2 V2t 2 

\j 2 Ci32-1 2 + 16 Ci32-2 2 + 8 Ci32-3 2 



|*)i33 = |4,o,o,r 4 , 3 ) 

= C133-1 (10202) - |2020)) 

+ Ci33-2 (|02du) - |02wd) + |0du2) - |0ud2) - |20dw) + |20ud) - |2du0) + |2ud0) 

— |d02u) + |d20u) + |du02) - |du20) + |u02d) - |u20d) - |ud02) + |wd20)) 
+ C133-3 (\0d2u) - \0u2d) + \2d0u) - \2u0d) - \d0u2) - \d2uQ) + \u0d2) + |u2d0» 
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Ci33_i = 



Ci33_2 

C*133-3 
N 133 



-1 



3V2 
-2\f2t 2 



(A 13 2 + 6A 13 J 

+ 9 J 2 - 72t 2 -6A 13 U - 18 JU 

+ 9 U 2 - 60 A 13 W - 180 J W + 180 U W + 900 W 2 ) 

( J + U - 2 W + y]~M cos(6 2 ) - V3 sin(9 2 )) 



2C 



133-1 



16 C 



133-2 



8C 



133-3 



|*)i34 = |4,o,2,r 4 , 3 > 



- (|02du> + \i)2ud) - \Mu2) - \Qud2) + \20du) + |20ud> - \2duQ) - \2ud0) 
+ |d02u) + \d2Qu) + \du02) + \du20) + \uQ2d) + \u2Qd) + |ud02) + |ud20)) 



I*) 



135 



|4,o,o,r 5l i) 

i (|02du) - \Q2ud) - \0du2) + \0ud2) - \20du) + \20ud) + \2du0) - \2ud0) 
+ \d02u) - \d20u) + \du02) - \du20) - \uQ2d) + \u20d) - \ud02) + \ud20)) 



|*)i36 = |4,0,2,r 5 ,i) 

= C136-1 (|02du) + |02ud> + \0du2) + \0ud2) + \20du) + \20ud) + \2du0) + \2ud0) 

-\dWZu) - |rf20u) + \du02) + \du20) - \u02d) - \u20d) + \ud02) + \ud2Q)) 
+ C136-2 (|0d2u) + \0u2d) - \2d0u) - \2u0d) - \d0u2) + \d2uQ) - \u0d2) + \u2dQ)) 
+ Ci36-3 (\dudu) — \udud)) 



Cl36-1 

Cl36-2 
Cl36-3 



N 



136 



3 V2 
-2 V2t 2 
-1 



J + U -2W + 2 JM cos 



(©3)) 



9\/2 



(A 12 2 - 72 1 2 + 6 J U - 3 U 2 + 60 J W - 132 U W - 1020 W 2 
+ 12^/MU cos(9 3 ) + 120 ^/MW cos(9 3 )) 



I6C136-I +8C136-2 +2Ci36-3 



|*) 137 - |4,0,2,r 5 ,i) 

= C137-1 (\02du) + \02ud) + \0du2) + \Qud2) + |20dw) + \2Qud) + \2du0) + \2ud0) 

-\dQ2u) - \d2Qu) + \du02) + \du20) - \u02d) - \u20d) + \ud02) + \ud20}) 
+ C137-2 (|0d2u) + \0u2d) - \2d0u) - \2u0d) - \d0u2) + \d2u0) - \u0d2) + |u2dO» 
+ Ci37_3 (\dudu) — \udud)) 
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Schumann: Extended Hubbard model on small clusters 



CW-i = jj= (J + U-2W- y[M cos(e 3 ) - V3 y[M sin(6 3 )) 

C137-2 = — 2v / 2f 2 

C137-3 = ^ (A 17 2 - 72t 2 -6A 17 (7 



+ 8 Ci37_2 + 2 C137- 



2 

3 



|*) 138 = |4,0,2,r B ,i) 

= C138-1 (|02du) + |02ud) + |0du2) + |0ud2) + |20dw) + \20ud) + \2du0) + \2ud0) 

-\dl02u) - \d20u) + \du02) + \du20) - \u02d) - \u20d) + \ud02) + \ud20}) 
+ C138-2 (\Qd2u) + \0u2d) - \2d0u) - \2u0d) - \d0u2) + \d2u0) - \u0d2) + \u2d0}) 
+ C138-3 (\dudu) — \udud)) 

Cias-i = —J= (J+U-2W- v^b cos(6 3 ) + V3 sin(6 3 )) 



3 V2 

C138-2 = -2V2t 2 

C 13S - 3 = ^= (A 15 2 - 72 t 2 - 6 A 15 (7 



+ 9 U 2 - 60 A15 + 180 U W + 900 W 2 ) 

+ 8Ci 3 8-2 +2Ci 3 8-3 



|*)i 39 = |4,o,o,r B , 2 ) 

= i (|0rf2u) - |0dw2) - |0u2d) + |0ud2) - |2d0u) + \2du0) + \2u0d) - \2ud0) 

-\d02u) + \d0u2) + \d2Qu) - \d2uQ) + |u02rf) - |u0d2) - |u20d) + |u2rf0)) 

|*>i4o - |4,o,2,r 5 , 2 ) 

= C140-1 (|02du) + |02ud> - \20du) - \2Qud) - \duQ2) + \du20) - \ud02) + \ud20)) 
+ C140-2 (|0rf2u) - |0rfu2) + \0u2d) - \0ud2) + \2d0u) - \2du0) + \2u0d) - \2ud0) 

-\d02u) + \d0u2) - \d20u) + |rf2u0) - |u02rf) + \u0d2) - \u20d) + \u2d0)) 
+ C140-3 (\dduu) — \uudd)) 

C140-1 = 2^/2t 2 
t 



C 140 -2 = 77-^= ( J + U - 2 W + 2 yfM cos(9 3 )) 

C140-3 = — i (A 12 2 - 72< 2 + 6JC/-3[/ 2 + 60JVK-132[/T^- 1020T^ 2 
9 v 2 

+ 12^/MU cos(9 3 ) + 120 ^/%W cos(6 3 )) 
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|*)ui - |4,0,2,r 5 , 2 ) 

= C 14 ,i-i (\02du) + \02ud) - \20du) - \20ud) - \du02) + \du20) - \ud02) + \ud20)) 
+ C141-2 (|0d2u) - \0du2) + \0u2d) - \0ud2) + \2d0u) - \2du0) + \2u0d) - \2ud0) 

-\d02u) + \d0u2) - \d20u) + \d2uQ) - \uQ2d) + \u0d2) - \u20d) + \u2d0)) 
+ Ci4i_3 (\dduu) — \uudd)) 



C141-1 

^141-2 



2y/2t 2 



^— (J + U-2W- y/M cos(6 3 ) - V3 t/M sin(9 3 )) 



C 



141-3 



2,^2 

-1 

971 



(A 17 2 - 72 1 2 - 6A 17 U 



N 



141 — 



+ 9 U 2 - 60 A 17 W + 180 U W + 900 W 2 ) 

+ 16Ci4i_2 +2Ci4i_3 



|*) U2 = |4,o,2,r 5 , 2 > 

= C142-1 (|02du) + \Q2ud) - \20du) - \2Qud) - \duQ2) + \du20) - \ud02) + \ud20)) 
+ C142-2 (|0d2u) - \0du2) + \Qu2d) - \0ud2) + |2rf0u) - |2du0) + |2u0d) - |2ud0) 

-\d02u) + \d0u2) - \d20u) + \d2uQ) - \uQ2d) + \u0d2) - \u20d) + |u2d0» 
+ C142-3 (\dduu) — \uudd)) 



C142-1 

Cl42-2 



2V2t 2 
t 



jj= ( J + U - 2 W - v^b cos(e 3 ) + V3 v^b sin(6 3 )) 



C 



142-3 — 



-1 

9 V2 



(A 15 2 - 72 1 2 -6 A 15 U 

+ 9 U 2 - 60 A 15 W + 180 U W + 900 W 2 ) 



N142 — \J 8 Ci42-1 2 + 16 Ci42-2 2 + 2 Ci42- 



|*)l43 = |4,0,0,I\3> 

= i (|02du) - \Q2ud) - \Qd2u) + \Qu2d) - \20du) + |20ud) + \2d0u) - \2u0d) 

+ \d0u2) - \d2u0) - \duQ2) + \du20) - \u0d2) + \u2d0) + \ud02) - \ud2Q)) 

|*) 144 = |4,0,2,r B , 3 ) 

= C144-1 (\02du) + \02ud) - \0d2u) - \0u2d) + \20du) + \20ud) - \2d0u) - \2u0d) 

+ \d0u2) + \d2u0) - \du02) - \du20) + \u0d2) + \u2d0) - \ud02) - \ud20)) 
+ C144-2 (\0du2) + \0ud2) - \2du0) - \2ud0) - \d02u) + \d20u) - \uQ2d) + |u20d» 
+ C144-3 (\duud) — \uddu)) 
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Schumann: Extended Hubbard model on small clusters 



Cl44-1 
C144-2 
C144-3 

iV 144 

|*)l46 



N 



145 



146 



C 



146-1 



-t 
371 

2 V2t 2 
-1 
9Vl 



(.J +U -2W + 2 cos(e 3 )) 



(Ai 2 2 - 72 t 2 + 6 J U - 3 1/ 2 + 60 J W - 132 U W - 1020 W 2 
+ 12^/MU cos(9 3 ) + 120 ^/MW cos(6 3 )) 



— \J 16 Ci44_i 2 + 8 Ci44_2 2 + 2 Ci44_3 2 

= |4,0,2,r B , 3 > 

= C145-1 (|02du) + |02ud) - |0d2u) - \Qu2d) + \20du) + \20ud) - \2d0u) - \2u0d) 

+ \d0u2) + |rf2u0) - \duQ2) - \du20) + \u0d2) + \u2d0) - \ud02) - \ud20)) 
+ C145-2 (|0du2) + \0ud2) - \2du0) - \2ud0) - \d02u) + \d20u) - \u02d) + \u20d)) 
+ Ci45_3 (\duud) — \uddu)) 



C145-1 — 

C145-2 = 
C145-3 = 



-t 



3 V2 
2V2t 2 
-1 
9V2 



(j + U-2W- ^J~A 2 cos(e 3 ) - V3 ^/M sin(6 3 )) 
2 

(A 17 2 - 72 1 2 -6A 17 U 
+ 9 U 2 - 60 A 17 W + 180 U W + 900 W 2 ) 



16 Ci 4 5_i 2 + 8 Ci45_2 2 + 2 Ci45_3 2 



= |4,0,2,r B , 3 > 

= Ci46-i (|02du) + |02ud) - |0d2u) - \0u2d) + \20du) + \20ud) - \2d0u) - \2u0d) 

+ \d0u2) + \d2u0) - \du02) - \du20) + \u0d2) + \u2d0) - \ud02) - \ud20)) 
+ Ci46-2 (|0du2) + |0wd2) - |2du0) - |2wd0) - |d02w) + |d20u) - |u02d) + |u20d)) 
+ C146-3 (\duud) — \uddu)) 



3 V2 

Cl46-2 = 2 V2 1 2 



l T (J + U-2W- ^fM cos(e 3 ) + V3 ^% sin(6 3 )) 



C 



146-3 



9 V2 



(A 15 2 - 72 1 2 - 6A 15 U 
+ 9 U 2 - 60 A 15 W + 180 U W + 900 W 2 ) 

Ni4Q = \J 16 Ci46-l 2 + 8 Ci46-2 2 + 2 Ci46- 3 2 

Unnormalized eigenvectors for N e = 4 and m s = 1. 



|*)u7 = |4,i,2,r 2 ) 

= — !-= (\Q2uu) - \Qu2u) + \0uu2) - \20uu) + \2u0u) - \2uu0) 
2 V3 

+ \uQ2u) - \u0u2) - \u2Qu) + \u2uQ) + \uu02) - \uu20)) 
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|f} 148 = |4,i,6,r 2 ) 



= - (\duuu) + \uduu) + \uudu) + \uuud)) 

|*) 148 = |4,i,2,r 3 ,i) 

= — !-= (\02uu) - \Quu2) - \2Quu) + \2uu0) - \u02u) + \u20u) + \uu02) - \uu20)) 
2 v 2 

|*)ibo = |4,l,2,r 3 ,2) 

= C150-1 (\02uu) + \0uu2) - \20uu) - \2uu0) + \u02u) - \u20u) + \uu02) - \uu20}) 

+ C 150 -2 (\0u2u) - \2u0u) + \u0u2) - \u2u0)) 



Cl50-1 
Cl50-2 
^150 



-1 

2 V6 



V6 



2V2(7i6o-i + C 



150-2 



151 = 



|4,i,2,r 4 ,i) 
i 



2 V2 



(|0u2u) + \0uu2) + \2u0u) + \2uu0) + \u02u) + \u0u2) + \u20u) + \u2u0)) 



|*)iB2 = |4,l,2,r 4 , 2 ) 

= — !-= {\Q2uu) + \0u2u) + \20uu) + \2u0u) - \u0u2) - \u2u0) - \uuQ2) - \uu2Q)) 
2 v 2 



I*)i53 = |4,i,2,r 4 , 3 > 



2 V2 



(|02uu) - \Quu2) + \20uu) - \2uu0) + \u02u) + \u20u) + \uu02) + \uu20}) 



154 



= |4,i,2,r 5 ,i) 

= C154-1 (|02uu) + |0ww2) + \20uu) + \2uu0) - \u02u) - \u20u) + \uuQ2) + \uu2Q)) 

+ C154-2 (|0u2u) - |2w0w) - \u0u2) + \u2u0)) 

+ C154-3 (\duuu) — \uduu) + \uudu) — \uuud)) 



C154-1 
C154-2 
C154-3 



- -I 
3 V 



J + U -2W + 2 JA 2 cos 



(©3)) 



-At 2 



72 



(-An 2 +6 A n J 



-9 J 2 + 288 t 2 + 12 A n U - 36 JU 

-36 U 2 + 120 A n W - 360 J W - 720 U W - 3600 W 2 ) 
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Schumann: Extended Hubbard model on small clusters 



|*} 155 - |4,i,2,r 5 ,i) 

= C155-1 (\02uu) + \0uu2) + \20uu) + \2uu0) - \u02u) - \u20u) + \uu02) + \uu20)) 
+ C155-2 (\0u2u) - \2u0u) - \u0u2) + \u2u0)) 
+ Ci55-3 (\duuu) — \uduu) + \uudu) — \uuud)) 

C1B5-1 = (2A 18 -3J-6U~mW) 


C155-2 = — 4t 2 

C155-3 - (-AA 1S 2 + 12A 18 J 

-9 J 2 + 288 t 2 + 24 Ais U — 36JU 

-36 [7 2 + 240 A 18 W — 360 JW — 720 UW - 3600 V^ 2 ) 

^155 = 2 ^ 2 Ci55_i 2 + Ci55_ 2 2 + Ci 5 5_ 3 2 

|*} 156 = |4,i,2,r B ,i) 

= C156-1 (|02uu) + |0uw2) + |20uu) + |2uw0) - |u02u) - |u20u) + \uu02) + |im20)) 
+ C156-2 (|0u2u) - |2m0m> - \uQu2) + \u2u0)) 
+ C156-3 (\duuu) — \uduu) + \uudu) — \uuud)) 

CiBe-i = \ (2Ai6-3 J-6U-60W) 

C*156-2 = — 4t 2 

Ci56-3 - ^ (-4A 16 2 + 12A 16 J 

—9 J 2 + 288 1 2 + 24 ^4i 6 U — 36JU 

-36 C/ 2 + 240 Ai 6 W - 360 J W - 720 U W - 3600 W 2 ) 

^156 = 2 ^ 2 Ci56-l 2 + Ci56-2 2 + Cl56-3 2 

|*) 1B7 = |4,i,2,r 5 , 2 ) 

= C157-1 (|02im) - |20uw) - |uu02) + \uu2Q)) 

+ Ci57_ 2 (|0tt2u) - |0uw2) + |2u0u) - |2um0) - |u02u) + |u0u2) - \u2Qu) + |u2u0}) 
+ Ci57_3 (\duuu) + \uduu) — \uudu) — \uuud)) 

<?157-1 = 4t 2 

C 157 -2 = ^(7+^-2^ + 2^008(63)) 

C157-3 - ^ (-A n 2 + 6 A n J 

—9 J 2 + 288 1 2 + 12 An U-36JU 

-36 U 2 + 120 An W - 360 J W - 720 U W - 3600 W 2 ) 

= 2 \J Ci57_i 2 + 2Ci57_2 2 + Ci57_3 2 
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|*)iB8 - |4,i,2,r 5 , 2 ) 

= Ci58_i (\02uu) - \20uu) - \uu02) + \uu20}) 

+ C158-2 (|0u2u) - |0uw2) + \2u0u) - \2uu0) - \uQ2u) + \u0u2) - \u20u) + |u2itO)) 
+ Ci58-3 (\duuu) + \uduu) — \uudu) — \uuud)) 

CiBs-i = 4t 2 

C 158 -2 = y (2^is-3J-6C/-60VK) 

Ci58- 3 = (-4A 18 2 + 12A 18 J 

-9 J 2 + 288 t 2 + 24 Ais U — 36JU 

-36 C/ 2 + 240 A 18 W — 360 JW — 720 UW - 3600 W 2 ) 



^158 — 2 ^ Ci58-i 2 + 2 Ci58-2 2 + Ci 58 -3 2 



|*) 159 = |4,i,2,r B , 2 ) 

= C159-1 (|02uu) - |20uw) - \uuQ2) + \uu20)) 

+ Ci 59 _2 (|0u2u) - |0uw2) + \2u0u) - \2uu0) - \u02u) + \u0u2) - \u20u) + \u2u0)) 
+ C159-3 (\duuu) + \uduu) — \uudu) — \uuud)) 

C159-1 = 4i 2 

C 159 - 2 = ^ {2A W --ZJ-§U -mw) 

C159-3 - ^ (-4A 16 2 + 12A 16 J 

—9 J 2 + 288 1 2 + 24 A w U — 36JU 

-36 U 2 + 240 A w W - 360 J W - 720 U W - 3600 VK 2 ) 



^159 — 2 \J Ci59_i 2 + 2 Ci5 9 _2 2 + C 



159-3 



|*)iao = |4,l,2,r B , 3 > 

= Cieo-i (|02im) - |0u2u) + |20um) - |2u0u) + |u0u2) + |u2u0) - |uw02) - \uu2Q)) 

+ Cieo-2 (\0uu2) - |2mm0) - \u02u) + \u20u)) 

+ C16O-3 (\duuu) — \uduu) — \uudu) + \uuud)) 

Cieo-i - I (j + U- 2W + 2^ cos(6 3 )) 

C160-2 = — 4t 2 

Cieo-3 - ^ (An 2 -6 A n J 

+ 9 J 2 - 288i 2 - 12^ii J7 + 36 JU 

+ 36 1/ 2 - 120 An W + 360 JW + 720 UW + 3600 VP 2 ) 

^160 = 2 \J 2 Ci6o-i 2 + C160-2 2 + C160-3 2 
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Schumann: Extended Hubbard model on small clusters 



|*)i6i = |4,i,2,r 5 , 3 > 

= Ciei-i (\02uu) - \0u2u) + \20uu) - \2u0u) + \uQu2) + \u2uQ) - \uu02) - \uu20)) 

+ Ciei-2 (|0m*2) - |2ww0) - \u02u) + \u20u)) 

+ C16I-3 (\duuu) — \uduu) — \uudu) + \uuud)) 



|*)i62 = |4,i,2,r B , 3 > 

= C162-1 (|02im) - \0u2u) + \20uu) - \2uQu) + \u0u2) + \u2u0) - \uuQ2) - \uu2Q)) 

+ C162-2 (\0uu2) - |2ww0) - \uQ2u) + \u20u)) 

+ C162-3 (\duuu) — \uduu) — \uudu) + \uuud)) 

C162-1 = y (2A W -3J-6U -60W) 

C162-2 = — 4t 2 

Ci 62 -3 = ^ (4^16 2 -12^1 6 J 




+ 9 J 2 - 288t 2 -24 A 18 U + 36 J U 



+ 36U 2 - 240 Ai 8 W + 360 JW + 720 UW + 3600 W 2 ) 




+ 9J 2 - 288t 2 -24Ai 6 J7 + 36J[7 

+ 36 U 2 - 240 Ai 6 W + 360 JW + 720 UW + 3600 V^ 2 ) 




Unnormalized eigenvectors for N e = 4 and m s = 2. 



|*>163 



|4,2,6,r 2 ) 
1 (|mt««)) 



Unnormalized eigenvectors for N e = 5 and m s = — |. 



I*) 164 



= - (\2ddd) - \d2dd) + \dd2d) - \ddd2)) 




|*}l65 



I r 3 15 

= l 5 '~2'T' M ' 
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I*)i66 = |5,-|^,r 5l2 ) 

= i (\2ddd) - \d2dd) - \dd2d) + \ddd2)) 



I*)i67 = |5,-|,^,r B ,3) 

= ^ (\2ddd) + \d2dd) - \dd2d) - \ddd2)) 



Unnormalized eigenvectors for N e = 5 and m s 



1 



— ^ (|022cf) + |02d2) + |0d22) + |202d) + |20rf2) + |220d) 
2 v 3 

+ |22d0) + |2d02) + |2d20) + |d022) + |rf202) + |d220» 



I*)i69 - |5,-1 ^,r 2 ) 



{\2ddu) + \2dud) + \2udd) - \d2du) - \d2ud) + \dd2u) 



2V3 

-\ddu2) + \du2d) - \dud2) - \u2dd) + \ud2d) - \udd2)) 

l*)i7o - |5,-^,r 3l i) 

= C170-1 (|022d) + |0d22) + |20rf2> + |220d) + |22d0) + |2d02) + |d022) + |d220}) 
+ C170-2 (|02d2) + |202rf) + |2d20) + |d202)) 

+ C170-3 {\2ddu) - \2udd) - \d2ud) - \dd2u) + \ddu2) + \du2d) + \u2dd) - \udd2)) 



Cl70-1 — 



2y/2 

Cl70-2 



V2 

C 170 _3 = jj= (\fT 7 + J + 2t + U-2W^ 



Nno = 2 \J 2 C170-1 2 + Ci 7 o- 2 2 + 2 C 



170-3 2 



l*>m = |5,-^,r 3 ,i) 

= Cm-i (|022d) + |0d22) + |20d2) + |220d) + |22d0) + |2d02) + |d022) + |d220}) 

+ C171-2 (|02d2) + |202d) + |2d20> + |d202)) 

+ C m _3 (|2ddu) - \2udd) - \d2ud) - \dd2u) + \ddu2) + \du2d) + \u2dd) - |udd2)) 
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Schumann: Extended Hubbard model on small clusters 



Cm-l = 
C171-2 = 
Ctfl-3 = 

N m = 

I*) 172 = 

+ 
+ 

C172-1 = 
C172-2 = 
C172-3 = 

^172 = 

I*) 173 = 

+ 
+ 

Cl73_l = 
Cl73-2 = 
Cl73-3 = 
^173 = 

I*) 174 = 

+ 
+ 



-t 

2 V2 
t 

71 



Ay/2 



J -2t-U + 2W) 



2y2Ci7i_i + Ci7i_2 +2Ci7i_3 Z 



' 5 ' ~2' 4' r3,2 ^ 

C172-1 (|022d) - |0d22) + |20d2) - |220d) - |22d0) + |2d02) - |d022) + |d220» 
Ci72~ 2 (\2ddu) + \2udd) - \d2ud) + \dd2u) - \ddu2) + \du2d) - \u2dd) - \udd2)) 
C172-3 {\2dud) - \d2du) - \dud2) + \ud2d)) 



+ J + 2t + U -2W 



-(#) 

2 

4V6 V 

— ^7= (y^+J + 2t + U-2w) 
2V6 V / 

2 ^ 2 C172-1 2 + 2 C172-2 2 + C; 



) 



172-3 



|5,--,-j-,r 3 , 2 ) 

C173-1 (|022d) - |0rf22) + |20d2> - |220rf> - |22rf0> + |2d02) - |d022) + |d220)) 
Ci73_ 2 (\2ddu) + \2udd) - \d2ud) + \dd2u) - \ddu2) + \du2d) - \u2dd) - \udd2)) 
C173-3 {\2dud) - \d2du) - \dud2) + \ud2d)) 



— ^= ( y/M - J -2t-U + 2W) 
4V6 V / 

— \= ( J~A~-, - J -2t-U + 2W) 
2V6 V v ) 

2 \j 2 Cn3-1 2 + 2 Ci73_2 2 + C; 



173-3 



|5,-^,r 4 ,i) 

C174-1 (|022d) + |02rf2) + |202d) + |20d2) - |2rf02) - |2d20) - |d202) - |d220)) 
C174-2 (|0d22) - |220rf) - |22d0> + |d022)) 

C174-3 (\2ddu) - \2dud) + \d2du) - \d2ud) - \du2d) - \dud2) + \ud2d) + \udd2)) 
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C174-1 
C174-2 
C174-3 



t 



2V2 
-t 

-1 
18 V2 



(j + U + U -2W + 2 y<% cos(6 5 )) 
( J - 12 i + U - 2 W + 2 cos(6 5 ) 



(A 22 2 - 3 A 22 f - 36 t 2 + 12 ^22 U - 18 1 [7 
+ 36 U 2 + 96 A 22 W - UAtW + 576U W + 230AW 2 ) 



N174 — 2 ^ 2 Ci74_i 2 + Ci74_2 2 + 2 C174- 



I*) 



175 — 



|5,-^,r 4 ,i) 



= C175-1 (|022d) + |02d2) + |202d) + |20d2) - |2d02) - |2d20) - |d202) - |d220}) 
+ C175-2 (|0d22) - |220rf) - |22d0) + |d022)) 

+ C175-3 (\2ddu) - |2dwd) + \d2du) - \d2ud) - \du2d) - \dud2) + \ud2d) + \udd2)) 



C175-1 
C175-2 
C175-3 



2^ 
-i 

3V2 
-1 

18 V2 



( J + 4f + {/ - 2 W - cos(e 5 ) - V3 sin(9 5 )) 
(j -I2t + U -2W - cos(9 5 ) - V3 a/Aj sin(6 5 )) 

(^24 2 - 3 Jt- 45 1 2 + 12 JC/ + 33 tU- 24 [7 2 + 96 JT4^ + 294 tW - 504 [7W 



^175 



-2496 VK 2 + 3 V^M cos(e 5 ) - 12 y/%U cos(9 5 ) - 96 ^/%W cos(6 5 ) 
+ 3V3 V^t sin(6 5 ) - 12 V3 ^/%U sin(6 5 ) - 96 V3 sin(6 5 )) 

2 ^/ 2 Ci75_i 2 + Ci75_2 2 + 2 Ci75_3 2 



I*)i76 = \5,~, ^,r 4jl ) 

= C176-1 (|022d) + |02rf2) + |202d) + |20d2) - |2rf02) - |2d20) - |d202) - |d220}) 

+ C176-2 (|0d22) - |220d) - |22d0> + |d022)) 

+ Ci 76 _3 (|2ddu) - \2dud) + \d2du) - \d2ud) - \du2d) - \dud2) + \ud2d) + |udd2)) 



Cl76-1 
Cl76-2 
Cl76-3 



t 



2V2 
-t 

3 V2 
-1 

18 V2 



\J + 4 1 + £/ -2W - cos(9 5 ) + V3 V^e sin(6 5 )^ 
(j -12t + U -2W - \[M cos(9 5 ) + V3 \[M sin(6 5 )) 
(A23 2 - 3 J t - 45 1 2 + 12 J £7 + 33 1 U - 24 U 2 + 96 J W + 294 1 W - 504 U W 



-2496 VK 2 + 3 V^6 i cos(6 5 ) - 12 ^A 6 U cos(9 5 ) - 96 y/A 6 W cos(6 5 ) 
-3V3 sin(6 5 ) + 12 V3 y^U sin(6 5 ) + 96 Vs VA5 W sin(6 5 )) 



^176 — 2 ^/ 2 Ci76-l 2 + Ci76-2 2 + 2 Ci 76- 
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Schumann: Extended Hubbard model on small clusters 




|*) 177 = 

+ 
+ 



C*177-l — 

Cl77-2 = 
Cl77-3 = 

N 177 = 

1 10178 = 

+ 
+ 



Cl78-1 — 



Cl78-2 = 
Cl78-3 = 
^178 = 

|*>179 = 

+ 
+ 



C177-1 (|022d) - |20rf2) - |2d02) + |rf220>) 

C177-2 (|02d2) + |0rf22) - |202d) - |220rf) + |22d0) + |2d20) - |d022) - |d202)) 
Ci 77 _3 (\2dud) - \2udd) + \d2du) + \dd2u) + \ddu2) - \dud2) - \u2dd) - \ud2d)) 

4 (-^22 2 + 3A 22 J+ 12 A 22 t -18 Jt 
18 

+ 45 1 2 - 9A 22 U + 18 JU + 54tU- 18U 2 - 102 A 22 W + 144 JW 
+ 612 1 W - 468 U W - 2592 W 2 ) 

I (j + 6t + U-2W-^Acos(e 5 )) 



C178-1 (|022d) - |20rf2) - |2d02) + |d220» 

C178-2 (|02d2) + |0rf22) - |202d) - |220rf> + |22d0) + |2d20) - |d022) - |d202)) 
Ci78_ 3 (\2dud) - \2udd) + \d2du) + \dd2u) + \ddu2) - \dud2) - \u2dd) - \ud2d)) 



— {-A 2 2 + 3 J 2 + 3Jt + 81t 2 -6JU -33W 
18 

+27 U 2 - 108 JW -2§UW + 4Q2UW + 2W8W 2 
-3 \/%J cos(9 5 ) - 12 t cos(6 5 ) + 9 \/A~qU cos(6 5 ) 
+102 yJ%W cos(6 5 ) - 3V3 \fA~6 J sin(6 5 ) - 12 V3 y/T Q t sin(6 5 ) 
+9V3 y/A^V sin(6 5 ) + 102 V3 y/A^W sin(6 5 )) 



C179-1 (|022d) - |20rf2) - |2d02) + |d220» 

C179-2 (|02d2) + |0rf22) - |202d> - |220rf> + |22rf0> + |2d20) - |d022) - |d202)) 
Ci79„ 3 (\2dud) - \2udd) + \d2du) + \dd2u) + \ddu2) - \dud2) - \u2dd) - \ud2d)) 
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Cl79-1 = 



Cl79-2 
Cl79-3 
^179 



I*) 



180 



+ 
+ 



Cl80-1 
Cl80-2 

Cl80-3 

iVi 8 o 



I*) 



181 



+ 

Cl81-1 = 
Cl81-2 = 



Cl81-3 

Nisi 



(-^2 3 2 + 3 J 2 + 3Ji + 81i 2 -6J[/-33iC/ 

18 

+27 U 2 - 108 JW — 294 tW + 492 UW + 2508 W 2 

-3 \/%J cos(9 5 ) - 12 t cos(6 5 ) + 9 y^t/ cos(6 5 ) 

+102 \/AW cos(6 5 ) + 3 V3 ■/ sin(9 5 ) + 12 V3 * sin(9 5 ) 

-9V3 y/AU sin(6 5 ) - 102 V3 sin(9 5 )) 

^ (2J+12i + 2C/-4W + cos(6 5 ) - V3 sin(6 5 )) 

^ (- J +12t-[/ + 2T^ + V^e cos(9 5 ) - V3 sin(6 5 )) 

2 \J Ci79_i 2 + 2 (Ci79_2 2 + Cl79-3 2 ) 
l 5 >-^^ r 4,3) 

Ciso-i (|022d) + |0d22) - |20d2) + |220rf) - |22d0) + |2d02) - |d022) - |d220}) 
Ciso-2 (|02d2) - |202rf) - |2d20) + |d202)) 

Ciso-3 (|2ddu) - |2udd) + |d2ud) - |dd2u) - |ddu2) + |du2rf) - \u2dd) + |udd2)) 



^A cos(9 5 )) 



- (J + 6t + U -2W 



— (-A 22 2 + 3 A 22 J +12 A 22 t- 18 J t 
18 



+ 45 1 2 -9A 22 U + 18 JU + 54tU- 18U 2 - 102 A 22 W + 144 JW 
+ 612 1 W - 468 U W - 2592 W 2 ) 



-t 
~6~ 



(-J +12t-U + 2W -2 y% cos(6 5 )) 
\J 2 Ciso-i 2 + Ciso-2 2 + 2 Ci 



'180-3 



l 5 ' _ 2'4' r4 ' 3 ^ 

Cisi-i (|022d) + |0d22) - |20d2) + |220rf) - |22rf0) + |2d02) - |d022) - |d220)) 
Cisi-2 (|02d2) - |202d) - |2d20) + |d202)) 

Cisi-3 (|2ddu) - |2wdd) + \d2ud) - |dd2u) - |ddw2) + |du2d) - \u2dd) + |udd2)) 
t 



- (2 J +12t + 2U -AW + y/A cos(6 5 ) + V3 y/A sin(9 5 )) 
(-^24 2 + 3 J 2 + 3Jf + 81t 2 -6J[/-33tC/ 

18 

+27 U 2 - 108 J 294 iVK + 492 C/VK + 250814^ 
-3 y/A J cos(9 5 ) - 12 y/~A~ 6 t cos(6 5 ) + 9 y/~AU cos(G 5 ) 
+102 y/AW cos(6 5 ) - 3\/3 sin(9 5 ) - 12\/3 sin(e 5 ) 

+9V3 y/A U sin(6 5 ) + 102 V3 y/AW sin(6 5 ) 



= y (-^+12i-?7 + 2VK+ yJ~A cos(e 5 ) + \/3v / ^sin(e5)) 



2 2 Cisi-i 2 + Cisi-2 2 + 2 C181-3 2 
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Schumann: Extended Hubbard model on small clusters 



I*)i82 = \5,~, |,r 4 , 3 > 

= C182-1 (|022d) + |0d22) - |20d2) + |220rf) - |22rf0) + |2d02) - |d022) - |d220}) 
+ C182-2 (|02d2) - |202rf) - |2d20) + |d202)) 

+ Ci 8 2-3 (\2ddu) - \2udd) + \d2ud) - \dd2u) - \ddu2) + \du2d) - \u2dd) + \udd2)) 
t 



Cl82-1 
Cl82-2 



= - (2J+12t + 2U-AW+ cos(9 5 ) - V3 sin(6 5 ) 



18 ^ 



Cl82-3 — 
^182 = 

I*) 183 = 

+ 

+ 

Cl83-1 = 
Cl83-2 = 

Cl83-3 = 

^183 = 

I*) 184 = 

+ 

+ 

Cl84-1 = 
Cl84-2 = 

Cl84-3 = 

^184 = 



- 3 J 2 + 3 J t + 81 i 2 - 6 J U - 33 tU 

+27U 2 - 108 JW -29UW + 492UW + 2508 W 2 

-3 \/%J cos(9 5 ) - 12 y/A^t cos(6 5 ) + 9 \/%U cos(6 5 ) 

+102 ifM W cos(6 5 ) + 3 V3 J sin(9 5 ) + 12 V3 y^t sin(e 5 ) 

-9V3 U sin(6 5 ) - 102V3^j%W sin(6 5 )) 

Y (- J +12t-U + 2W + cos(6 5 ) - V3 y<% sin(6 5 )) 

2 \J 2 Ci82-1 2 + Ci82-2 2 + 2 Ci82-3 2 

|5,-^,r 5>1 ) 

Ciss-i (|022d) - |0rf22) - |20d2) - |220rf) + |22d0) + |2d02) + |d022) - |d220}) 
Ci83- 2 (\2ddu) + \2udd) + \d2ud) + \dd2u) + \ddu2) + \du2d) + \u2dd) + \udd2)) 
Ci83_ 3 (|2dud) + \d2du) + \dud2) + \ud2d)) 



4V6 
-1 

2 V6 



^= (v^U + J -2t + U -2W 
(yCAi + J -2t + U -2W^ 



2 \/2 Ci83-1 2 + 2 Ci83-2 2 + Ci83-3 2 



'183- 

|5,-^,r 5>1 ) 

Ci84-i (|022d) - |0d22> - |20d2) - |220d) + |22d0) + |2d02) + |d022) - |d220)) 
Cis4-2 (|2ddu) + \2udd) + \d2ud) + \dd2u) + \ddu2) + \du2d) + \u2dd) + \udd2}) 
Ci84- 3 (|2dud) + \d2du) + \dud2) + \ud2d)) 



(yM- J + 2t-U + 2W^j 

— ^= (\fA~i- J + 2t-U + 2W) 
2V6 V / 

2 \J 2 Ci84-1 2 + 2 Ci84-2 2 + C*i84-3 2 
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185 = 



' 5 ' ~2~' T' 1 " 5 ' 1 ^ 



1 



2 V3 



(\2ddu) + \2dud) + \2udd) + \d2du) + \d2ud) + \dd2u) 
\ddu2) + \du2d) + \dud2) + \u2dd) + \ud2d) + \udd2)) 



186 



= Cise-i (|022d) - |02rf2) - |202d) + |20rf2> - |2rf02) + |2d20) + |d202) - |d220)) 
+ C186-2 (\2ddu) + \2dud) - \d2du) - \d2ud) - \du2d) + \dud2) - \ud2d) + \udd2)) 
+ C186-3 (\2udd) - \dd2u) + \ddu2) - \u2dd)) 



C186-1 
C186-2 

C186-3 

iVl86 



+ J-2t + U-2W 



) 



4V6 

— ^7= ( \fM + J-2t + U-2w) 
2V6 V / 

2 \J 2 Ci86-i 2 + 2 C186-2 2 + C : 



186-3 



Cl87-1 
Cl87-2 

Cl87-3 





15, -i. 


4,r Bl 2) 




Cl87-1 


(|022d) - 


+ 


Cl87-2 


(\2ddu) ■ 


+ 


Cl87-3 


(\2udd) - 




V 2 






2 






1 

4V6 






-1 
2 V6 





A4- J + 2t-U + 2W^j 
~A 4 - J + 2t-U + 2W^j 

Nis7 = 2 \J 2 Ci87-1 2 + 2 Ci87-2 2 + C'i87-3 2 



I*) 



188 



1 15 
2'T' 



r B ,2> 



2 V3 



(\2ddu) + \2dud) + \2udd) - \d2du) - \d2ud) - \dd2u) 



I*) 



189 



+ \ddu2) - \du2d) + \dud2) - \u2dd) - \ud2d) + \udd2)) 
= |5,-^, j,r 5)3 ) 

= Cigg-i (|02d2) - |0rf22) - |202d) + |220rf) - |22d0) + |2d20) + |d022) - |d202)) 
+ C189-2 (\2ddu) + \d2ud) - \du2d) - \udd2)) 

+ Ci 89 _3 (\2dud) + \2udd) + \d2du) - \dd2u) - \ddu2) - \dud2) + \u2dd) - \ud2d)) 
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Schumann: Extended Hubbard model on small clusters 



Ci89-1 
Ci89-2 
Ci89-3 
^189 
|*>190 



1 



4^6 



^ (yM + J -2t + U -2W^j 
+ J -2t + U -2W) 



2 \J 2 Ci89-1 2 + Ci89-2 2 + 2 Ci89-3 Z 



= l 5 '~2' 4' r5 ' 3 ^ 

= Cigo-i (|02d2) - |0rf22) - |202d) + |220rf> - |22rf0) + |2d20) + |d022) - |d202)) 

+ C190-2 {\2ddu) + \d2ud) - \du2d) - \udd2)) 

+ C190-3 {\2dud) + \2udd) + \d2du) - \dd2u) - \ddu2) - \dud2) + \u2dd) - \ud2d)) 



Cl90-1 
Cl90-2 
Cl90-3 
^190 
I*)l91 



( 



%- J + 2t-[/ + 2W 



2 V6 

— ^ f v / ^-^ + 2t-C/ + 2T^ N ) 
4\/6 V v / 

2 ^/ 2 Cigo-l 2 + Cigo-2 2 + 2 C; 

1 15 



190-3 



= 15,- 



2' 4 



r 5 , 3 > 



2 V3 



(|2ddu> + |2dwd) + |2udd) + |d2rfw) + |rf2ud) - \dd2u) 



192 



-|rfdw2) - \du2d) - \dud2) + \u2dd) - \ud2d) - \udd2)) 

Unnormalized eigenvectors for N e = 5 and m s = i. 

,13 
= 15,^^} 



1 



2 V3 



(|022u) + |02u2) + |0m22) + |202u) + |20u2) + |220u) 



I*) 



193 



+ |22u0) + |2u02) + |2u20) + |u022) + |«202) + |u220}) 

|5,i^,r a > 



1 



2 V3 



(\2duu) + \2udu) + \2uud) - \d2uu) + \du2u) - \duu2) 



194 — 



— \u2du) — \u2ud) + \ud2u) — \udu2) + \uu2d) — \uud2)) 

|5,^,r 3ll ) 



= C l94 -i (|022m) + |0u22) + |20u2) + |220u) + |22u0) + |2u02) + |u022) + |u220)) 
+ C194-2 (|02it2) + |202u) + |2w20) + |u202» 

+ Ci94_3 (\2duu) — \2uud) — \d2uu) + \duu2) + \u2du) — \ud2u) + \uu2d) — \uud2)) 
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C194-1 = 
C194-2 = 

C194-3 = 

N 194 = 

I*} 195 = 

+ 
+ 

C195-I = 
Cl95-2 = 
Cl95-3 = 
^195 = 

|*>196 = 

+ 
+ 

Cl96-1 = 
Cl96-2 = 
Cl96-3 = 
^196 = 

|*>197 = 

+ 
+ 



-t 

2 V2 
t 

71 



— ^7= (jA^+J + 2t + U-2w) 

2 \ 2 Ci94_i 2 + Ci94_2 2 + 2 Ci94_3 Z 



I 5 ) 2'i' r 3'!) 

Ci 95 _i (|022m) + |0u22) + |20w2) + |220u) + |22u0) + |2u02) + |u022) + |it220)) 
C195-2 (|02u2) + |202u) + |2w20) + |u202» 

Ci 95 _3 (\2duu) - \2uud) - \d2uu) + \duu2) + \u2du) - \ud2u) + \uu2d) - \uud2)) 



-t 

2V2 
t 

-1 



_ ( JA~i - J -2t-U + 2W) 
4V2 V v J 

2 \J 2 Ci95_i 2 + Ci95_2 2 + 2 C 



195-3 



Ci 96 _i (|022m) - |0u22) + |20u2) - |220u) - |22u0) + |2w02) - |u022) + |u220}) 
C196-2 (\2duu) + \2uud) — \d2uu) — \duu2) — \u2du) + \ud2u) + \uu2d) — \uud2)) 
Ci 96 _3 (\2udu) + \du2u) - \u2ud) - \udu2)) 



2 V6 



^ (V^7 + J + 2t + U -2W^j 
J + 2t + U -2W^j 



4\/6 



2 \J 2 Ci96-l 2 + 2 Ci96-2 2 + C; 



196-3 



|5,^, j,r 3i2 ) 

Ci 97 _i (|022w) - |0u22> + |20u2) - |220u) - |22u0) + |2u02) - |u022) + |u220)) 
C197-2 (\2duu) + \2uud) — \d2uu) — \duu2) — \u2du) + \ud2u) + \uu2d) — \uud2)) 
Ci 97 _3 (\2udu) + \du2u) - \u2ud) - \udu2}) 
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Schumann: Extended Hubbard model on small clusters 



Cl97-1 — 
Cl97-2 = 

Cl97-3 = 

^197 = 

|*>198 = 

+ 
+ 

Cl98-1 = 
Cl98-2 = 
Cl98-3 = 

^198 = 

|*>199 = 

+ 
+ 

Cl99_l = 
Cl99-2 = 
Cl99-3 = 



^199 
|*>200 



+ 
+ 



— ^= f - J -2t-U + 2w) 
— \= ( JTi - J -2t-U + 2w) 

2 \J 2 Ci97_i 2 + 2 Ci97_2 2 + Ci97_3 2 

.13 

' ' 2' 4' A ' 

Ci 98 _i (|022u) + |02u2> + |202u) + |20u2) - |2u02) - |2u20> - |u202) - |u220» 
C198-2 (|0u22) - |220u) - |22u0) + |u022)) 

C198-3 (\2udu) - \2uud) - \du2u) - \duu2) + \u2du) - \u2ud) + \ud2u) + \udu2)) 



t 



2V2 



(j + At + U -2W + 2 cos(e 5 )) 



3^2 
-1 

18 V2 



J- I2t+U -2W + 2WA e cos 



(©5)) 



(A 22 2 - 3 A 22 1 - 36 t 2 + 12 A 22 U - 18 1 U 



+ 36 U 2 + 96A 22 W- UUW+ 576 UW + 2304 W 2 ) 



2 \J 2 Ci98-l 2 + Ci98- 2 2 + 2 C 



198-3 



I 5 ! 2' 4' r4 ^) 

C199-1 (|022m) + |02u2) + |202u) + |20u2) - |2u02) - |2u20) - |u202) - |u220)) 
C199-2 (|0u22) - |220u) - |22u0) + |u022» 

Ci 99 _3 (\2udu) - \2uud) - \du2u) - \duu2) + \u2du) - \u2ud) + \ud2u) + \udu2}) 



2V2 



(j + At + U -2W - v / ^cos(e 5 )-\/3V^sin(e 5 )) 
l -j= (j -12t + U -2W - ^/A~q cos(9 5 ) - VH \[K sin(6 5 )) 



3-\/2 
-1 

18 V2 



(A 24 2 - 3 J t - 45 1 2 + 12 J U + 33 1 U - 24 U 2 + 96 J W + 294 1 W - 504 U W 

-2496 W 2 + 3 i/Tfit cos(e 5 ) - 12 V^t 7 cos(9 5 ) - 96 y/AW cos(6 5 ) 
+ 3V3 sin(6 5 ) - 12VSyj%U sin(6 5 ) - 96 VSy/AW sin(6 5 )) 



2 \/2 Cigg_i 2 + Cigg_2 2 + 2 Cigg_3 2 



|5,^,r 4il ) 

C200-1 (|022m) + |02u2) + |202u) + |20u2) - |2u02) - |2u20) - |u202) - |u220)) 
C200-2 (|0u22) - |220u) - |22u0) + |u022» 

C200-3 (\2udu) - \2uud) - \du2u) - \duu2) + \u2du) - \u2ud) + \ud2u) + \udu2)) 
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C200-1 
6200-2 
C200-3 



t 



- (j + 4t + U-2W-M%cos(9 5 ) + V3M%sm(G 5 )) 
(j -12t + U -2W - cos(9 5 ) + VS x/% sin(6 5 )) 



18 V2 



(A 23 2 - 3 Jt- 45 1 2 + 12 JU + mtU - 24[/ 2 + 96 JW + 2MtW - 504UW 

-2496 W 2 + 3 M%t cos(6 5 ) - 12 M%U cos(9 5 ) - 96 M% W cos(6 5 ) 
-3V3 M%t sin(6 5 ) + 12 V3 V^e^ sin(6 5 ) + 96 V3 M%W sin(6 5 )) 



-^200 — 2\/2C2oo-i +6200-2 +2(7200-3^ 



I*) 



201 — 



|5,^, j,r 4)2 ) 



= C201-1 (|022u) - |20u2> - |2w02) + |u220)) 

+ C201-2 (|02u2) + |0u22) - |202u) - |220u) + |22u0) + |2w20) - |u022) - |u202» 
+ C201-3 {\2duu) - \2udu) + \d2uu) + \du2u) - \u2ud) + \udu2) - \uu2d) - \uud2}) 



C201- 



1 = 



C201-2 
C201-3 



(-A22 2 + 3A 22 J +12 A 22 t- 18 J t 

+ 45 1 2 - 9A 22 U + 18 JU + 54tU- 18U 2 - 102 A 22 W + 144 JW 
+ 612 1 W - 468 U W - 2592 W 2 ) 

I (j + 6t +U- 2 W- VAcos(9 5 )) 
(- J +12t-[/ + 2T^-2 cos(6 5 )) 



f 
6 

2M C. 



201-1 



+ 2 (C : 



201-2 



201-3 



l*) 2 02 - |5,-,-,r 4 , 2 ) 

= C202-1 (|022m) - |20u2) - |2u02) + |u220)) 

+ C202-2 (|02u2) + |0u22) - |202u) - |220u) + |22u0) + |2u20) - |u022) - |u202)) 

+ C202-3 {\2duu) - \2udu) + \d2uu) + \du2u) - \u2ud) + \udu2) - \uu2d) - \uud2}) 



a 



202-1 



C202-2 
C202-3 

N 2 o 2 



— (-A 24 2 + 3 J 2 + 3 J t + 81 1 2 - 6 J U - 33 1 U + 27 U 2 - 108 J W - 294 1 W + 492 U W 



18 



+ 2508 W 2 - 3 M% J cos(9 5 ) - 12 M%t cos(6 5 ) + 9 M%U cos(6 5 ) + 102 M%,W cos(6 5 ) 
—3 \[2>\j~A§J sin(6 5 ) - 12 V3 M%t sin(6 5 ) + 9\/3 M%U sin(6 5 ) + 102 V3 M%W sin(6 5 )) 



- [2J+l2t + 2U-AW + M% cos(9 5 ) + V3 M% sin(6 5 
6 V 



- (- J + 12 1 - U + 2 W + M% cos(9 5 ) + V3 sin(6 5 )) 
2 ^/ C202-1 2 + 2 (C202-2 2 + C202-3 2 ) 
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Schumann: Extended Hubbard model on small clusters 



1^)203 = 

+ 
+ 



C203-I - 

C203-2 = 
C203-3 = 
-^203 = 

I*) 204 = 

+ 
+ 

C*204-l = 
C204-2 = 
C*204-3 = 

N 2 m = 



|*>205 = 

+ 
+ 



|5, -, ^,r 4: 2) 

C203-1 (|022m) - |20u2) - |2m02) + |u220)) 

C203-2 (|02u2) + |0u22) - |202u) - |220u) + |22u0) + |2u20) - |u022) - |w202)) 
C203-3 {\2duu) — \2udu) + \d2uu) + \du2u) — \u2ud) + \udu2) — \uu2d) — \uud2)) 



\ {-A 2 z 2 + 3 J 2 + 3 J t + 81 t 2 - 6 J U - 33 1 U + 27 U 2 - 108 J W - 294 1 W + 492 U W 
18 

+ 2508 W 2 - 3 \/~A§J cos(6 5 ) - 12 \/A~ e t cos(9 5 ) + 9 y/AU cos(6 5 ) + 102 y/A^W cos(6 5 ) 
+ 3V3 y%J sin(6 5 ) + 12 V3 \/I^t sin(6 5 ) - 9 V3 y^U sin(6 5 ) - 102 V3 y^W sin(9 5 )) 

- (2J +12t + 2U -4W + JA~ & cos(9 5 ) - V3 JA~ sin(6 5 
6 V 

l - (- J + 12 1 - U + 2 W + y^T 6 cos(e 5 ) - V3 y% sin(6 5 )) 



6 

2 \J C203-1 2 + 2 (C203-2 2 + C203-3 2 ) 



l 5 >2'4' r4 . 3 ) 

C204-1 (|022u> + |0u22) - |20u2) + |220u) - |22u0) + |2u02> - |u022) - |u220» 
C204-2 (|02u2) - |202u) - |2w20) + |u202)) 

C204-3 {\2duu) — \2uud) + \d2uu) — \duu2) — \u2du) — \ud2u) + \uu2d) + \uud2)) 



-t 



2V2 
t 



(j + U + U -2W + 2 y% cos(6 5 )) 
( J - 12 1 + U - 2 W + 2 yfU cos(9 5 )) 



3 V2 

— ^= (A 22 2 - 3 A 22 1 - 36 1 2 + 12 ^22 U - 18 1 U 
18 V2 v 

+ 36 U 2 + 96A 22 W- UAtW + 576 UW + 2304 W 2 ) 



2 \J 2 C204-1 2 + C204-2 2 + 2C\ 



204-3 



|5,^,r 4i3 ) 

C205-1 (|022m) + |0u22) - |20u2) + |220u) - |22u0) + |2u02) - |u022) - |u220)) 
C205-2 (|02it2) - |202u) - |2w20) + |u202» 

C205-3 {\2duu) — \2uud) + \d2uu) — \duu2) — \u2du) — \ud2u) + \uu2d) + \uud2}) 
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C205-1 
C205-2 
C205-3 



-t 

2 V2 
t 

-1 
18 V2 



(^J + it + U - 2W - \Z~Aq cos(6 5 ) - V3 \/~A% sin(0 5 )^ 
(j -12t + U -2W - cos(9 5 ) - V3 \/~A % sin(6 5 )) 

{A 2 2 - 3 Jt- Aht 2 + 12 JU + 33W- 24[/ 2 + 96 JW + 2MtW - 50AUW 

-2496 W 2 + 3 y/%t cos(6 5 ) - 12 y/%U cos(9 5 ) - 96 W cos(6 5 ) 
+ 3V3 sin(6 5 ) - 12 V3 sin(6 5 ) - 96 V3 a/Aj W sin(6 5 )) 



iV 2 05 



2 \ 2 C205-I + C205-2 + 2 C205-3 



l*) 206 - |5,^,|,r 4 , 3 ) 

= C206-1 (|022u) + |0u22) - |20u2) + |220u) - |22u0) + |2u02) - |u022) - |u220)) 

+ C206-2 (|02u2) - |202u> - |2w20) + |u202)) 

+ C206-3 (\2duu) - \2uud) + \d2uu) - \duu2) - \u2du) - \ud2u) + \uu2d) + \uud2}) 



C206-1 
C206-2 
C206-3 



2V2 
t 

3 V2 
-1 

18 V2 



(j + At + U -2W - yf~A % cos(9 5 ) + V3 \fA~a sin(6 5 )) 
(j -I2t + U -2W - cos(9 5 ) + VS \fM sin(6 5 )) 

(A23 2 - 3 Jt- 45 1 2 + 12 JU + 33 tU- 24 U 2 + 96 JW + 294 tW - 504 UW 

-2496 W 2 + 3 x/%t cos(6 5 ) - 12 a/I^C/ cos(9 5 ) - 96 \/AW cos(6 5 ) 
-3V3 /^t sin(6 5 ) + 12 V3 y/A^U sin(6 5 ) + 96 V3 W sin(6 5 )) 



^206 



2 A/2 C206-I + C206-2 + 2 C206-3^ 



l*>207 = |5,|,|,r 5;1 > 

= C207-1 (|022u) - |0u22) - |20w2) - |220u) + |22u0) + |2u02) + |u022) - |u220)) 

+ C207-2 {\2duu) + \2uud) + \d2uu) + \duu2) + \u2du) + \ud2u) + \uu2d) + \uud2}) 

+ C207-3 {\2udu) + \du2u) + \u2ud) + \udu2)) 



C207-1 
C207-2 

C207-3 

N207 



X ->= (yM + J -2t + U -2W^j 



4V6 



— \= ( \fA~i + J -2t + U -2W) 
2V6 V / 

2 \J 2 C207-1 2 + 2 C207-2 2 + C207-3 
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Schumann: Extended Hubbard model on small clusters 



l*>208 = |5,|,|,r 5! i> 

= C208-1 (|022u) - |0u22) - |20w2) - |220u) + |22u0) + |2u02) + |u022> - |u220)) 

+ C208-2 (\2duu) + \2uud) + \d2uu) + \duu2) + \u2du) + \ud2u) + \uu2d) + \uud2)) 

+ C208-3 (\2udu) + \du2u) + \u2ud) + \udu2)) 



C208-1 
C208-2 
C208-3 
N208 



M4- J + 2t-U + 2W 



) 



4V6 

= — ^= (JM- J + 2t-U + 2W 
2V6 V 

= 2 \J 2 C208-1 2 + 2 C208-2 2 + C*208-3 Z 



1 1 ^ 

l*>209 - |5,-, T ,r M > 



— !-= (\2duu) + \2udu) + \2uud) + \d2uu) + \du2u) + \duu2) 
2 v 3 

+ \u2du) + \u2ud) + \ud2u) + \udu2) + \uu2d) + \uud2)) 



l*>2io = |5,^|,r B , 2 > 

= C210-1 (|022u) - |02u2) - |202u) + |20u2) - |2u02) + |2w20) + |u202) - |u220)) 

+ C210-2 (\2duu) - \d2uu) - \uu2d) + \uud2)) 

+ C210-3 (\2udu) + \2uud) - \du2u) + \duu2) - \u2du) - \u2ud) - \ud2u) + \udu2}) 



C210-1 
C210-2 

C210-3 

^210 



v/fi (' 



7 A 4 + J -2t + U -2W 



) 



2 V6 

—^7= f V^I + J- 2i + £7-2W > ) 



2\/2C2io-i + C210-2 +2C < 2io-3 / 



l*>2n = |5,|,|,r B>2 ) 

= C211-1 (|022w) - |02u2) - |202u) + |20u2) - |2u02) + |2u20) + |u202) - |u220)) 

+ C211-2 (\2duu) - \d2uu) - \uu2d) + \uud2)) 

+ C211-3 (\2udu) + \2uud) - \du2u) + \duu2) - \u2du) - \u2ud) - \ud2u) + \udu2)) 
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C211-1 — 
C211-2 = 

C211-3 = 



1* 



2 

^= (v^4l- J + 2t-C/ + 2T4^) 
(yM- J + 2t-U + 2W^j 
+ C211-2 +2C211-3 



I*>2i2 = |5,^,r 5 ,2> 



— -= (\2duu) + \2udu) + \2uud) - \d2uu) - \du2u) + \duu2) 
2 V 3 

-\u2du) - \u2ud) - \ud2u) + \udu2) - \uu2d) + \uud2)) 



|*}2i3 = |5,^^r 5 , 3 



2' 4' 

= C213-1 (|02u2) - |0u22) - |202u) + |220u) - |22u0) + |2u20) + |u022) - |u202)) 
+ C213-2 (\2duu) + \2udu) + \d2uu) - \du2u) + \u2ud) - \udu2) - \uu2d) - \uud2)) 
+ C213-3 (\2uud) - \duu2) + \u2du) - \ud2u)) 



C213-1 — 
C213-2 = 

C213-3 = 



-AM 

2 

— \= ( \[A± + J-2t + U-2w) 

— ^7= ( \[A± + J-2t + U-2w) 
2VQ V ' 



N213 — 2 \J 2 C213-1 2 + 2 C213-2 2 + C*2i3- 



2 

3 



|*>2i4 = |5,^?,r B , 3 



2' 4' 

= C214-1 (|02u2> - |0u22) - |202u) + |220u) - |22u0) + |2u20) + |u022) - |u202» 
+ C214-2 (\2duu) + \2udu) + \d2uu) — \du2u) + \u2ud) — \udu2) — \uu2d) — \uud2)) 
+ C214-3 (\2uud) - \duu2) + \u2du) - \ud2u)) 



C214-1 
C214-2 



2 

(yM- J + 2t-U + 2W^j 
C214-3 = 0=(\/A^-J + 2t-U + 2W^) 
N214 = 2 \J 2 C214-1 2 + 2 C214-2 2 + C*214-3 2 
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Schumann: Extended Hubbard model on small clusters 



I*>2i5 = |5,^,r Bl3 > 

= — !-= (\2duu) + \2udu) + \2uud) + \d2uu) - \du2u) - \duu2) 
2 v 3 

+ \u2du) + \u2ud) - \ud2u) - \udu2) - \uu2d) - \uud2)) 
Unnormalized eigenvectors for N e = 5 and m s = |. 

I*>2i6 = |5,|f,r 2 ) 

= - {\2uuu) — \u2uu) + \uu2u) — \uuu2)) 

I*>2i7 = |5,|,^,r M ) 

= 1 (\2uuu) + \u2uu) + \uu2u) + \uuu2)) 

I*>2i8 = |5,|^,r 5 , 2 ) 

= 1 (\2uuu) — \u2uu) — \uu2u) + \uuu2)) 

I*>2i9 = |5,^,r 5i3 > 

= 1 (\2uuu) + \u2uu) — \uu2u) — \uuu2)) 

Unnormalized eigenvectors for N e = 6 and m s = —1. 

|*) 22 o = |6,-l,2,r 4 ,i) 

= 1(| 2d2d) + 1 2dd2) + \d22d) + | d2d2) ) 

|*>22i = |6,-i,2,r 4 , 2 ) 

= 1 (|22dd) + \2d2d) - \d2d2) - \dd22)) 



|*) 2 22 = |6,-l,2,r 4 , 3 > 

= 1 (|22dd) - \2dd2) + \d22d) + |dd22)) 

|*) 223 = |6,-l,2,r 6 ,i) 

= 1 (|22dd) + \2dd2) - \d22d) + |dd22)) 



|*>2 24 - |6,-l,2,r 5j2 ) 

= 1 (\2d2d) - \2dd2) - \d22d) + \d2d2)) 

|*) 225 = |6,-l,2,r 5i3 ) 

= 1 (|22dd) - \2d2d) + \d2d2) - \dd22)) 
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Unnormalized eigenvectors for N e = 6 and m s = 0. 



226 



C226-1 
C226-2 

N226 

|*>227 



= 16,0,0,^) 

= C 226 -i (|0222) + |2022) + |2202) + |2220>) 
+ C 226 - 2 (\22du) - \22ud) + \2d2u) + \2du2) - \2u2d) - \2ud2) 
+ \d22u) + \d2u2) + \du22) - \u22d) - \u2d2) - \ud22)) 



= 2 \J C226-1 2 + 3 C226-2 2 



C227-1 
C227-2 

N227 

|*) 2 28 



= 16,0,0,^) 

= C227-1 (|0222) + |2022) + |2202) + |2220>) 

+ C227-2 (\22du) - \22ud) + \2d2u) + \2du2) - \2u2d) - \2ud2) 
+ \d22u) + \d2u2) + \du22) - \u22d) - \u2d2) - \ud22)) 

= Vst 

= — \= (JA~ & - J-U-U + 2w) 

= 2 \J C227-1 2 + 3 C227-2 2 

= |6,o,o,r 3 ,i) 

= C228-1 {\22du) - \22ud) + \2du2) - \2ud2) + \d22u) + \du22) - \u22d) - \ud22)) 

+ C228-2 {\2d2u) - \2u2d) + \d2u2) - \u2d2)) 



C228-1 
C228-2 
N228 
|*) 229 



2 V6 
1 

= 2 \J 2 C228-1 2 + C; 



230 



C230-1 
C230-2 

^230 



228-2 



= |6,o,o,r 3 , 2 ) 

= — !-= {\22du) - \22ud) - \2du2) + \2ud2) - \d22u) + \du22) + \u22d) - \ud22)) 
2 v 2 

= |6,o,o,r 4 ,i) 

= C230-1 (|0222) + |2022) - |2202) - |2220)) 

+ C230-2 {\22du) - \22ud) - \du22) + \ud22)) 

= -t 



X - (^M + J + U -2W) 
2 \J C230-1 2 + C230-2 2 
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Schumann: Extended Hubbard model on small clusters 



|*>23i = |6,o,o,r 4 ,i> 

= C231-1 (|0222) + |2022> - |2202> - |2220>) 

+ C231-2 (|22du) - \22ud) - \du22) + \ud22)) 

C231-1 = —t 

C 23 i-2 = \(-^fM + J + U-2W) 



N231 



2 \l C231-1 2 + C231-2 2 



|*) 232 = |6,0,2,r 4 ,i) 

= — !-= (\2d2u) + \2du2) + \2u2d) + \2ud2) + \d22u) + \d2u2) + \u22d) + |u2d2)) 
2 V2 

l*) 233 = |6,o,o,r 4 , 2 ) 

= C233-1 (|0222) - |2022) - |2202> + |2220>) 

+ C 23 3-2 (|2du2) - |2wd2) - |d22u) + \u22d)) 

C 2 33-l = —t 

C233-2 = \ (^M + J + U -2W) 

N233 = 2 \J C 2 33_i 2 + C233-2 2 

l*) 234 - |6,o,o,r 4 , 2 ) 

= C234-1 (|0222) - |2022) - |2202) + |2220)) 

+ C234-2 (|2du2) - |2wd2) - |d22u) + \u22d)) 



C234-1 — 

C234-2 = 



1 A x + J + U - 2 W 



N234 — 2 \J C234-1 2 + C; 
|*) 235 



234-2 



I*) 



236 



|6,0,2,r 4 , 2 

1 



2 V2 



(|22du) + \22ud) + \2d2u) + \2u2d) - \d2u2) - \du22) - \u2d2) - \ud22)) 



= |6,o,o,r 4 , 3 ) 

= C 236 -i (|0222) - |2022) + |2202) - |2220» 
+ C 236 -2 (|2d2u) - |2u2d) - |d2u2) + |u2d2)) 



C 2 36-l 
C236-2 

^236 



-t 



X - (^M + J + U -2W) 
2 \J C236-1 2 + C% 



^236-2 
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|*>237 = |6,o,o,r 4 , 3 } 

= C237-1 (|0222) - |2022> + |2202> - |2220)) 

+ C237-2 (|2d2u) - \2u2d) - \d2u2) + \u2d2)) 

C237-1 = —t 

C 237 -2 = ^(-^M + J + U-2W^ 

N237 = 2 \J C237-1 2 + C237-2 2 

|*) 23 8 - |6,o,2,r 4 , 3 > 

= — !-= {\22du) + \22ud) - \2du2) - \2ud2) + |rf22u) + \du22) + \u22d) + |ud22)) 
2 v 2 

|*) 23 9 = |6,0,2,r 6 ,i) 

= — ^ {\22du) + \22ud) + \2du2) + \2ud2) - \d22u) + \du22) - \u22d) + \ud22)) 
2 v 2 

|*) 24 o - |6,0,2,r 6 ,2) 

= -^-= (\2d2u) - \2du2) + \2u2d) - \2ud2) - \d22u) + \d2u2) - \u22d) + \u2d2}) 
2 v 2 

|*) 24 i = |6,0,2,r 5 , 3 ) 

= -^-j= (\22du) + \22ud) - \2d2u) - \2u2d) + \d2u2) - \du22) + \u2d2) - \ud22)) 

Unnormalized eigenvectors for N e = 6 and m s = 1. 

|*) 242 = |6,i,2,r 4 ,i> 

= ]^{\2u2u) + \2uu2) + \u22u) + \u2u2)) 

|*) 243 = |6,l,2,r 4 , 2 > 

= i (\22uu) + \2u2u) - \u2u2) - \uu22)) 

|*) 244 = |6,l,2,r 4 , 3 ) 

= i (|22uu) - |2um2) + \u22u) + |uu22)) 

|*) 245 = |6,l,2,r 5il ) 

= i (|22mt) + |2uw2) - \u22u) + |uu22)) 

|*) 246 = |6,l,2,r 6 , 2 > 

= i (|2u2u) - \2uu2) - |u22u) + |u2u2» 

|*) 2 47 = |6,l,2,r 6 , 3 ) 

= i (|22uu) - \2u2u) + \u2u2) - \uu22)) 
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Unnormalized eigenvectors for N e = 7 and m s = — |. 

I*>248 - l 7 '4'I' ri) 

= i (|222d) + |22rf2) + |2d22) + |d222)) 

l*>249 = |7,-^,r 4 ,i) 

= i (|222d) + |22rf2) - |2d22) - |d222» 

l*>25o = |7,-^,r 4 , 2 ) 

= i (|222d) - |22d2) - |2d22) + |rf222}) 

l*>25i = |7,-^,r 4 ,3) 

= i (|222d) - |22rf2) + |2d22) - |d222» 
Unnormalized eigenvectors for N e = 7 and m s = i. 

I*) 252 = I7,^,rx) 

= i (|222u) + |22u2) + |2u22) + |u222» 

l*>253 = |7,^,r 4 ,i) 

= i (|222u) + |22u2) - |2u22) - |u222)) 

l*> 25 4 - |7,^,|,r 4 , 2 ) 

= i (|222u) - |22u2) - |2u22) + |u222» 

l*>255 = |7,|,|,r 4 , 3 > 

= i (|222u> - |22u2) + |2u22) - |u222)) 

Unnormalized eigenvectors for N e = 8 and m s = 0. 

I*) 256 = |8,o,o,ri) 

= 1 (|2222>) 
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Fig. 1 The canonical spectrum (scaled to primed values) of the triangle in dependence on U' (left panel), 
J' (middle panel) , and W' (right panel) for two, three, and four respectively electrons on the cluster. The 
values of the other parameters are indicated within the figures. 



N=2, UA=4, J/tsO N=3, U/t=4, W/teO N=4, Urt=4, W/t=0 




Fig. 2 The crossings of the canonical groundstates in dependence of the exchange parameter J for the 
triangular cluster. The insets sketch the J-dependencies for the complete canonical spectra. For details 
see explicit formula in appendix A. 



104 



Schumann: Extended Hubbard model on small clusters 
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-1 -0.5 0.5 1 



Fig. 4 Left panel: N(fj,) of the triangle for different values of the on-site correlation U. Right panel: 
The electron occupation for the related fi-U plane scaled to primed variables. The meaning of the colours 
is according to the palette given in Fig. 5 




-2 -0.5 0.5 1 



Fig. 5 Left panel: The electron occupation and S z eigenvalues of the groundstate of the triangular 
cluster for the complete jj,-h plane (scaled to primed variables) . Right panel: The palette used here and in 
Figs. 4, 6, and 7 to differ the groundstates with respect of their electron occupation and S z eigenvalues. 
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Schumann: Extended Hubbard model on small clusters 




Fig. 7 Left panel: N(/i) of the triangle for different values of the nearest-neighbour Coulomb interaction 
W for U/t = 4. Right panel: The electron occupation and S z eigenvalues of the groundstate for the 
related [i-W plane scaled to primed variables. The meaning of the colours is the same as in Fig. 5. 
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Fig. 8 The lowest part of the grand canonical spectrum for the triangle with and without a magnetic 
field (left panel) and with ferromagnetic and antiferromagnetic exchange taken into account (right panel). 
The energylevels are labeled by the magnetic quantum numbers and the irreducible representation of the 
point group. 
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Schumann: Extended Hubbard model on small clusters 
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Fig. 9 Left: The complete canonical spectrum for the tetrahedron (scaled to primed values) in depen- 
dence on U' for h/t — 0, W/t — and J/t = 0. Right: The groundstates for the complete h — £/-plane 
(scaled to primed values) for W/t = and J/t — 0. The indicated numbers, referring to the appendix 
B show which states are degenerated. The colours contain the information about the irr. representation 
of the point group and their intensity shows the value of S z as depicted in the palette in the left lower 
corner. Since the system is symmetrically with respect to a simultaneous sign change of magnetic field 
and spin projection we gave the numbers of the right half-plane only. 
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Fig. 10 Left: The complete canonical spectrum for the tetrahedron (scaled to primed values) in depen- 
dence on J' for h/t = 0, U/t = 4 and W/t = 0. Right: The groundstates for the complete J — £/-plane 
(scaled to primed values) for W/t = and h/t = 0. The numbers and colours have the same meaning as 
in Fig. 9. 
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Schumann: Extended Hubbard model on small clusters 




Fig. 11 The ground state levels for three electrons 
on the tetrahedron in dependence of the exchange 
parameter J for h/t = 0, U/t = 4 and W/t = 0. 



Fig. 12 The ground state levels for four electrons 
on the tetrahedron in dependence of the exchange 
parameter J for h/t = 0, U/t = 4 and W/t = 0. 
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Fig. 13 Left: The complete canonical spectrum for the tetrahedron (scaled to primed values) in depen- 
dence on W . The values of the other parameters are indicated in the figures. Right: The groundstates for 
the complete W — [/-plane (scaled to primed values) for J/t = and h/t — 0. The numbers and colours 
have the same meaning as in Fig. 9 
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Schumann: Extended Hubbard model on small clusters 





0.5 



-0.5 



Fig. 14 Left: N(fi) for the tetrahedron for different values of the on-site correlation U. Right: The 
electron occupation for the whole yAJ plane (scaled to primed values). The curves with U/t — —4, 
U/t — 0, and U/T = 4 in the left picture correspond to horizontal cuts across the right one for U' = —0.5, 
U' — 0, and U' = 0.5 respectively. 



h' 




Fig. 15 Left panel: The electron occupation iV and the 5 z -eigenvalue for the ground state in dependence 
on the chemical potential and applied magnetic field (both scaled to primed values) for the tetrahedron 
and U/t = 4, J/t = 0, and W/t = 0. Right panel. The palette used here and in Figs. 14,19, and 20 to 
differ the groundstates with respect of their electron occupation and Sz eigenvalues. 
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Fig. 16 The grand-canonical spectrum for the tetrahedron at n = — t. The other parameters are U/t — 4, 
J/t = 0, W/t = 0, and h/t = 0. 
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Fig. 17 The grand-canonical spectrum for the tetrahedron at = 2.4056 t. The other parameters are 
U/t = 4, J ft = 0, W/t = 0, and h/t = 0. 
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Schumann: Extended Hubbard model on small clusters 





J' 




Fig. 19 Left panel: N{fj) of the tetrahedron for the on-site correlation fixed to U = At for an antifer- 
romagnetic (J/t — 1) and a ferromagnetic exchange (J = —It). Right panel: The electron occupation 
iV and the ^-eigenvalue for the ground state in dependence on the chemical potential and exchange 
parameter J (both scaled to primed values) for the tetrahedron and U/t — 4, h/t = 0, and W/t — 0. The 
curves with J = —t and J = t in the left picture correspond to horizontal cuts across the right one for 
J' = —0.5 and J' = 0.5 respectively. The meaning of the colours is according to the palette given in Fig. 
15. 
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Fig. 20 Left panel: N(^) of the tetrahedron for the on-site correlation fixed to U/t = 4 for a repulsive 
(W/t = 1) and an attractive (W/t — —0.1) nn Coulomb interaction. Right panel: The electron occupation 
N and the S z -eigenvalue for the ground state in dependence on the chemical potential and nn-Coulomb 
interaction parameter W (both scaled to primed values) for the tetrahedron and U/t = 4, h/t = 0, and 
J/t = 0. The curves with W/t = 1 and W/t = — 1 in the left picture correspond to horizontal cuts across 
the right one for W' = 0.5 and W' = —0.5 respectively. The meaning of the colours is according to the 
palette given in Fig. 15. 



